Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,148,944 C→T 6.7% P424L (CCG→CTG)  yegD → Hsp70 chaperone family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,148,9440CT 6.7% 252.3 / 2.7 84P424L (CCG→CTG) yegDHsp70 chaperone family protein
Reads supporting (aligned to +/- strand):  ref base C (37/41);  new base T (3/3);  total (40/44)
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.98e-01

GAAAGCGCCCTCAGCCAGCCACTGACGCGGATTCTGGAACAGGTGCAACTGGCCCTGGATAACGCCCAGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCGTTAATTGAAAAATGGCTCAGCAGCGGCTTAAATCTGCCTGCGCCTGAGCCAGTTTTTCTGCTGCCTG  >  NC_000913/2148796‑2149089
                                                                                                                                                    |                                                                                                                                                 
gAAAGCGCCCTCAGCCAGCCACTGACGCGGATTCTGGAACAGGTGCAACTGGCCCTGGATAACGCCCAGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGcc                                                                                                                                                   >  1:232766/1‑149 (MQ=255)
 aaaGCGCCCTCAGCCAGCCACTGACGCGGATTCTGGAACAGGTGCAACTGGCCCTGGATAACGCCCAGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCg                                                                                                                                                  <  2:672543/149‑1 (MQ=255)
   aGCGCCCTCAGCCAGCCACTGACGCGGATTCTGGAACAGGTGCAACTGGCCCTGGATAACGCCCAGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCGGGCAGAACAGTTGCCggg                                                                                                                                                >  2:1115475/1‑149 (MQ=255)
     cgcCCTCAGCCAGCCACTGACGCGGATTCTGGAACAGGTGCAACTGGCCCTGGATAACGCCCAGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCa                                                                                                                                              <  2:1330310/149‑1 (MQ=255)
     cgcCCTCAGCCAGCCACTGACGCGGATTCTGGAACAAGTGCAACTGGCCCTGGATAGCGCCCACGAAAAACCGGGCGTTAGCTACCTGACCGGCGGTAGCGCCCGGTCTCCGCTGAGTAAAAAAGGCGTGGGAGGACAGTTGCCggcga                                                                                                                                              >  1:354971/1‑146 (MQ=255)
      gcCCTCAGCCAGCCACTGACGCGGATTCTGGAACAGGTGCAACTGGCCCTGGATAACGCCCAGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCGGAACAGTGGCCGGGCAt                                                                                                                                             >  2:1292787/1‑149 (MQ=255)
      gcCCTCAGCCAGCCACTGACGCGGATTCTGGAACAGGTGCAACTGGCCCTGGATAACGCCCAGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCAt                                                                                                                                             >  1:504133/1‑149 (MQ=255)
        ccTCAGCCAGCCACTGACGCGGATTCTGGAACAGGTGCAACTGGCCCTGGATAACGCCCAGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTc                                                                                                                                           >  1:803596/1‑149 (MQ=255)
         cTCAGCCAGCCACTGACGCGGATTCTGGAACAGGTGCAACTGGCCCTGGATAACGCCCAGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTcc                                                                                                                                          >  1:1152437/1‑149 (MQ=255)
            agccagccACTGACGCTGATTCTGGAACAGGTGCAACTGGCCCTGGATAACGCCCACGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGAt                                                                                                                                       <  2:33118/149‑1 (MQ=255)
             gccagccaCTGACGCGGATTCTGGAACAGGGGCAACTGGCCCTGGATCACGCCCAGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGAtt                                                                                                                                      <  2:1290614/149‑1 (MQ=255)
              ccagccaCTGACGCGGATTCTGGAACAGGTGCAACTGGCCCTGGATAACGCCCAGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTg                                                                                                                                     >  1:1239257/1‑149 (MQ=255)
                agcccCTGACGCTGATTGTGGAACAGGGGCAACTGGCCCCGGATAACGCCCAGGAAAAACCGGACGTTATCTCCCTGACCGGCTGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCa                                                                                                                                   <  2:974291/149‑1 (MQ=255)
                agccaCTGACGCGGATTCTGGAACAGGTGCAACTGGCCCTGGATAACGCCCAGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCa                                                                                                                                   <  2:404539/149‑1 (MQ=255)
                agccaCTGACGCGGATTCTGGAACAGGTGCAACTGGCCCTGGATAACGCCCAGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCa                                                                                                                                   >  1:114294/1‑149 (MQ=255)
                agccaCTGACGCGGATTCTGGAACAGGTGCAACTGGCCCTGGATAACGCCCAGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCa                                                                                                                                   >  2:79895/1‑149 (MQ=255)
                agccaCTGACGCGGATTCTGGAACAGGTGCAACTGGCCCTGGATAACGCCCAGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCa                                                                                                                                   <  1:1123628/149‑1 (MQ=255)
                    aCTGACGCGGATTCTGGAACAGGTGCAACTGGCCCTGGATAACGCCCAGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATCCCGATTGCAggcg                                                                                                                               >  2:800792/1‑149 (MQ=255)
                      tGACGCGGATTCTGGAACAGGTGCAACTGGCCCTGGATAACGCCCAGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAggcggc                                                                                                                             >  1:1352778/1‑149 (MQ=255)
                      tGACGCGGATTCTGGAACAGGTGCAACTGGCCCTGGATAACGCCCAGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAggcggc                                                                                                                             >  1:41376/1‑149 (MQ=255)
                       gACGCGGATTCTGGAACAGGTGCAACTGGCCCTGGATAACGCCCAGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGgcggcg                                                                                                                            <  1:145921/149‑1 (MQ=255)
                        aCGCGGATTCTGGAACAGGTGCAACTGGCCCTGGATAACGCCCAGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGa                                                                                                                           <  1:464374/149‑1 (MQ=255)
                         cgcgTATTCTGGAACAGGTGCACCTGGCCCTGGATAACGCCCAGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGAt                                                                                                                          <  2:1053077/149‑1 (MQ=255)
                          gcgGATTCTGGAACAGGTGCAACTGGCCCTGGATAACGCCCAGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAGAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATg                                                                                                                         >  1:1118277/1‑149 (MQ=255)
                            ggATTCTGGAACAGGTGCAACTGGCCCTGGATAACGCCCAGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAAGAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGAc                                                                                                                       >  1:752497/1‑149 (MQ=255)
                               ttCTGGAACAGGTGCAACTGGCCCTGGATAACGCCCAGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACttt                                                                                                                    <  1:208012/149‑1 (MQ=255)
                                  tGGAACAGGTGCAACTGGCCCTGGATAACGCCCAGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGCCTTTGGc                                                                                                                 >  2:83142/1‑149 (MQ=255)
                                  tGGAACAGGTGCAACTGACCCTGGATAACGCCCAGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGc                                                                                                                 <  1:362964/149‑1 (MQ=255)
                                   ggAACAGGTGCAACTGGCCCTGGATAACGCCCAGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCt                                                                                                                <  2:1067442/149‑1 (MQ=255)
                                      aCAGGTGCAACTGGCCCTGGATAACGCCCAGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGCCTCCg                                                                                                             >  1:1127642/1‑149 (MQ=255)
                                         ggTGCAACTGGCCCTGGATAACGCCCAGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTcc                                                                                                          >  1:475105/1‑148 (MQ=255)
                                              aaCTGGCCCTGGATAACGCCCAGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGGTTAAAAAAGCGCTGGCAGAAAAGTTGCCGGGCATTCCGATTGCAGGCGGCGCTGACTTTGGCTCCGTCCccggc                                                                                                     >  2:8155/1‑147 (MQ=255)
                                                    cccTGGATAACGCCCAGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTg                                                                                               <  1:79895/149‑1 (MQ=255)
                                                              cGCCCAGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTTATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTggg                                                                                     <  1:980209/149‑1 (MQ=255)
                                                                  cacgAAAAACCGGACTTTAGCTACCGGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGa                                                                                 <  2:879613/146‑1 (MQ=255)
                                                                   aGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGaa                                                                                <  2:114294/149‑1 (MQ=255)
                                                                   aGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGaa                                                                                <  2:1094850/149‑1 (MQ=255)
                                                                   aGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGaa                                                                                <  1:948989/149‑1 (MQ=255)
                                                                   aGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGaa                                                                                <  1:651564/149‑1 (MQ=255)
                                                                        aaaCCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAgtggt                                                                           >  1:818469/1‑149 (MQ=255)
                                                                        aaaCCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAgtggt                                                                           >  1:1197580/1‑149 (MQ=255)
                                                                            cGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGttt                                                                       <  2:41376/149‑1 (MQ=255)
                                                                            cGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGttt                                                                       <  2:1352778/149‑1 (MQ=255)
                                                                            cGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCCGGGCGGAAGTGGTGttt                                                                       >  1:340451/1‑149 (MQ=255)
                                                                             ggACGTTATCTACCTGACCGGCGGTCGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTc                                                                      >  1:675273/1‑149 (MQ=255)
                                                                                cGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCGtt                                                                   <  2:923503/149‑1 (MQ=255)
                                                                                    aTCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCGTTAAtt                                                               <  1:887397/149‑1 (MQ=255)
                                                                                       tatcTGACCGGCGGTAGCGCCCGTGCTCCGCGGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCGTTAATTGaa                                                            <  2:578923/146‑1 (MQ=255)
                                                                                         ccTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCGTTAATTGaaaa                                                          >  2:274298/1‑149 (MQ=255)
                                                                                          cTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCGTTAATTGaaaaa                                                         <  1:817270/149‑1 (MQ=255)
                                                                                          cTGACCGGCGGGAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGGTGCCGGGCATTCCGCTTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCGTTAATTGaaaaa                                                         <  1:134813/149‑1 (MQ=255)
                                                                                            gACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCGTTAATTGAAAAATg                                                       <  2:321275/149‑1 (MQ=255)
                                                                                               cggcggTAGCGCCCGTTCCCCGCTGATTAAAAAAGCGCGGGCAGAACGGTTGCCGGGGATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCGTTAATTGAAAAATGGCt                                                    <  2:1118277/149‑1 (MQ=255)
                                                                                                   ggTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCGGTAATTGAAAAATGGCTcagc                                                >  2:791962/1‑149 (MQ=255)
                                                                                                   ggTAGCGCCCGTGCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCGTTAATTGAAAAATGGCTcagc                                                <  1:524015/149‑1 (MQ=255)
                                                                                                    gTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCTGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTGCGTTAATTGAAAAGTGGCTcggca                                               >  2:717155/1‑149 (MQ=255)
                                                                                                    gTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCTGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGGTTCGTTAATTGAAAAATggacccgca                                               >  2:717132/1‑142 (MQ=255)
                                                                                                    gTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCTGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGAACGCTGGGCGGAAGTGGTGTTTCGTTAATTGAACAATGGGAcagcg                                               >  2:717170/1‑148 (MQ=255)
                                                                                                    gTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGTATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCGTTAATTGAAAAATGGCttatca                                               <  2:1037205/149‑6 (MQ=255)
                                                                                                     tAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCGTTAATTGAAAAATGGCTcagcag                                              >  1:1134569/1‑149 (MQ=255)
                                                                                                      cgcgcCCGTTCCCCGCTGATTAACAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCGTTAATTGACAAATGGCTCagcagc                                             <  2:752497/148‑1 (MQ=255)
                                                                                                         gcCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCCTTAATTGAAAAATGGCTCAGCAGCGGc                                          >  1:1118138/1‑149 (MQ=255)
                                                                                                         gcCCGGTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCGTTAATTGAAAAATGGCTCAGCAGCGGc                                          <  1:569469/149‑1 (MQ=255)
                                                                                                             gTTCTCTGCTGATTAACAAAGCGCTGGCAGAACTGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCGTTAATTGAAAAATGGCTCAGCAGCGGCTTaa                                      <  1:413388/149‑1 (MQ=255)
                                                                                                             gTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGATGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACGGCCGGGGTGGCACGCTGGGCGGAAGTGGTGTTTCGTTAATGGAAAAGTGGCTCAGCGGCGCCTTaa                                      >  2:24025/1‑149 (MQ=255)
                                                                                                                ctcCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCGTTAATTGAAAAATGGCTCAGCAGCGGCTTAAATc                                   >  2:644042/1‑149 (MQ=255)
                                                                                                                   cGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCGTTAATTGAAAAATGGCTCAGCAGCGGCTTAAATcttc                                <  2:232766/149‑3 (MQ=255)
                                                                                                                   cGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCGTTAATTGAAAAATGGCTCAGCAGCGGCTTAAATctgc                                >  1:871960/1‑149 (MQ=255)
                                                                                                                   cGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCGTTAATTGAAAAATGGCTCAGCAGCGGCTTAAATctgc                                <  1:1131375/149‑1 (MQ=255)
                                                                                                                    gCTGATTAAAAACGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGCTGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGGGTTTCGTTAATTGAAAAATGGATCAGCAGCGGCTTAAATctgcc                               >  2:176860/1‑149 (MQ=255)
                                                                                                                    gCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCGTTAATTGAAAAATGGCTCAGCAGCGGCTTCAATcggcc                               >  2:1071061/1‑149 (MQ=255)
                                                                                                                        aTTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGATGGCACGCTGGGCGGAAGGTGTGTTTCTTTCATTGAAAAATGGCTCAGCCGCGGCTTAAATCTGCCTgcg                           >  2:150562/1‑149 (MQ=255)
                                                                                                                        aTTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGAGGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGACGTGGGGTTTCGCTAATTGAAAAACGGCTCAGCAGCGGCTTAACTccgcccggg                           >  2:839830/1‑145 (MQ=255)
                                                                                                                        aTTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGAGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCGTTAATTGAAAAATGGCTCAGCAGCGGCTTAAATCTGCCTgcg                           >  1:930369/1‑149 (MQ=255)
                                                                                                                            aaaaaGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCGTTAATTGCAAAATGGCTCAGCAGCGGCTTAAATCTGCCTGCGCCTg                       >  2:392873/1‑149 (MQ=255)
                                                                                                                               aaGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCGTTAATTGAAAAATGGCTCAGCAGCGGCTTAAATCTGCCTGCGCCTGAGc                    <  2:442579/149‑1 (MQ=255)
                                                                                                                                    cTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCGTTAATTGACAAATGGCTCAGCAGCGGCTTAAATCTGCCTGCGCCTGAGCCAGtt               <  2:1087298/149‑1 (MQ=255)
                                                                                                                                       gCAGAACAGGTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCGTTAATTGAAAAATGGCTCAGCAGCGGCTTAAATCTGCCTGCGCCTGAGCCAGttttt            <  1:336768/149‑1 (MQ=255)
                                                                                                                                        cAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCGTTAATTGAAAAATGGCTCAGCAGCGGCTTAAATCTGCCTGCGCCTGAGCCAGTTTTTc           <  1:800792/149‑1 (MQ=255)
                                                                                                                                         aGAACAGTTGCTGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCGTTTATTGAAAAATGGCTCAGCAGCGGCTTAAACCTGCCTGCGCCTGAGCCAGTTTTTct          <  1:717170/149‑1 (MQ=255)
                                                                                                                                         aGAACAGTTGCTGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTCCGTTAATTGAAAAATGGCTCAGCAGCGGCTTAAATCTGCCTGCGCCTGAGCCAGTTTTTct          <  1:717132/149‑1 (MQ=255)
                                                                                                                                         aGAACAGTTGCTGGGCATTCCGAATGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCGTTAATTGAAAAATGGCTCAGCAGCGGCTTAAATCTGCCTGCGCCTGAGCCAGTTTTTct          <  1:717155/149‑1 (MQ=255)
                                                                                                                                         aGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCGTTAATTGAAAAATGGCTCAGCAGCGGCTTAAATCTGCCTGCGCCTGGGCCAGTTTTTct          >  1:223837/1‑149 (MQ=255)
                                                                                                                                         aGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCCCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCGTTAATTGAAAAATGGCTCAGCAGCGGCTTAAATCTGCCTGCGCCTGAGCCAGTTTTTct          <  1:33634/149‑1 (MQ=255)
                                                                                                                                                 tGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCGTTAATTGAAAAATGGCTCAGCAGCGGCTTAAATCTGCCTGCGCCTGAGCCAGTTTTTCTGCTGCCTg  >  2:1344180/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GAAAGCGCCCTCAGCCAGCCACTGACGCGGATTCTGGAACAGGTGCAACTGGCCCTGGATAACGCCCAGGAAAAACCGGACGTTATCTACCTGACCGGCGGTAGCGCCCGTTCTCCGCTGATTAAAAAAGCGCTGGCAGAACAGTTGCCGGGCATTCCGATTGCAGGCGGCGATGACTTTGGCTCCGTCACCGCCGGGCTGGCACGCTGGGCGGAAGTGGTGTTTCGTTAATTGAAAAATGGCTCAGCAGCGGCTTAAATCTGCCTGCGCCTGAGCCAGTTTTTCTGCTGCCTG  >  NC_000913/2148796‑2149089

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