Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,449,624 Δ1 bp 7.7% coding (426/822 nt) yfcO ← DUF2544 family putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,449,6240C. 7.7% 314.7 / 21.6 104coding (426/822 nt)yfcODUF2544 family putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (47/49);  new base . (4/4);  total (51/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

Only 100 of 111 total aligned reads displayed.

TACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAGGACAACTTGCTTTGTT  >  NC_000913/2449475‑2449773
                                                                                                                                                      |                                                                                                                                                     
tACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATAc                                                                                                                                                       <  1:1536635/150‑1 (MQ=255)
tACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATAc                                                                                                                                                       <  1:133829/150‑1 (MQ=255)
    ttGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGGGCCTGATA‑GTCAc                                                                                                                                                  <  1:1207757/150‑1 (MQ=255)
     tGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCAc                                                                                                                                                  <  1:1376728/150‑1 (MQ=255)
     tGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAAACGTGACGTGAC‑TATAAGGCCCCTGATACGTAAc                                                                                                                                                  >  2:334563/1‑150 (MQ=255)
     tGAAAGCGTTTTCTCAACCACACCGCCCGGCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCAc                                                                                                                                                  <  1:1370043/150‑1 (MQ=255)
      gAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACg                                                                                                                                                 >  1:1170038/1‑150 (MQ=255)
        aaGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGcc                                                                                                                                               >  1:1581997/1‑150 (MQ=255)
         aGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGccc                                                                                                                                              >  2:1402592/1‑150 (MQ=255)
          gCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGcccc                                                                                                                                             <  2:793768/150‑1 (MQ=255)
           cGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCACgacacg                                                                                                                                            <  1:1170852/150‑6 (MQ=255)
           cGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGccccc                                                                                                                                            >  1:1095889/1‑150 (MQ=255)
             ttttCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTg                                                                                                                                          <  1:188447/150‑1 (MQ=255)
                  tcAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGagag                                                                                                                                     >  1:880274/1‑150 (MQ=255)
                   cAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAgg                                                                                                                                    >  1:898811/1‑150 (MQ=255)
                    aaCCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGc                                                                                                                                   <  1:1642446/150‑1 (MQ=255)
                          acCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCaga                                                                                                                             >  1:705577/1‑150 (MQ=255)
                          acCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACATCACGCCCCCTGGAGAGGCTTCaga                                                                                                                             <  1:38697/150‑1 (MQ=255)
                           ccGCCCGTCGTCTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCagag                                                                                                                            <  2:162046/150‑1 (MQ=255)
                           ccGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCagag                                                                                                                            >  1:175012/1‑150 (MQ=255)
                            cGCCCGGCGAATTCAGGCTGAGCAGTTTTGACTTAATGACATAATTCGCATTCCATGAAACGTCATAAGGTTGAACCGACACGGATGGACAAATGGGGTCGTGAC‑TATCAGGCCCATGATACGTCACGCCCCCTGGAGAGGCTTCAGAGt                                                                                                                           <  2:289366/150‑1 (MQ=255)
                               ccGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGCACCCTGCAACGTCACGCCCCCGGGAGAGGCTTCAGAGtcct                                                                                                                        >  1:461199/1‑148 (MQ=255)
                                 gTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcattc                                                                                                                      <  1:1648071/150‑1 (MQ=255)
                                   cGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcattcat                                                                                                                    >  1:1383415/1‑150 (MQ=255)
                                    gACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCata                                                                                                                   <  1:448792/150‑1 (MQ=255)
                                      cTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATTAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCataat                                                                                                                 <  1:931323/150‑1 (MQ=255)
                                      cTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCataat                                                                                                                 <  2:185021/150‑1 (MQ=255)
                                      cTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCataat                                                                                                                 <  1:975799/150‑1 (MQ=255)
                                      cTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCataat                                                                                                                 >  1:359173/1‑150 (MQ=255)
                                          aGGGTGAGCAGTTTTGACTTACTGCCATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGACCCCTGGAGAGGATTCAGAGTCATTCATAATaacc                                                                                                             >  1:268842/1‑150 (MQ=255)
                                          aGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaacc                                                                                                             <  2:880274/150‑1 (MQ=255)
                                            ggTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaaccaa                                                                                                           <  2:1345806/150‑1 (MQ=255)
                                             gTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATTGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaaccaac                                                                                                          <  1:320628/150‑1 (MQ=255)
                                             gTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaaccaac                                                                                                          <  2:673052/150‑1 (MQ=255)
                                              tGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaaccaacc                                                                                                         >  1:153237/1‑150 (MQ=255)
                                                 gCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACATATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATg                                                                                                      <  2:546108/150‑1 (MQ=255)
                                                 gCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATg                                                                                                      >  2:465322/1‑150 (MQ=255)
                                                        tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATCACCAACCATGGACAAGa                                                                                               >  2:486396/1‑150 (MQ=255)
                                                        tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGa                                                                                               >  2:1023422/1‑150 (MQ=255)
                                                         gACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGaa                                                                                              >  2:1636805/1‑150 (MQ=255)
                                                          aCTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCACCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGaaa                                                                                             <  1:220887/150‑1 (MQ=255)
                                                            ttACTGACCTAGTTGTCATGCCATGAACAGTGATCAGGTTGATCCGAACCGGCGGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTGCAGAGTCTTTGATATTAACCTACCATGGACAAGAAAAg                                                                                           <  2:131530/150‑1 (MQ=255)
                                                            ttACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAg                                                                                           <  1:1314856/150‑1 (MQ=255)
                                                                   cATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACTCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGAtt                                                                                    >  1:569540/1‑150 (MQ=255)
                                                                     tAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCa                                                                                  <  2:603266/150‑1 (MQ=255)
                                                                      aaTTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCaa                                                                                 >  1:749594/1‑150 (MQ=255)
                                                                            cATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTg                                                                           >  1:550415/1‑150 (MQ=255)
                                                                            cATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATAAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCAATCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTg                                                                           >  1:131213/1‑150 (MQ=255)
                                                                            cATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGGGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTg                                                                           <  2:1242356/150‑1 (MQ=255)
                                                                             aTTCCATGCAACGTCATAAGGTTGAACAGGTGCGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCGCGGTGGAGAGGCTTCTGAGTCATTGATTATAACCAGCCTGGGCCTAGATATGTTAGGTTCAGTCGCTgg                                                                          >  2:279128/1‑149 (MQ=255)
                                                                             aTTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc                                                                          >  1:366820/1‑150 (MQ=255)
                                                                              ttCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCa                                                                         <  1:930700/150‑1 (MQ=255)
                                                                               tCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCag                                                                        <  1:373393/150‑1 (MQ=255)
                                                                                ccATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGTCCACTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCata                                                                       >  2:174697/1‑148 (MQ=255)
                                                                                 cctGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCagag                                                                      <  2:1435488/148‑1 (MQ=255)
                                                                                 cctGAAACGTCATAACGTGCCACCTACACGGACGTACAACTCGGGGCGTGAC‑TATCAGGCCCCGGATACGTCACGCCCCCGGTAGAGGCTGCAGAGTCATTCATAATAACCAACCATGGACAAGAAATGTAAGATTCAATCTCTGCagcg                                                                      <  2:120733/148‑3 (MQ=255)
                                                                                 cATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCagag                                                                      <  1:477923/150‑1 (MQ=255)
                                                                                  aTGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAgaga                                                                     >  2:928517/1‑150 (MQ=255)
                                                                                   tGAAACGTCATAAGGTTTAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGGCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGAc                                                                    <  1:1072671/150‑1 (MQ=255)
                                                                                   tGAAACGTCATAAGGTTGAACAGAAACGGATTGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGAc                                                                    <  2:1051849/150‑1 (MQ=255)
                                                                                   tGAAACGTCATAAGGATGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCATCACAGTCATTCATAATAACCAACCAAGGACAAGAAAAATAATATACAATCTCTGCAGacac                                                                    >  1:351443/1‑150 (MQ=255)
                                                                                     aaaCGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACaa                                                                  <  1:1550585/150‑1 (MQ=255)
                                                                                         gTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGTAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAc                                                              <  2:830136/150‑1 (MQ=255)
                                                                                         gTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAc                                                              <  2:546270/150‑1 (MQ=255)
                                                                                          tCATAAGGTTGAAGAGAATCTGATTGACAAATCGTGTCGTGAC‑TATCAGGCACCTGATACGTCACGCACCCTGTAGAGGCTTAAGAGTCATTCTTAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAcc                                                             <  1:255224/150‑1 (MQ=255)
                                                                                          tCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAcc                                                             >  2:707630/1‑150 (MQ=255)
                                                                                          tCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAcc                                                             <  1:1636805/150‑1 (MQ=255)
                                                                                            aTAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTc                                                           <  2:175012/150‑1 (MQ=255)
                                                                                             aaaGGTTGAACATAAACGGAATGTAAAATCGTGTAGTTTC‑TAATATGGCCCTGATACCTCACGGCCCCTGGAGATGCTTCACATTGATTTAAAAAAACCAAACATTGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCa                                                          <  2:882532/149‑1 (MQ=255)
                                                                                                  ttGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGt                                                     >  1:203768/1‑150 (MQ=255)
                                                                                                   tGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGta                                                    >  1:417676/1‑150 (MQ=255)
                                                                                                   tGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGta                                                    >  1:1176550/1‑150 (MQ=255)
                                                                                                    gAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCGTGCATAATAACCCACCATGGACAAGAAAAGTAAGATTAAATCTCTGCAGAGACAATTAACTCAAATGtct                                                   >  2:320663/1‑148 (MQ=255)
                                                                                                      aCAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTataa                                                 >  1:172192/1‑150 (MQ=255)
                                                                                                       cAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTataat                                                <  2:898811/150‑1 (MQ=255)
                                                                                                       cAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGAGACGTCACGCCCCCTGGAGAGGCTGCAGAGTCATTAATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAACTCAAATGTataat                                                >  2:1455826/1‑150 (MQ=255)
                                                                                                           aaCGGATGGACATATCGTGTCGTGAC‑TATCAGGCACCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATCGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATAccc                                            <  1:334563/150‑1 (MQ=255)
                                                                                                             cGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGc                                          <  1:664651/150‑1 (MQ=255)
                                                                                                              ggATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCa                                         <  1:1346813/150‑1 (MQ=255)
                                                                                                                   gACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCACCCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTcc                                    >  1:379151/1‑150 (MQ=255)
                                                                                                                   gACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTcc                                    >  2:1026294/1‑150 (MQ=255)
                                                                                                                     cAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCAc                                  >  2:467843/1‑150 (MQ=255)
                                                                                                                        aTCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACAAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCa                               <  2:1581997/150‑1 (MQ=255)
                                                                                                                               cgtgACGTATCAGGGGCCTGATA‑GTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTa                        >  1:1627113/1‑150 (MQ=255)
                                                                                                                               cgtgACGTATCAGGGGCCTGATA‑GTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTa                        >  1:1370314/1‑150 (MQ=255)
                                                                                                                               cgtgACGTATCAGGGGCCTGATA‑GTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTa                        >  1:1373630/1‑150 (MQ=255)
                                                                                                                                gtgacgtATCAGGGGCCTGATA‑GTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCCCCCAGACGGTaa                       >  2:920667/7‑150 (MQ=255)
                                                                                                                                gtgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTaa                       <  2:1095889/150‑1 (MQ=255)
                                                                                                                                  gAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATc                     <  2:203768/150‑1 (MQ=255)
                                                                                                                                      tATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATCCCCGCAGTTCCACCAAGACGGTAATCCAg                  >  1:1142239/1‑150 (MQ=255)
                                                                                                                                      tATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCCGACGGTAATCCAg                  >  1:1105073/1‑150 (MQ=255)
                                                                                                                                       aTCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAgg                 >  2:248026/1‑150 (MQ=255)
                                                                                                                                          aGGCCCCTGATACGTCACTCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAGGACa              <  2:569540/150‑1 (MQ=255)
                                                                                                                                           ggCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCAAGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAACTCAAATGTATTATAACAGCAGTTACACCCAGACGGTTATGCAGGACaa             >  2:936887/1‑150 (MQ=255)
                                                                                                                                             ccccTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACGATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAGGACAACt           >  1:101827/1‑150 (MQ=255)
                                                                                                                                                      cGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAGGACAACTTGCTttgtt  <  1:1561755/150‑1 (MQ=255)
                                                                                                                                                      cGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAGGACAACTTGCTttgtt  <  1:710172/150‑1 (MQ=255)
                                                                                                                                                      cGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAGGACAACTTGCTttgtt  <  1:1158179/150‑1 (MQ=255)
                                                                                                                                                      cGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAGGACAACTTGCTttgtt  >  1:467963/1‑150 (MQ=255)
                                                                                                                                                      cGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACAATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAGGACAACTTGCTttgtt  <  1:174697/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAGGACAACTTGCTTTGTT  >  NC_000913/2449475‑2449773

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑