Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,004,142 A→C 10.8% T66P (ACG→CCG)  ycbF → putative periplasmic pilin chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,004,1420AC10.8% 227.8 / 3.6 116T66P (ACG→CCG) ycbFputative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base A (66/35);  new base C (3/12);  total (69/47)
Fisher's exact test for biased strand distribution p-value = 1.33e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 118 total aligned reads displayed.

TATCGCTGCTGGTGGCGGGTGAATTGCAGGCAGGAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGATCGCGAAACGCTATTT  >  NC_000913/1003993‑1004291
                                                                                                                                                     |                                                                                                                                                     
tATCGCTGCTGGTGGCGGGTGAATTGCAGGCAGGAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCaa                                                                                                                                                       >  1:1686467/1‑150 (MQ=255)
 aTCGCTGCTGGTGGCGGGTGAATTGCAGGCAGGAGTGGTGGGCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAc                                                                                                                                                      <  2:862905/150‑1 (MQ=255)
    gctgctGGTGGCGGGTGAATTGCAGGCAGGAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTTCTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTATAATCACCCGCCCAAcgcg                                                                                                                                                   >  1:987568/1‑150 (MQ=255)
    gctgctGGTGGCGGGTGAATTGCAGGCAGGAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAAcgcg                                                                                                                                                   >  1:18241/1‑150 (MQ=255)
      cgctgGTGCCGTGTGAATTGCAGGCAGGTGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTTACTAATTCCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCACCGCGtt                                                                                                                                                 <  2:311960/149‑1 (MQ=255)
        ctgGTGGCGGGTGAATTGCAGGCAGGAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCACCGCGTTgg                                                                                                                                               >  2:1846185/1‑150 (MQ=255)
        ctgGTGGCGGGTGAATTGCAGGCAGGAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTgg                                                                                                                                               <  2:1662603/150‑1 (MQ=255)
         tggtggCGGGGGAATTGCAGGCAGGAGTGGGGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTggg                                                                                                                                              <  2:1528672/150‑1 (MQ=255)
          ggtggCGGGTGAATTGCAGGCAGGAGCGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGCATCATAGCTTATTAACAGTAAAATAAACCGCCCAACCCGGTGGGc                                                                                                                                             >  2:353446/1‑150 (MQ=255)
             ggCGGGTGAATTGCAGGCAGGAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCgg                                                                                                                                          >  2:1089612/1‑150 (MQ=255)
                gggTGAATGGCAGGCAGGAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCggggg                                                                                                                                       <  2:238856/150‑1 (MQ=255)
                    gAATTGCAGGCAGGAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGaa                                                                                                                                   >  1:992738/1‑150 (MQ=255)
                     aaTTGCAGGCAGGTGTGGTGGTCGGCGGACCGCGATTTATCTTTCCGGCAGACAGACAATCGATATCTATTTTACTTACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCACCCGCCCAACGCGTGGGGCCGGGGGGGAAg                                                                                                                                  <  2:1742132/150‑1 (MQ=255)
                      aTTGCAGGCAGGAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATAAACCGCCCAACGCGATTGGCCGGGGGGGAAGc                                                                                                                                 >  2:1134204/1‑150 (MQ=255)
                        tGCAGGCAGGAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGTAAGCGTc                                                                                                                              <  2:930141/150‑1 (MQ=255)
                         gcaggcagGAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTc                                                                                                                              >  1:1546583/1‑150 (MQ=255)
                          caggcaggTGTGGTGGTCGGCGGTACGCGATTTATCTTTCCGGCAGACAGAGAATCCATATCTAATTTACTGACTAATCCCAGTCAGGAATCCTGGCTTATAAACAGAAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCg                                                                                                                             <  2:750403/150‑1 (MQ=255)
                          caggcaggAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGAGGGGGAAGCGTCg                                                                                                                             <  2:763543/150‑1 (MQ=255)
                                ggAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTg                                                                                                                       >  1:1350118/1‑150 (MQ=255)
                                 gAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTAGTTTACTGACTTATACCAGTCAGGAATCCTGGCTTATTTACAGTAAAATCAACCGCCCAACGGGTTGGGCCGGGGGGGAAGCGTCGACAGTGc                                                                                                                      >  1:693250/1‑150 (MQ=255)
                                    tggtggTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGAAAGTGCCAg                                                                                                                   >  1:1303629/1‑150 (MQ=255)
                                     ggtggtCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGc                                                                                                                  >  1:1537793/1‑150 (MQ=255)
                                     ggtggtCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGAGTCGACAGTGCCAGc                                                                                                                  >  2:995311/1‑150 (MQ=255)
                                        ggtCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                                                                                                               >  2:1982966/1‑150 (MQ=255)
                                        gggcggcggAACGCGATTTATCTTTCAGGCAGACAGAGAATCTATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAAAAGTAAAATAAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAcc                                                                                                               <  1:1100472/147‑1 (MQ=255)
                                           gtgctgAACGCGATTTATCTTTCCGGCATACAGAGAATCTATATCTATTTAACATACTAATACCAGCCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCAtt                                                                                                            <  1:1421363/145‑1 (MQ=255)
                                           cggcggAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCCGCACCAtt                                                                                                            >  2:866138/1‑150 (MQ=255)
                                             gcggAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTAc                                                                                                          >  2:454678/1‑150 (MQ=255)
                                               ggAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTg                                                                                                        >  2:1710573/1‑150 (MQ=255)
                                                 aaCGCGATTTATCTTTCCGGAAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGc                                                                                                      <  1:1080103/150‑1 (MQ=255)
                                                 aaCGAGATTTATCTTTCAGGCAGACAGAGAATCGAAATCTATTTTACTGACTATTACCAGTCAGGAATCCTGGCTTATTAACAGTAAACTCAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGc                                                                                                      <  2:451760/150‑1 (MQ=255)
                                                  aCGCGATTTATATTTCCGGCAGACATAGAATCGATATCTATTTTTCTTCCTAATACCTGTCAGGAATCGTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGcc                                                                                                     <  2:1914162/150‑1 (MQ=255)
                                                   cgcgATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCg                                                                                                    >  2:1054878/1‑150 (MQ=255)
                                                   cgcgATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCg                                                                                                    >  2:1208769/1‑150 (MQ=255)
                                                   cgcgATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCg                                                                                                    >  1:260256/1‑150 (MQ=255)
                                                   cgcgATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGACGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCg                                                                                                    >  2:307803/1‑150 (MQ=255)
                                                    gcgATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAATTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGACCGCACCATTACTGGCCga                                                                                                   >  2:911234/1‑149 (MQ=255)
                                                    gcgATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCACCACCATTACTGGCCGc                                                                                                   >  2:1657460/1‑150 (MQ=255)
                                                    gcgATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAAGCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCCTTACTGGCCGc                                                                                                   >  2:1296965/1‑150 (MQ=255)
                                                       aTTTATCTTTCCGGCAGACAGACAATCGATAACTATTTTACTGACTAATCCCAGTCAGGAATCCTGGCTTATTAACACTAAAATCAACCGCCCCACGCGTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTcc                                                                                                <  1:1955229/150‑1 (MQ=255)
                                                          tATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGGTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGcc                                                                                             <  2:778992/150‑1 (MQ=255)
                                                            tCTTTCCGGCACACAGATAAGCGATATCTATTTTACTGTCTAATTCGAGTCTGGAATCCTTGCTTATTAACAGTAAAAGCAACCGCCCAACGCGTTGGTCCGGGTGGTATGCGTCGACAGTGCCAGAACCATTACTGGCCGCTCCGCCAc                                                                                           <  2:127927/150‑1 (MQ=255)
                                                               tacgggCAGACAGAGAATCGATATCTATTTTACTGACTAATCCCAGTCAGCAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTa                                                                                        <  1:1070214/146‑1 (MQ=255)
                                                                   ggCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCCCTCCGCCACTTATtca                                                                                    >  2:1401338/1‑149 (MQ=255)
                                                                     cagacagaGAATCGATATCTATTTTACTGACAAATCCCAGTCAGTAATCCTGGCTTATACACAGTAAAATCACCCGCCCACCGCGGTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTcc                                                                                  <  1:1233175/150‑1 (MQ=255)
                                                                       gacagaTAATCGATTTCTATTTTACTGACAAAAACCAGTCAGGAATCCTGGCTTATTAACAGAAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTg                                                                                <  2:136366/150‑1 (MQ=255)
                                                                         cagaGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCTCTCAGCCACTTTATCTCCTGaa                                                                              >  1:1741804/1‑150 (MQ=255)
                                                                         cagaGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGaa                                                                              >  1:1626945/1‑150 (MQ=255)
                                                                         cagaGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCAACTTAATCTCCTGaa                                                                              >  2:1955873/1‑150 (MQ=255)
                                                                              aaTCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCg                                                                         >  2:1182507/1‑150 (MQ=255)
                                                                               atctatatctatTTAACTGACAAATCCCACTCAGGAATCCTGGCTTAAAAACAGAAAAACCACCCCCCCAACGCGTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCgg                                                                        <  2:676264/150‑1 (MQ=255)
                                                                               aTCGATATCTATTTTACTGACAAATCCCCGGCAGGAATCCTGGCTTATTAACAGTAAAATCACCCCCCCACCGCGTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCgg                                                                        <  1:726073/150‑1 (MQ=255)
                                                                                    tatCTATTTTACTGACTAATACCAGTCAGGAATCCTGTCATATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCTTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGa                                                                   <  1:1326916/150‑1 (MQ=255)
                                                                                     atCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGAc                                                                  >  2:1414342/1‑150 (MQ=255)
                                                                                      tCTATTTTACTGACTACTACCAGGAAGGAATCCTGGCTTATTATCAGTAAAATCAACCGCCCAACGCGTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACa                                                                 <  1:47809/150‑1 (MQ=255)
                                                                                      tCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACGTTGTCTCCTGAAGCCCGGTACGAcc                                                                 >  2:1744522/1‑149 (MQ=255)
                                                                                       cTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAg                                                                >  1:1932295/1‑150 (MQ=255)
                                                                                        tATTTTACTGACTACTCCGAGTCAGGAATCCGGGCTTATTAACAGTAAAATCTACCGCCCAACGCGGTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAgg                                                               <  1:1207997/150‑1 (MQ=255)
                                                                                        tATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCACCGCCACTTATTCTCCTGAAGCCCGGTACGACAgg                                                               >  1:833496/1‑150 (MQ=255)
                                                                                         aTTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGc                                                              >  2:1553225/1‑150 (MQ=255)
                                                                                         aTTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGc                                                              >  2:468863/1‑150 (MQ=255)
                                                                                         aTTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGATGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCAGGTACGACAGGc                                                              >  1:770024/1‑150 (MQ=255)
                                                                                          ttttACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCa                                                             <  2:1437696/150‑1 (MQ=255)
                                                                                            ttACTGACTAATACCAGTCAGTAATCCTGGCTTATTAACATTAAAATCAACCCCCCCACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACg                                                           <  1:833752/150‑1 (MQ=255)
                                                                                               cTGACTAAAACCTCTCATGAATCCTGTCTTAAAAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTg                                                        <  2:1468677/150‑1 (MQ=255)
                                                                                                tGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGACGCCCGGTACGACAGGCCCGTTgc                                                       >  1:167337/1‑150 (MQ=255)
                                                                                                  aCTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTgcgc                                                     >  1:1594413/1‑150 (MQ=255)
                                                                                                     aaTACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCACCGCGTGTGGCCGGGGGGAAGCGTCGACAGTGCCCGCAGCCTTACAGGCCGCTCCGCCAGATATTCGACTGCCGCCCGGTGCGACCGGCCCGTGgcgagtgg                                                 >  1:1015063/1‑145 (MQ=255)
                                                                                                       tACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCt                                                >  2:437025/1‑150 (MQ=255)
                                                                                                       tACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCt                                                >  2:8641/1‑150 (MQ=255)
                                                                                                       aaCCAGTCTGGAATCCTGGATTTTTAACATAAAAATCACCCCCCCACCGCGGTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCt                                                <  1:1862298/149‑1 (MQ=255)
                                                                                                        acaaGTCAGTAATCGTGGCTTATTAAAAGTAAAATCAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTg                                               <  2:613445/147‑1 (MQ=255)
                                                                                                         ccAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCACCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTga                                              <  2:1959203/150‑1 (MQ=255)
                                                                                                          cAGTCAGGATTCCTGGCTTATTAACAGTAAATTCACCCTCCCTACGCGGTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCTCGGTACGACAGGCACGTTGCGCTTGCTgag                                             <  1:1213744/150‑1 (MQ=255)
                                                                                                            gTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGaa                                           >  2:1134279/1‑150 (MQ=255)
                                                                                                              cAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCCCCATTACTGGCCGCTCAGCCACTTATTCTACTGAAGCCCGGTACGAAAGCCACGTTCCGCTTGCTGAGAACg                                         >  2:1161204/1‑150 (MQ=255)
                                                                                                              cAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACg                                         >  2:1135970/1‑150 (MQ=255)
                                                                                                               aGGAATCCTGGCTTATTAACAGTAAAATCACCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTGATTCTCCTCAAGCCCGGTACGACCGGCCCGTTGCGCTTGCTGAGAACgg                                        >  1:228604/1‑150 (MQ=255)
                                                                                                                ggAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGGCGACAGGGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGa                                       <  2:1213837/150‑1 (MQ=255)
                                                                                                                    tCCTGGCTGATTAACAGTAAACTCACCCGCCCAACGCGGTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGc                                   <  2:283653/150‑1 (MQ=255)
                                                                                                                      cTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACTGAAAGCGa                                 >  2:1722185/1‑150 (MQ=255)
                                                                                                                       tGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGAc                                >  2:973295/1‑150 (MQ=255)
                                                                                                                          ataattaACATTAAAATAAACCCCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCAAAATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATc                             <  2:538714/147‑1 (MQ=255)
                                                                                                                             ataaaCAGAAAACCCACCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTg                          <  1:1846185/147‑1 (MQ=255)
                                                                                                                               taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCCCTTATTCTCCGGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGTACGGAAAGCGACATCTTGcc                        >  2:994851/1‑150 (MQ=255)
                                                                                                                               taACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGCAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGcc                        >  1:739276/1‑150 (MQ=255)
                                                                                                                               aaaCAGAAAAACCCACCGCCCCCCGCGTGGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGcc                        <  1:955737/149‑1 (MQ=255)
                                                                                                                                aaCAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATACTACTGAAGCCCGGTACGACAGGCACGAAGCCCTAGCTGAGAACGGAAAGCTACATCTTGCCtg                      >  2:109089/1‑150 (MQ=255)
                                                                                                                                  cAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCtgt                     >  2:1711866/1‑150 (MQ=255)
                                                                                                                                   aGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCtgtg                    >  2:215794/1‑150 (MQ=255)
                                                                                                                                     aaaaaTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACACTGACAGCACCATTACTGGCCGCACCGCCACTTAATCTCCTGAAGCCCGGTACGACAGCCACGTTGCCCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGa                  >  1:1092348/2‑150 (MQ=255)
                                                                                                                                       aaaTCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGATc                >  2:1930338/1‑150 (MQ=255)
                                                                                                                                         aTCCACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGATcgc              <  1:298222/150‑1 (MQ=255)
                                                                                                                                         aTCAACCGCCGACGGCGGTGGCGCGGGGGGGAAGGGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGATcgc              <  2:894388/150‑1 (MQ=255)
                                                                                                                                          tCAACCGCCCAACGGGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGATcgcg             <  2:1186624/150‑1 (MQ=255)
                                                                                                                                             aCCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGATCGCGaaa          >  1:1143885/1‑150 (MQ=255)
                                                                                                                                              cggCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGATCGCGATAc         >  2:858334/3‑150 (MQ=255)
                                                                                                                                               cGCCCACCGCGTTTGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGATCGCGAAACg        <  2:1380165/150‑1 (MQ=255)
                                                                                                                                                  ccAACGCGTTGGGCCGTGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTATCCTGAAGCCAGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGTAAGCGACATCATGCCTGTGGATAGCGAACAGCTa     >  2:962398/1‑150 (MQ=255)
                                                                                                                                                     aCGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGATCGCGAAACGCTAttt  >  1:634712/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TATCGCTGCTGGTGGCGGGTGAATTGCAGGCAGGAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATACCAGTCAGGAATCCTGGCTTATTAACAGTAAAATCAACCGCCCAACGCGTTGGGCCGGGGGGGAAGCGTCGACAGTGCCAGCACCATTACTGGCCGCTCCGCCACTTATTCTCCTGAAGCCCGGTACGACAGGCACGTTGCGCTTGCTGAGAACGGAAAGCGACATCTTGCCTGTGGATCGCGAAACGCTATTT  >  NC_000913/1003993‑1004291

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: