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.

TACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAGGACAACTTGCTTTGTT  >  NC_000913/2449475‑2449773
                                                                                                                                                      |                                                                                                                                                     
tACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATAc                                                                                                                                                       <  1:1536635/150‑1 (MQ=255)
tACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATAc                                                                                                                                                       <  1:133829/150‑1 (MQ=255)
    ttGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGGGCCTGATAGTCAc                                                                                                                                                  <  1:1207757/150‑1 (MQ=255)
     tGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATAGTCAc                                                                                                                                                  <  1:1376728/150‑1 (MQ=255)
     tGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAAACGTGACGTGACTATAAGGCCCCTGATACGTAAc                                                                                                                                                  >  2:334563/1‑150 (MQ=255)
     tGAAAGCGTTTTCTCAACCACACCGCCCGGCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATAGTCAc                                                                                                                                                  <  1:1370043/150‑1 (MQ=255)
      gAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACg                                                                                                                                                 >  1:1170038/1‑150 (MQ=255)
        aaGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGcc                                                                                                                                               >  1:1581997/1‑150 (MQ=255)
         aGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGccc                                                                                                                                              >  2:1402592/1‑150 (MQ=255)
          gCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGcccc                                                                                                                                             <  2:793768/150‑1 (MQ=255)
           cGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATAGTCACgacacg                                                                                                                                            <  1:1170852/150‑6 (MQ=255)
           cGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGccccc                                                                                                                                            >  1:1095889/1‑150 (MQ=255)
             ttttCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTg                                                                                                                                          <  1:188447/150‑1 (MQ=255)
                  tcAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGagag                                                                                                                                     >  1:880274/1‑150 (MQ=255)
                   cAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAgg                                                                                                                                    >  1:898811/1‑150 (MQ=255)
                    aaCCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGc                                                                                                                                   <  1:1642446/150‑1 (MQ=255)
                          acCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCaga                                                                                                                             >  1:705577/1‑150 (MQ=255)
                          acCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACATCACGCCCCCTGGAGAGGCTTCaga                                                                                                                             <  1:38697/150‑1 (MQ=255)
                           ccGCCCGTCGTCTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCagag                                                                                                                            <  2:162046/150‑1 (MQ=255)
                           ccGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCagag                                                                                                                            >  1:175012/1‑150 (MQ=255)
                            cGCCCGGCGAATTCAGGCTGAGCAGTTTTGACTTAATGACATAATTCGCATTCCATGAAACGTCATAAGGTTGAACCGACACGGATGGACAAATGGGGTCGTGACTATCAGGCCCATGATACGTCACGCCCCCTGGAGAGGCTTCAGAGt                                                                                                                           <  2:289366/150‑1 (MQ=255)
                               ccGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGCACCCTGCAACGTCACGCCCCCGGGAGAGGCTTCAGAGtcct                                                                                                                        >  1:461199/1‑148 (MQ=255)
                                 gTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcattc                                                                                                                      <  1:1648071/150‑1 (MQ=255)
                                   cGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcattcat                                                                                                                    >  1:1383415/1‑150 (MQ=255)
                                    gACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCata                                                                                                                   <  1:448792/150‑1 (MQ=255)
                                      cTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATTAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCataat                                                                                                                 <  1:931323/150‑1 (MQ=255)
                                      cTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCataat                                                                                                                 <  2:185021/150‑1 (MQ=255)
                                      cTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCataat                                                                                                                 <  1:975799/150‑1 (MQ=255)
                                      cTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCataat                                                                                                                 >  1:359173/1‑150 (MQ=255)
                                          aGGGTGAGCAGTTTTGACTTACTGCCATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGACCCCTGGAGAGGATTCAGAGTCATTCATAATaacc                                                                                                             >  1:268842/1‑150 (MQ=255)
                                          aGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaacc                                                                                                             <  2:880274/150‑1 (MQ=255)
                                            ggTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaaccaa                                                                                                           <  2:1345806/150‑1 (MQ=255)
                                             gTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATTGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaaccaac                                                                                                          <  1:320628/150‑1 (MQ=255)
                                             gTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaaccaac                                                                                                          <  2:673052/150‑1 (MQ=255)
                                              tGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaaccaacc                                                                                                         >  1:153237/1‑150 (MQ=255)
                                                 gCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACATATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATg                                                                                                      <  2:546108/150‑1 (MQ=255)
                                                 gCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATg                                                                                                      >  2:465322/1‑150 (MQ=255)
                                                        tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATCACCAACCATGGACAAGa                                                                                               >  2:486396/1‑150 (MQ=255)
                                                        tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGa                                                                                               >  2:1023422/1‑150 (MQ=255)
                                                         gACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGaa                                                                                              >  2:1636805/1‑150 (MQ=255)
                                                          aCTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCACCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGaaa                                                                                             <  1:220887/150‑1 (MQ=255)
                                                            ttACTGACCTAGTTGTCATGCCATGAACAGTGATCAGGTTGATCCGAACCGGCGGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTGCAGAGTCTTTGATATTAACCTACCATGGACAAGAAAAg                                                                                           <  2:131530/150‑1 (MQ=255)
                                                            ttACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAg                                                                                           <  1:1314856/150‑1 (MQ=255)
                                                                   cATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACTCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGAtt                                                                                    >  1:569540/1‑150 (MQ=255)
                                                                     tAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCa                                                                                  <  2:603266/150‑1 (MQ=255)
                                                                      aaTTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCaa                                                                                 >  1:749594/1‑150 (MQ=255)
                                                                            cATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTg                                                                           >  1:550415/1‑150 (MQ=255)
                                                                            cATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATAAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCAATCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTg                                                                           >  1:131213/1‑150 (MQ=255)
                                                                            cATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGGGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTg                                                                           <  2:1242356/150‑1 (MQ=255)
                                                                             aTTCCATGCAACGTCATAAGGTTGAACAGGTGCGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCGCGGTGGAGAGGCTTCTGAGTCATTGATTATAACCAGCCTGGGCCTAGATATGTTAGGTTCAGTCGCTgg                                                                          >  2:279128/1‑149 (MQ=255)
                                                                             aTTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc                                                                          >  1:366820/1‑150 (MQ=255)
                                                                              ttCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCa                                                                         <  1:930700/150‑1 (MQ=255)
                                                                               tCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCag                                                                        <  1:373393/150‑1 (MQ=255)
                                                                                ccATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGTCCACTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCata                                                                       >  2:174697/1‑148 (MQ=255)
                                                                                 cctGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCagag                                                                      <  2:1435488/148‑1 (MQ=255)
                                                                                 cctGAAACGTCATAACGTGCCACCTACACGGACGTACAACTCGGGGCGTGACTATCAGGCCCCGGATACGTCACGCCCCCGGTAGAGGCTGCAGAGTCATTCATAATAACCAACCATGGACAAGAAATGTAAGATTCAATCTCTGCagcg                                                                      <  2:120733/148‑3 (MQ=255)
                                                                                 cATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCagag                                                                      <  1:477923/150‑1 (MQ=255)
                                                                                  aTGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAgaga                                                                     >  2:928517/1‑150 (MQ=255)
                                                                                   tGAAACGTCATAAGGTTTAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGGCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGAc                                                                    <  1:1072671/150‑1 (MQ=255)
                                                                                   tGAAACGTCATAAGGTTGAACAGAAACGGATTGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGAc                                                                    <  2:1051849/150‑1 (MQ=255)
                                                                                   tGAAACGTCATAAGGATGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCATCACAGTCATTCATAATAACCAACCAAGGACAAGAAAAATAATATACAATCTCTGCAGacac                                                                    >  1:351443/1‑150 (MQ=255)
                                                                                     aaaCGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACaa                                                                  <  1:1550585/150‑1 (MQ=255)
                                                                                         gTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGTAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAc                                                              <  2:830136/150‑1 (MQ=255)
                                                                                         gTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAc                                                              <  2:546270/150‑1 (MQ=255)
                                                                                          tCATAAGGTTGAAGAGAATCTGATTGACAAATCGTGTCGTGACTATCAGGCACCTGATACGTCACGCACCCTGTAGAGGCTTAAGAGTCATTCTTAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAcc                                                             <  1:255224/150‑1 (MQ=255)
                                                                                          tCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAcc                                                             >  2:707630/1‑150 (MQ=255)
                                                                                          tCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAcc                                                             <  1:1636805/150‑1 (MQ=255)
                                                                                            aTAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTc                                                           <  2:175012/150‑1 (MQ=255)
                                                                                             aaaGGTTGAACATAAACGGAATGTAAAATCGTGTAGTTTCTAATATGGCCCTGATACCTCACGGCCCCTGGAGATGCTTCACATTGATTTAAAAAAACCAAACATTGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCa                                                          <  2:882532/149‑1 (MQ=255)
                                                                                                  ttGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGt                                                     >  1:203768/1‑150 (MQ=255)
                                                                                                   tGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGta                                                    >  1:417676/1‑150 (MQ=255)
                                                                                                   tGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGta                                                    >  1:1176550/1‑150 (MQ=255)
                                                                                                    gAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCGTGCATAATAACCCACCATGGACAAGAAAAGTAAGATTAAATCTCTGCAGAGACAATTAACTCAAATGtct                                                   >  2:320663/1‑148 (MQ=255)
                                                                                                      aCAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTataa                                                 >  1:172192/1‑150 (MQ=255)
                                                                                                       cAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTataat                                                <  2:898811/150‑1 (MQ=255)
                                                                                                       cAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGAGACGTCACGCCCCCTGGAGAGGCTGCAGAGTCATTAATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAACTCAAATGTataat                                                >  2:1455826/1‑150 (MQ=255)
                                                                                                           aaCGGATGGACATATCGTGTCGTGACTATCAGGCACCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATCGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATAccc                                            <  1:334563/150‑1 (MQ=255)
                                                                                                             cGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGc                                          <  1:664651/150‑1 (MQ=255)
                                                                                                              ggATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCa                                         <  1:1346813/150‑1 (MQ=255)
                                                                                                                   gACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCACCCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTcc                                    >  1:379151/1‑150 (MQ=255)
                                                                                                                   gACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTcc                                    >  2:1026294/1‑150 (MQ=255)
                                                                                                                     cAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCAc                                  >  2:467843/1‑150 (MQ=255)
                                                                                                                        aTCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACAAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCa                               <  2:1581997/150‑1 (MQ=255)
                                                                                                                               cgtgACGTATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTa                        >  1:1627113/1‑150 (MQ=255)
                                                                                                                               cgtgACGTATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTa                        >  1:1370314/1‑150 (MQ=255)
                                                                                                                               cgtgACGTATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTa                        >  1:1373630/1‑150 (MQ=255)
                                                                                                                                gtgacgtATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCCCCCAGACGGTaa                       >  2:920667/7‑150 (MQ=255)
                                                                                                                                gtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTaa                       <  2:1095889/150‑1 (MQ=255)
                                                                                                                                  gACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATc                     <  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)
                                                                                                                                                      |                                                                                                                                                     
TACGTTGAAAGCGTTTTCTCAACCACACCGCCCGTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTATAATACCCGCAGTTCCACCCAGACGGTAATCCAGGACAACTTGCTTTGTT  >  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: