Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,167,962 G→T 17.4% D98Y (GAT→TAT)  nusG → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,167,9620GT17.4% 773.2 / 100.3 262D98Y (GAT→TAT) nusGtranscription termination factor
Reads supporting (aligned to +/- strand):  ref base G (85/131);  new base T (23/23);  total (108/154)
Fisher's exact test for biased strand distribution p-value = 1.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 265 total aligned reads displayed.

AGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGT  >  CP009273/4167814‑4168110
                                                                                                                                                    |                                                                                                                                                    
 gTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGa                                                                                                                                                     <  2:1837337/149‑1 (MQ=255)
  tGGTTGAAATCCGGGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGAt                                                                                                                                                    <  1:4120396/149‑1 (MQ=255)
       gAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTcc                                                                                                                                               <  1:3107283/149‑1 (MQ=255)
       gAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTcc                                                                                                                                               <  1:3613926/149‑1 (MQ=255)
        aaaTCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCt                                                                                                                                              <  2:4516207/149‑1 (MQ=255)
          aTCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTgc                                                                                                                                            >  2:3522002/1‑149 (MQ=255)
              gTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCa                                                                                                                                        >  2:3103565/1‑149 (MQ=255)
               tGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCaa                                                                                                                                       <  1:795454/149‑1 (MQ=255)
                 gCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATc                                                                                                                                     <  1:3021420/149‑1 (MQ=255)
                 gCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTCGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATc                                                                                                                                     <  2:3577044/149‑1 (MQ=255)
                 gCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCACTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATc                                                                                                                                     <  1:4286635/149‑1 (MQ=255)
                   ggTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCTATCGTCCTGCGCCAATCAg                                                                                                                                   >  1:1009033/1‑149 (MQ=255)
                     tCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCg                                                                                                                                 >  2:2363655/1‑149 (MQ=255)
                       aGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGAt                                                                                                                               <  1:3454360/149‑1 (MQ=255)
                             gCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAaagaa                                                                                                                         >  2:1259636/1‑149 (MQ=255)
                             gCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAaagaa                                                                                                                         <  1:588416/149‑1 (MQ=255)
                             gCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAaagaa                                                                                                                         >  1:341171/1‑149 (MQ=255)
                             gCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAaagaa                                                                                                                         <  2:932760/149‑1 (MQ=255)
                             gCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAaagaa                                                                                                                         <  2:3322291/149‑1 (MQ=255)
                             gCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAaagaa                                                                                                                         >  1:1169769/1‑149 (MQ=255)
                               aaaaGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGt                                                                                                                       <  2:2224644/149‑1 (MQ=255)
                               aaaaGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGt                                                                                                                       <  2:3189272/149‑1 (MQ=255)
                                aaaGCGAACGTAAATTATTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTc                                                                                                                      >  1:4073626/1‑149 (MQ=255)
                                  aGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGa                                                                                                                    <  2:81405/149‑1 (MQ=255)
                                    cGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATg                                                                                                                  <  2:2917681/149‑1 (MQ=255)
                                      aaCGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCTATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCg                                                                                                                >  1:2368171/1‑149 (MQ=255)
                                      aaCGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCg                                                                                                                >  1:4557020/1‑149 (MQ=255)
                                       aCGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCTATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGa                                                                                                               >  2:1325253/1‑149 (MQ=255)
                                       aCGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGa                                                                                                               <  2:2482939/149‑1 (MQ=255)
                                        cGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGGGGGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGAt                                                                                                              >  2:3765194/1‑149 (MQ=255)
                                          tAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTa                                                                                                            >  1:2927732/1‑149 (MQ=255)
                                           aaaTTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTAt                                                                                                           <  1:419204/149‑1 (MQ=255)
                                           aaaTTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTAt                                                                                                           <  1:3198486/149‑1 (MQ=255)
                                              ttcttcCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGaa                                                                                                        >  2:1413360/1‑149 (MQ=255)
                                                cttcCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAAcc                                                                                                      <  1:2270665/149‑1 (MQ=255)
                                                    ccTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCt                                                                                                  >  2:1740681/1‑149 (MQ=255)
                                                        gCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTgcag                                                                                              >  1:963408/1‑149 (MQ=255)
                                                           aCGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCTATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGcagcag                                                                                           <  2:1503266/149‑1 (MQ=255)
                                                           aCGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCTATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGcagcag                                                                                           <  2:1504032/149‑1 (MQ=255)
                                                           aCGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGcagcag                                                                                           >  1:4001436/1‑149 (MQ=255)
                                                           aCGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGcagcag                                                                                           >  1:2618832/1‑149 (MQ=255)
                                                             gTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGt                                                                                         <  2:3064122/149‑1 (MQ=255)
                                                             gTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGt                                                                                         >  1:1067443/1‑149 (MQ=255)
                                                              tCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGtt                                                                                        <  2:2570196/149‑1 (MQ=255)
                                                              tCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGtt                                                                                        <  2:4608106/149‑1 (MQ=255)
                                                                 tCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGt                                                                                     >  1:856235/1‑149 (MQ=255)
                                                                 tCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGt                                                                                     >  1:3245653/1‑149 (MQ=255)
                                                                    ttCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCTATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGAt                                                                                  >  2:1805980/1‑149 (MQ=255)
                                                                      cAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATaa                                                                                <  1:4631509/149‑1 (MQ=255)
                                                                        gATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCTATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGc                                                                              >  2:895296/1‑149 (MQ=255)
                                                                        gATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGc                                                                              <  2:2810045/149‑1 (MQ=255)
                                                                        gATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGc                                                                              <  2:3240630/149‑1 (MQ=255)
                                                                           ggTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCcgc                                                                           <  1:1580108/149‑1 (MQ=255)
                                                                             tgatgaACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGt                                                                         >  1:1282414/1‑149 (MQ=255)
                                                                               atgaACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTcc                                                                       <  1:2770770/149‑1 (MQ=255)
                                                                                  aaCGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGaa                                                                    <  2:1136753/149‑1 (MQ=255)
                                                                                   acgacgCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCTATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGaaa                                                                   >  1:1471022/1‑149 (MQ=255)
                                                                                    cgacgCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCTATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGaaaa                                                                  <  2:365339/149‑1 (MQ=255)
                                                                                        gcgAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCt                                                              >  2:4104633/1‑149 (MQ=255)
                                                                                        gcgAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCt                                                              >  2:3407874/1‑149 (MQ=255)
                                                                                          gAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGt                                                            >  1:4454435/1‑149 (MQ=255)
                                                                                           aGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGtt                                                           <  1:4751622/149‑1 (MQ=255)
                                                                                              tGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGa                                                        <  2:2775142/149‑1 (MQ=255)
                                                                                               ggCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGaa                                                       >  1:561505/1‑149 (MQ=255)
                                                                                               ggCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGaa                                                       <  2:693370/149‑1 (MQ=255)
                                                                                                    cTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCggg                                                  <  1:1382724/149‑1 (MQ=255)
                                                                                                    cTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCggg                                                  <  1:3771973/149‑1 (MQ=255)
                                                                                                         gcgcAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGaaa                                             <  2:1127945/149‑1 (MQ=255)
                                                                                                         gcgcAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGaaa                                             <  1:2177936/149‑1 (MQ=255)
                                                                                                         gcgcAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGaaa                                             <  1:2675687/149‑1 (MQ=255)
                                                                                                          cgcAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAAt                                            <  1:1952459/149‑1 (MQ=255)
                                                                                                           gcAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATg                                           <  2:77732/149‑1 (MQ=255)
                                                                                                           gcAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATg                                           <  1:4558078/149‑1 (MQ=255)
                                                                                                           gcAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATg                                           <  1:1766388/149‑1 (MQ=255)
                                                                                                            cAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCACTCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATgg                                          <  1:173747/149‑1 (MQ=255)
                                                                                                              gCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCTATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTc                                        >  1:115483/1‑149 (MQ=255)
                                                                                                              gCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTc                                        <  2:4111089/149‑1 (MQ=255)
                                                                                                               cGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTcc                                       <  1:1605233/149‑1 (MQ=255)
                                                                                                                 tACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCgt                                     <  1:3183085/149‑1 (MQ=255)
                                                                                                                  aCCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGGCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCgtg                                    <  2:3272034/149‑1 (MQ=255)
                                                                                                                   ccGCGTGTGATGGGCTTCATCGGCGGTACTTCCTATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCgtgt                                   <  2:1077921/149‑1 (MQ=255)
                                                                                                                        tgtgATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAatg                              >  2:2082642/1‑149 (MQ=255)
                                                                                                                        tgtgATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAatg                              <  1:3526544/149‑1 (MQ=255)
                                                                                                                         gtgATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAatga                             >  1:635718/1‑149 (MQ=255)
                                                                                                                         gtgATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGAGTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACTCTGTGTGAACCGGGTGAAATGGTCCGTGTTAatgg                             >  1:1224588/1‑148 (MQ=255)
                                                                                                                          tgATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAatgat                            >  1:4094890/1‑149 (MQ=255)
                                                                                                                          tgATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAatgat                            <  2:4361095/149‑1 (MQ=255)
                                                                                                                          tgATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAatgat                            <  2:4361003/149‑1 (MQ=255)
                                                                                                                           gATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAatgatg                           >  1:2157505/1‑149 (MQ=255)
                                                                                                                                gCTTCATCGGCGGTACTTCCTATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCg                      >  2:839220/1‑149 (MQ=255)
                                                                                                                                gCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCg                      >  2:221788/1‑149 (MQ=255)
                                                                                                                                 cTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCAGGTGAAATGGTCCGTGTTAATGATGGTCCGt                     >  2:4775691/1‑149 (MQ=255)
                                                                                                                                  ttCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGtt                    <  1:1895580/149‑1 (MQ=255)
                                                                                                                                    cATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCg                  >  2:4338608/1‑149 (MQ=255)
                                                                                                                                     aTCGGCGGTACTTCCTATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGc                 <  2:2258860/149‑1 (MQ=255)
                                                                                                                                         gcggTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGAc             <  2:1996085/149‑1 (MQ=255)
                                                                                                                                           ggTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACtt           <  2:3818364/149‑1 (MQ=255)
                                                                                                                                            gTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTc          >  1:2542450/1‑149 (MQ=255)
                                                                                                                                              aCTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCaa        <  2:455962/149‑1 (MQ=255)
                                                                                                                                                 tCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACgg     <  2:1181895/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGTGGTTGAAATCCGTGGCGGTCAGCGTCGCAAAAGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGT  >  CP009273/4167814‑4168110

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: ‑