Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,396,598 A→T 87.4% E43V (GAA→GTA)  rnr → exoribonuclease R, RNase R

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,396,5980AT87.4% 983.8 / 105.6 299E43V (GAA→GTA) rnrexoribonuclease R, RNase R
Reads supporting (aligned to +/- strand):  ref base A (22/16);  new base T (162/99);  total (184/115)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 301 total aligned reads displayed.

TGACAACGGAGGAACCGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGAAAGGTACCGTTATTGGCCAC  >  CP009273/4396452‑4396746
                                                                                                                                                  |                                                                                                                                                    
tGACAACGGAGGAACCGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTAc                                                                                                                                                    >  2:2008783/1‑149 (MQ=255)
 gACAACGGAGGAACCGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACt                                                                                                                                                   >  1:4249931/1‑149 (MQ=255)
    aaCGGAGGAACCGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGca                                                                                                                                                >  2:4383506/1‑149 (MQ=255)
    aaCGGAGGAACCGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAAAATTTAAACAAACGTTAAAAAACGGCCATCCGTGATGAGCTGGCGGTAGTACTGca                                                                                                                                                >  1:3693621/1‑149 (MQ=255)
      cGGAGGAACCGAGATGTCACAAGATCCTTTCCAGGAACTCTAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGcaca                                                                                                                                              <  1:952053/149‑1 (MQ=255)
       ggaggaACCGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAAGCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTACCCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACAt                                                                                                                                             <  1:3995004/149‑1 (MQ=255)
         aggaACCGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTg                                                                                                                                           >  2:1565399/1‑149 (MQ=255)
            aaCCGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAg                                                                                                                                        >  2:294315/1‑149 (MQ=255)
             aCCGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAgg                                                                                                                                       >  1:796899/1‑149 (MQ=255)
               cGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCg                                                                                                                                     >  1:4466981/1‑149 (MQ=255)
                gagaTGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGGATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGa                                                                                                                                    >  2:2979893/1‑149 (MQ=255)
                 agaTGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGaa                                                                                                                                   >  1:334769/1‑149 (MQ=255)
                  gaTGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAag                                                                                                                                  <  2:280376/149‑1 (MQ=255)
                  gaTGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAag                                                                                                                                  <  2:2620666/149‑1 (MQ=255)
                   aTGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAaga                                                                                                                                 <  1:2123918/149‑1 (MQ=255)
                   aTGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAaga                                                                                                                                 >  1:4516171/1‑149 (MQ=255)
                     gTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGagc                                                                                                                               >  1:1592056/1‑149 (MQ=255)
                      tCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGagca                                                                                                                              >  2:1760284/1‑149 (MQ=255)
                       cacaAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGagcag                                                                                                                             <  2:1491050/149‑1 (MQ=255)
                         caAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCt                                                                                                                           >  1:3159738/1‑149 (MQ=255)
                         caAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCt                                                                                                                           >  1:4677454/1‑149 (MQ=255)
                          aaGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCtt                                                                                                                          >  2:4532111/1‑149 (MQ=255)
                           aGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTg                                                                                                                         >  1:4012265/1‑149 (MQ=255)
                           aGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTg                                                                                                                         >  2:2640374/1‑149 (MQ=255)
                              tCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAg                                                                                                                      <  1:980450/149‑1 (MQ=255)
                               ccTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAgg                                                                                                                     >  2:4433798/1‑149 (MQ=255)
                                 tttCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGcc                                                                                                                   >  2:734791/1‑149 (MQ=255)
                                 tttCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGcc                                                                                                                   >  2:1035600/1‑149 (MQ=255)
                                    ccAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGc                                                                                                                <  1:2282472/149‑1 (MQ=255)
                                       ggAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTc                                                                                                             >  1:3253188/1‑149 (MQ=255)
                                       ggAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTc                                                                                                             >  1:4606445/1‑149 (MQ=255)
                                       ggAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTc                                                                                                             >  1:2344153/1‑149 (MQ=255)
                                             cgAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCgccgcc                                                                                                       <  1:1806771/149‑1 (MQ=255)
                                              gAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCt                                                                                                      <  1:2011550/149‑1 (MQ=255)
                                                agatGAAAAATACGCGAATACCATCCCAAGTCGGGAATTTATCCTCGAACATTTAACCAAAAGAGAAAAACCGGCCAGCCGAGATGAGCTGGAGGTAGTACTGAACATTGAAGGCGAAGAGCAGCTAGAAGGCATGCGTCTCCGCCAgc                                                                                                    >  2:2138266/4‑149 (MQ=255)
                                                aGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTgc                                                                                                    >  2:497550/1‑149 (MQ=255)
                                                 gCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTgcg                                                                                                   >  1:2894059/1‑149 (MQ=255)
                                                     aaaaaTACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGcgcgcg                                                                                               >  1:284200/1‑149 (MQ=255)
                                                     aaaaaTACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGcgcgcg                                                                                               >  2:1304833/1‑149 (MQ=255)
                                                          tACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGa                                                                                          >  2:1874711/1‑149 (MQ=255)
                                                              cgAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAGGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAgcgc                                                                                      <  1:1847291/149‑1 (MQ=255)
                                                              cgAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAgcgc                                                                                      >  1:737109/1‑149 (MQ=255)
                                                               gAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGcgcg                                                                                     >  1:2064594/1‑149 (MQ=255)
                                                               gAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGAGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGcgcg                                                                                     >  1:3846704/1‑149 (MQ=255)
                                                                aaTCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGa                                                                                    >  1:883405/1‑149 (MQ=255)
                                                                aaTCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGa                                                                                    <  2:4104746/149‑1 (MQ=255)
                                                                   cccatcccTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATgg                                                                                 <  2:2269988/149‑1 (MQ=255)
                                                                   cccatcccTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATgg                                                                                 >  1:1035584/1‑149 (MQ=255)
                                                                   cccatcacTAGTCGTTAATTTATCATCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCTAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATgg                                                                                 <  2:3358282/149‑1 (MQ=255)
                                                                    ccatcccTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGt                                                                                >  1:1155532/1‑149 (MQ=255)
                                                                       tcccTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCaa                                                                             >  2:3789955/1‑149 (MQ=255)
                                                                        cccTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAAc                                                                            <  1:2489282/149‑1 (MQ=255)
                                                                        cccTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAAc                                                                            <  1:135073/149‑1 (MQ=255)
                                                                           tAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTgg                                                                         >  1:1254337/1‑149 (MQ=255)
                                                                            aGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGt                                                                        <  1:650104/149‑1 (MQ=255)
                                                                            aGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGt                                                                        <  2:233288/149‑1 (MQ=255)
                                                                             gTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTc                                                                       >  2:4243706/1‑149 (MQ=255)
                                                                                  gAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCAc                                                                  >  2:2958392/1‑149 (MQ=255)
                                                                                  gAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCAc                                                                  >  1:4312418/1‑149 (MQ=255)
                                                                                    aTTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACtc                                                                >  2:1634052/1‑149 (MQ=255)
                                                                                      ttATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACtcgt                                                              >  1:4168938/1‑149 (MQ=255)
                                                                                        aTCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACtcgtcg                                                            >  2:2790805/1‑149 (MQ=255)
                                                                                         tCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTcgtcgt                                                           >  2:1122304/1‑149 (MQ=255)
                                                                                          ccTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCgtcgtc                                                          <  1:2645553/149‑1 (MQ=255)
                                                                                           cTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCa                                                         <  1:3947212/149‑1 (MQ=255)
                                                                                             cGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCACCTGGTCTTCACTCGTCGTCAGt                                                       >  2:3560065/1‑149 (MQ=255)
                                                                                             cGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGt                                                       <  2:4227095/149‑1 (MQ=255)
                                                                                                 cATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTa                                                   <  2:3417635/149‑1 (MQ=255)
                                                                                                  aTTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTAt                                                  <  1:2570653/149‑1 (MQ=255)
                                                                                                  aTTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTAt                                                  >  1:4223711/1‑149 (MQ=255)
                                                                                                       aCCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCt                                             >  2:3239505/1‑149 (MQ=255)
                                                                                                         cAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGc                                           >  2:3653147/1‑149 (MQ=255)
                                                                                                          aaaCGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGcc                                          >  1:864705/1‑149 (MQ=255)
                                                                                                          aaaCGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGcc                                          >  2:1871669/1‑149 (MQ=255)
                                                                                                            aCGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCgg                                        >  2:2112915/1‑149 (MQ=255)
                                                                                                               tGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAAc                                     <  2:2906076/149‑1 (MQ=255)
                                                                                                               tGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAAc                                     <  1:2217105/149‑1 (MQ=255)
                                                                                                                gAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACg                                    >  1:3773184/1‑149 (MQ=255)
                                                                                                                  aaaaCCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGcc                                  <  2:2560271/149‑1 (MQ=255)
                                                                                                                  aaaaCCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGcc                                  <  2:1119056/149‑1 (MQ=255)
                                                                                                                    aaCCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTc                                <  1:4375750/149‑1 (MQ=255)
                                                                                                                     aCCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCg                               >  2:897347/1‑149 (MQ=255)
                                                                                                                     aCCGGCCAGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCg                               >  1:4038862/1‑149 (MQ=255)
                                                                                                                      ccGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGa                              >  1:3887988/1‑149 (MQ=255)
                                                                                                                         gCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCt                           >  1:1560694/1‑149 (MQ=255)
                                                                                                                            aGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGt                        >  1:4082364/1‑149 (MQ=255)
                                                                                                                            aGCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGt                        >  1:2407400/1‑149 (MQ=255)
                                                                                                                             gCCGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTg                       <  1:3460941/149‑1 (MQ=255)
                                                                                                                              ccGTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCAGGAACGCCTCGACCTGGTGa                      <  2:2983327/149‑1 (MQ=255)
                                                                                                                              ccGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGa                      <  2:618294/149‑1 (MQ=255)
                                                                                                                                gTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGaaa                    >  1:321326/1‑149 (MQ=255)
                                                                                                                                gTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGaaa                    >  1:1888266/1‑149 (MQ=255)
                                                                                                                                gTGATGAGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGaaa                    >  1:3542397/1‑149 (MQ=255)
                                                                                                                                 tgatgaGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGAAAg                   <  1:1004297/149‑1 (MQ=255)
                                                                                                                                  gatgaGCTGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGAAAgg                  <  2:405852/149‑1 (MQ=255)
                                                                                                                                         tGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGAAAGGTACCGtt           <  1:4084657/149‑1 (MQ=255)
                                                                                                                                         tGGCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGAAAGGTACCGtt           >  1:598028/1‑149 (MQ=255)
                                                                                                                                          ggCGGTAGTACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGAAAGGTACCGTTa          >  1:4142087/1‑149 (MQ=255)
                                                                                                                                              gtagtaCTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGAAAGGTACCGTTATTgg      <  2:330879/149‑1 (MQ=255)
                                                                                                                                                  aaCTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGAAAGGTACCGTTATTGGCCAc  >  2:2397619/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TGACAACGGAGGAACCGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGAAAGGTACCGTTATTGGCCAC  >  CP009273/4396452‑4396746

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