Predicted mutation
evidence seq id position mutation annotation gene description
RA AP012306 1,587,482 A→G S34P (TCC→CCC)  cspC ← stress protein, member of the CspA‑family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123061,587,4820AG99.8% 1719.9 / ‑4.5 423S34P (TCC→CCC) cspCstress protein, member of the CspA‑family
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (235/187);  total (235/188)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

Only 100 of 429 total aligned reads displayed.

ATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAG  >  AP012306/1587335‑1587628
                                                                                                                                                   |                                                                                                                                                  
  tGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAg                                                                                                                                                <  2:2142861/150‑1 (MQ=255)
   gAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGCTAGCGGGGAAgt                                                                                                                                               <  2:4881631/150‑1 (MQ=255)
   gAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAgt                                                                                                                                               >  1:3456356/1‑150 (MQ=255)
    aaGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAgtg                                                                                                                                              >  1:838512/1‑150 (MQ=255)
      gCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAgtgca                                                                                                                                            >  2:3380585/1‑148 (MQ=255)
      gCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTa                                                                                                                                            >  2:173300/1‑150 (MQ=255)
      gCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTa                                                                                                                                            <  1:3843244/150‑1 (MQ=255)
          aTTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGaa                                                                                                                                        >  1:1461612/1‑150 (MQ=255)
             cacaCTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAacac                                                                                                                                     <  1:2006490/150‑1 (MQ=255)
               cacTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACAt                                                                                                                                   >  2:1660217/1‑150 (MQ=255)
                 cTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCt                                                                                                                                 <  2:1421680/150‑1 (MQ=255)
                 cTTCAGAGCAGTGGATTCGATCAGATGGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCt                                                                                                                                 <  2:1632005/150‑1 (MQ=255)
                  ttCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCtt                                                                                                                                <  2:822664/150‑1 (MQ=255)
                  ttCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCtt                                                                                                                                >  1:2520357/1‑150 (MQ=255)
                    cAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTtg                                                                                                                              <  1:1682526/150‑1 (MQ=255)
                      gATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTtgct                                                                                                                            <  1:2823876/150‑1 (MQ=255)
                          aGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCa                                                                                                                        >  1:2219281/1‑150 (MQ=255)
                             ggATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCa                                                                                                                     >  2:2546428/1‑150 (MQ=255)
                             ggATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCa                                                                                                                     >  2:1058941/1‑150 (MQ=255)
                               aTTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGc                                                                                                                   >  1:927186/1‑150 (MQ=255)
                                  cGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCgg                                                                                                                >  1:2041019/1‑150 (MQ=255)
                                    aTCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAg                                                                                                              >  1:4550887/1‑150 (MQ=255)
                                     tCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTTGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCTTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGt                                                                                                             <  2:873978/150‑1 (MQ=255)
                                       aGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTaa                                                                                                           >  2:2668834/1‑150 (MQ=255)
                                       aGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTaa                                                                                                           <  2:4658082/150‑1 (MQ=255)
                                        gATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTTCCATCAGCCGGAGTAAt                                                                                                          >  2:1925339/1‑150 (MQ=255)
                                          tAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGa                                                                                                        >  1:3393122/1‑150 (MQ=255)
                                           aGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGaa                                                                                                       <  2:3540288/150‑1 (MQ=255)
                                            gCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAg                                                                                                      >  1:2149426/1‑150 (MQ=255)
                                             cTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGc                                                                                                     >  2:1191813/1‑150 (MQ=255)
                                                ttACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCaa                                                                                                  <  1:382752/150‑1 (MQ=255)
                                                 tACGTTAAGAGCTGCCGTACCTTTCTGGCCGTCATGAATTTCGAACTCAGCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCGGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCGGGAGTAATGAAGCCaaa                                                                                                 <  2:3815016/150‑1 (MQ=38)
                                                  aCGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCaaaa                                                                                                >  1:4455487/1‑150 (MQ=255)
                                                     ttAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCt                                                                                             >  1:3896013/1‑150 (MQ=255)
                                                      tAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCtt                                                                                            >  2:416771/1‑150 (MQ=255)
                                                       aaCAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCttt                                                                                           <  1:4620350/150‑1 (MQ=255)
                                                       aaCAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCttt                                                                                           >  1:2109008/1‑150 (MQ=255)
                                                        aCAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTa                                                                                          <  2:2804259/150‑1 (MQ=255)
                                                          aGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGa                                                                                        <  1:3199306/150‑1 (MQ=255)
                                                           gCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGAc                                                                                       >  2:1994129/1‑150 (MQ=255)
                                                            cTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACt                                                                                      >  1:1916389/1‑150 (MQ=255)
                                                                cGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGtt                                                                                  <  1:1636192/150‑1 (MQ=255)
                                                                cGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGtt                                                                                  >  2:3679314/1‑150 (MQ=255)
                                                                 ggACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTg                                                                                 >  2:570101/1‑150 (MQ=255)
                                                                  gACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGa                                                                                <  1:548270/150‑1 (MQ=255)
                                                                   aCCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGaa                                                                               >  1:481220/1‑150 (MQ=255)
                                                                   aCCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGaa                                                                               >  2:2326126/1‑150 (MQ=255)
                                                                    ccTTTCTGGCCGTCCTGATTTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAAc                                                                              <  2:4275842/150‑1 (MQ=255)
                                                                           ggCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTaa                                                                       >  2:4353007/1‑150 (MQ=255)
                                                                             ccGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAAcc                                                                     <  1:2928145/150‑1 (MQ=255)
                                                                              cGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCt                                                                    >  1:3039751/1‑150 (MQ=255)
                                                                                tCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGa                                                                  <  1:128991/150‑1 (MQ=255)
                                                                                tCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGa                                                                  >  1:178121/1‑150 (MQ=255)
                                                                                 ccTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGAc                                                                 >  2:449064/1‑150 (MQ=255)
                                                                                     aaTTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCttt                                                             <  1:1390632/150‑1 (MQ=255)
                                                                                     aaTTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCttt                                                             >  1:630701/1‑150 (MQ=255)
                                                                                      aTTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTa                                                            <  2:2447792/150‑1 (MQ=255)
                                                                                      aTTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTa                                                            >  1:2167914/1‑150 (MQ=255)
                                                                                       tttCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTaa                                                           >  1:2308113/1‑150 (MQ=255)
                                                                                         tCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATc                                                         >  1:369437/1‑150 (MQ=255)
                                                                                          cGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCt                                                        >  2:2886662/1‑150 (MQ=255)
                                                                                            aaCTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCttt                                                      <  2:3650389/150‑1 (MQ=255)
                                                                                            aaCTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCttt                                                      <  2:1352614/150‑1 (MQ=255)
                                                                                              cTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGTCCTTTAATCTTTGc                                                    <  2:1524847/150‑1 (MQ=255)
                                                                                                 aaCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt                                                 >  1:4695739/1‑150 (MQ=255)
                                                                                                 aaCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt                                                 <  2:628876/150‑1 (MQ=255)
                                                                                                  aCGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAtt                                                <  2:3246814/150‑1 (MQ=255)
                                                                                                    gTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTg                                              <  1:1012563/150‑1 (MQ=255)
                                                                                                     ttCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGa                                             <  1:3232495/150‑1 (MQ=255)
                                                                                                       cTGCCCTTCAGCCAGAGTTTTGAAGCCATTCCCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAACCCTTTAGACTCGTTGAACCACTTAACCTGCCCTTTAATCTTTGCCATTTGaaa                                           <  2:2753480/150‑1 (MQ=39)
                                                                                                        tGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaa                                          >  2:1571539/1‑150 (MQ=255)
                                                                                                          aCCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAt                                        <  2:3290832/150‑1 (MQ=255)
                                                                                                           ccTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAtt                                       <  2:3326004/150‑1 (MQ=255)
                                                                                                             ttCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTcc                                     >  1:1041106/1‑150 (MQ=255)
                                                                                                              tCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCt                                    <  2:2973398/150‑1 (MQ=255)
                                                                                                                 gCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAg                                 <  1:3552642/150‑1 (MQ=255)
                                                                                                                 gCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAg                                 >  2:4737204/1‑150 (MQ=255)
                                                                                                                  ccAGAGTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGa                                <  1:3966743/150‑1 (MQ=255)
                                                                                                                      agTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGt                            >  2:4542106/1‑150 (MQ=255)
                                                                                                                       gTTTTGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGtt                           <  1:739466/150‑1 (MQ=255)
                                                                                                                        ttttGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGttt                          >  1:298801/1‑150 (MQ=255)
                                                                                                                          ttGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTttc                        >  2:1244756/1‑150 (MQ=255)
                                                                                                                           tGAAGCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTttct                       >  2:2110112/1‑150 (MQ=255)
                                                                                                                               gCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGc                   <  1:1087562/150‑1 (MQ=255)
                                                                                                                               gCCATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGc                   >  1:3585126/1‑150 (MQ=255)
                                                                                                                                 cATTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTNCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGccc                 <  2:2084363/150‑1 (MQ=255)
                                                                                                                                  aTTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCg                <  1:3108690/150‑1 (MQ=255)
                                                                                                                                  aTTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCg                <  1:4180845/150‑1 (MQ=255)
                                                                                                                                  aTTACCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCg                >  2:3060936/1‑150 (MQ=255)
                                                                                                                                     aCCCTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAg             >  2:3951062/1‑150 (MQ=255)
                                                                                                                                      cccTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAgg            >  1:1569225/1‑150 (MQ=255)
                                                                                                                                      cccTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAgg            >  2:3163556/1‑150 (MQ=255)
                                                                                                                                      cccTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAgg            >  2:1889066/1‑150 (MQ=255)
                                                                                                                                       ccTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGc           <  2:3851778/150‑1 (MQ=255)
                                                                                                                                        cTGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCa          >  2:1492928/1‑150 (MQ=255)
                                                                                                                                         tGGATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCAt         >  1:4069433/1‑150 (MQ=255)
                                                                                                                                          ggATAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATa        >  2:43959/1‑150 (MQ=255)
                                                                                                                                             tAGCGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACa     >  2:717536/1‑150 (MQ=255)
                                                                                                                                                cGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAg  >  1:2645888/1‑150 (MQ=255)
                                                                                                                                                cGGGGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAg  >  2:2843969/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATTGAAGCGTATTCACACTTCAGATCAGTGGATTCGATCAGATAGCTGTTACGTTAACAGCTGCCGGACCTTTCTGGCCGTCCTGAATTTCGAACTCAACGTTCTGACCTTCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATTCCTTAGATTGTTTTCTTCGCCCGCAGGCATAACATAG  >  AP012306/1587335‑1587628

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑