Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,417,999 C→A 6.3% noncoding (2144/2904 nt) rrlD ← 23S ribosomal RNA of rrnD operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,417,9990CA 6.3% 917.3 / 33.8 294noncoding (2144/2904 nt)rrlD23S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base C (139/136);  new base A (9/10);  total (148/146)
Fisher's exact test for biased strand distribution p-value = 8.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 298 total aligned reads displayed.

TCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCGGGTGGA  >  CP009273/3417862‑3418148
                                                                                                                                         |                                                                                                                                                     
tcctccGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCt                                                                                                                                           <  1:5002605/150‑1 (MQ=11)
  ctccGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTaa                                                                                                                                         <  2:4729670/150‑1 (MQ=32)
     cGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGcc                                                                                                                                      >  1:3578727/1‑150 (MQ=32)
     cGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGcc                                                                                                                                      >  1:906629/1‑150 (MQ=32)
       ttACTCTTTCTGAGGAGACCGCGCCAGTAAATCTACGCACCAGACACTGTCCGCAACCAGGATTTCGGGTCAACTTTAGAACTTCAAACATTAAAGGGTGGTATTTCAAGTTCGGCTCCATGCATACTGGCGTCCACTCTTCTAATCCTc                                                                                                                                    <  2:2848787/150‑1 (MQ=11)
         aCTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTccc                                                                                                                                  >  2:98751/1‑150 (MQ=32)
          ctctTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCa                                                                                                                                 <  2:1003592/150‑1 (MQ=32)
           tctTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCAc                                                                                                                                <  1:472908/150‑1 (MQ=32)
            ctTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCAcc                                                                                                                               <  1:2097847/150‑1 (MQ=32)
            ctTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCAcc                                                                                                                               >  1:2800879/1‑150 (MQ=32)
             tttAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCt                                                                                                                              <  2:4466629/150‑1 (MQ=32)
                  gaggagACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCt                                                                                                                         >  1:4614906/1‑150 (MQ=32)
                   aggagACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTa                                                                                                                        <  2:2626820/150‑1 (MQ=32)
                   aggagACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGAGTCCACACTTCTAAGCCTCCCACCTATCCTa                                                                                                                        <  2:4877242/150‑1 (MQ=17)
                     gagaCCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTaca                                                                                                                      <  1:799875/150‑1 (MQ=32)
                     gagaCCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTaca                                                                                                                      <  2:1593956/150‑1 (MQ=32)
                       gaCCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTAcaca                                                                                                                    <  1:4903301/150‑1 (MQ=32)
                            ccccAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACTTCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAg                                                                                                               <  1:3467470/150‑1 (MQ=17)
                            ccccAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAg                                                                                                               <  2:1208563/150‑1 (MQ=32)
                             cccAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAgg                                                                                                              <  2:1764130/150‑1 (MQ=32)
                              ccAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGc                                                                                                             >  1:2719028/1‑150 (MQ=32)
                               cAGTCAAACTACCCACCAGACACTTTCCGCAAACCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCAATCCTCCACATCAAGGCt                                                                                                            >  1:37204/1‑150 (MQ=255)
                                 gTCAAACTACCCCCCAGACCCTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCACACATTAAAGGGTGGTATGTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCa                                                                                                          <  2:3620571/150‑1 (MQ=255)
                                  tCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCaa                                                                                                         >  2:3303778/1‑150 (MQ=32)
                                        tACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCa                                                                                                   >  2:1875975/1‑150 (MQ=32)
                                        tACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCa                                                                                                   >  1:2652227/1‑150 (MQ=32)
                                        tACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTGCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCa                                                                                                   >  2:3698993/1‑150 (MQ=17)
                                         aCCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAg                                                                                                  <  2:656393/150‑1 (MQ=32)
                                         aCCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAg                                                                                                  >  2:2822246/1‑150 (MQ=32)
                                          cccACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAgt                                                                                                 <  2:2743606/150‑1 (MQ=32)
                                            caccaGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAgtgt                                                                                               <  1:520973/150‑1 (MQ=32)
                                             accaGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTACAGGGTGGTATTTCAAGGTCGGCTCCATGCACACTGGCTTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTc                                                                                              <  1:4236136/150‑1 (MQ=12)
                                             accaGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTc                                                                                              <  2:4653758/150‑1 (MQ=32)
                                               caGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCaa                                                                                            <  2:724008/150‑1 (MQ=32)
                                               caGACACGGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCaa                                                                                            <  2:4418970/150‑1 (MQ=17)
                                                   cacTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCta                                                                                        <  2:1133185/150‑1 (MQ=32)
                                                   cacTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGAGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCta                                                                                        >  2:2350435/1‑150 (MQ=17)
                                                      tGTCCGCAAACCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGAAGACTGGCGTCCACACTTCTAAGCCTACCACCTATCCTACACAACAACGCTCAATGTTCAGAGTCAAGCtatat                                                                                     >  1:203988/1‑149 (MQ=11)
                                                       gTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGt                                                                                    >  2:3907590/1‑150 (MQ=32)
                                                        tCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTa                                                                                   >  2:939150/1‑150 (MQ=32)
                                                         ccGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTaa                                                                                  <  2:4223918/150‑1 (MQ=32)
                                                         ccGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTaa                                                                                  <  1:2510471/150‑1 (MQ=32)
                                                          cGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTaaa                                                                                 >  1:3435391/1‑150 (MQ=32)
                                                            cAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAgg                                                                               >  2:3582350/1‑150 (MQ=32)
                                                                  ggATGACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACg                                                                         <  2:893558/150‑1 (MQ=17)
                                                                      tACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGt                                                                     >  1:3479152/1‑150 (MQ=32)
                                                                      tACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGt                                                                     >  2:2272803/1‑150 (MQ=32)
                                                                      tACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGt                                                                     <  1:325510/150‑1 (MQ=32)
                                                                      tACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGt                                                                     <  1:1337855/150‑1 (MQ=32)
                                                                      tACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGAGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGt                                                                     <  1:1951093/150‑1 (MQ=17)
                                                                       aCGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTc                                                                    >  1:1859120/1‑150 (MQ=32)
                                                                       aCGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGCTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTc                                                                    <  2:3144300/150‑1 (MQ=17)
                                                                        cGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCt                                                                   <  1:3340694/150‑1 (MQ=32)
                                                                         gggTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCtt                                                                  >  1:2007380/1‑150 (MQ=32)
                                                                         gggTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCtt                                                                  >  1:349019/1‑150 (MQ=32)
                                                                          ggTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCttt                                                                 <  1:1680927/150‑1 (MQ=32)
                                                                           gTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTc                                                                <  1:1550443/150‑1 (MQ=32)
                                                                            tCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTcc                                                               >  1:882936/1‑150 (MQ=32)
                                                                             cAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCg                                                              <  2:2990627/150‑1 (MQ=32)
                                                                              aaCGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGt                                                             >  2:396942/1‑150 (MQ=32)
                                                                                cGTTCGAACATCAAACATGACAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCt                                                           <  2:3376502/150‑1 (MQ=12)
                                                                                  ttAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTg                                                         <  1:1415075/150‑1 (MQ=32)
                                                                                    aGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGcc                                                       >  2:4001036/1‑150 (MQ=32)
                                                                                    aGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGcc                                                       >  1:4760283/1‑150 (MQ=32)
                                                                                     gAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCcg                                                      >  2:3464562/1‑150 (MQ=32)
                                                                                      taCATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCcgc                                                     >  1:4140593/2‑150 (MQ=17)
                                                                                      aaCATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCcgc                                                     >  2:3392541/1‑150 (MQ=32)
                                                                                          tCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGt                                                 <  1:4263132/150‑1 (MQ=32)
                                                                                            aaaCATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTac                                               >  2:1084435/1‑150 (MQ=32)
                                                                                            aaaCATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTac                                               >  1:950133/1‑150 (MQ=32)
                                                                                             aaCATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTaca                                              <  2:768026/150‑1 (MQ=32)
                                                                                               cATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACt                                            >  1:1056338/1‑150 (MQ=32)
                                                                                               cATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACt                                            >  2:3229060/1‑150 (MQ=32)
                                                                                                     aGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCt                                      >  1:656472/1‑150 (MQ=32)
                                                                                                       ggtggtATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTc                                    >  2:203930/1‑150 (MQ=32)
                                                                                                       ggtggtATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTc                                    <  2:5106632/150‑1 (MQ=32)
                                                                                                          ggtATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTcaca                                 <  2:2133628/150‑1 (MQ=32)
                                                                                                           gtATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAg                                <  1:3325958/150‑1 (MQ=32)
                                                                                                            tATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGc                               >  2:1939957/1‑150 (MQ=32)
                                                                                                            tATTTCAAGGTCGGCTCCATGCAGACTGGAGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGc                               <  1:4529989/150‑1 (MQ=17)
                                                                                                             aTTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCg                              >  1:391085/1‑150 (MQ=32)
                                                                                                               ttCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAg                            >  2:2808058/1‑150 (MQ=32)
                                                                                                                 cAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGtt                          >  1:1503376/1‑150 (MQ=32)
                                                                                                                  aaGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTc                         >  2:1802544/1‑150 (MQ=32)
                                                                                                                  aaGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTc                         >  2:5050559/1‑150 (MQ=32)
                                                                                                                   aGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCa                        >  1:3948337/1‑150 (MQ=32)
                                                                                                                     gTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCGCCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAAt                      >  1:2344001/1‑150 (MQ=17)
                                                                                                                        ggCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTc                   >  1:102972/1‑150 (MQ=32)
                                                                                                                        ggCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTc                   >  2:4786587/1‑150 (MQ=32)
                                                                                                                            ccATGCAGACTGGAGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAGTGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTg               >  1:3280162/1‑150 (MQ=17)
                                                                                                                                 cAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGtct          <  2:3055619/150‑1 (MQ=32)
                                                                                                                                  aGACTGGAGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAAATTCACTGAGtctc         >  1:1109314/1‑150 (MQ=17)
                                                                                                                                     cTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCggg      <  1:5038172/150‑1 (MQ=32)
                                                                                                                                      tGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCGGGt     >  1:3091319/1‑150 (MQ=32)
                                                                                                                                      tGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCGGGt     >  2:857319/1‑150 (MQ=32)
                                                                                                                                       ggCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCGGGTg    <  1:2970868/150‑1 (MQ=32)
                                                                                                                                        gCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTTCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCGGGTgg   >  1:979051/1‑150 (MQ=17)
                                                                                                                                         cGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCGGGTGga  >  2:4952581/1‑150 (MQ=32)
                                                                                                                                         cGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCGGGTGga  >  2:3488618/1‑150 (MQ=32)
                                                                                                                                         cGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCGGGTGga  <  2:2505740/150‑1 (MQ=32)
                                                                                                                                         |                                                                                                                                                     
TCCTCCGTTACTCTTTAGGAGGAGACCGCCCCAGTCAAACTACCCACCAGACACTGTCCGCAACCCGGATTACGGGTCAACGTTAGAACATCAAACATTAAAGGGTGGTATTTCAAGGTCGGCTCCATGCAGACTGGCGTCCACACTTCTAAGCCTCCCACCTATCCTACACATCAAGGCTCAATGTTCAGTGTCAAGCTATAGTAAAGGTTCACGGGGTCTTTCCGTCTTGCCGCGGGTACACTGCATCTTCACAGCGAGTTCAATTTCACTGAGTCTCGGGTGGA  >  CP009273/3417862‑3418148

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: