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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,167,9950GT 6.0% 1035.3 / 24.1 267D109Y (GAT→TAT) nusGtranscription termination factor
Reads supporting (aligned to +/- strand):  ref base G (94/157);  new base T (5/11);  total (99/168)
Fisher's exact test for biased strand distribution p-value = 7.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 269 total aligned reads displayed.

AGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTGAAGAAGTGGATTACGAGAAATCTCGTCT  >  CP009273/4167848‑4168143
                                                                                                                                                   |                                                                                                                                                    
aGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCTATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGa                                                                                                                                                     <  1:2183488/149‑1 (MQ=255)
aGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGa                                                                                                                                                     <  2:2473778/149‑1 (MQ=255)
aGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGa                                                                                                                                                     <  2:1784949/149‑1 (MQ=255)
  cGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATg                                                                                                                                                   <  2:526118/149‑1 (MQ=255)
     aCGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGa                                                                                                                                                <  2:2482939/149‑1 (MQ=255)
      cGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGGGGGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGAt                                                                                                                                               >  2:3765194/1‑149 (MQ=255)
       gTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGAtt                                                                                                                                              <  2:2965663/149‑1 (MQ=255)
       gTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGAtt                                                                                                                                              >  1:3978329/1‑149 (MQ=255)
        tAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTa                                                                                                                                             >  1:2927732/1‑149 (MQ=255)
         aaaTTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTAt                                                                                                                                            <  1:419155/149‑1 (MQ=255)
         aaaTTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTAt                                                                                                                                            <  1:3198486/149‑1 (MQ=255)
          aaTTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATg                                                                                                                                           >  2:1963228/1‑149 (MQ=255)
           aTTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGa                                                                                                                                          >  1:113714/1‑149 (MQ=255)
            ttcttcCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGaa                                                                                                                                         >  1:3478641/1‑149 (MQ=255)
                tcCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGc                                                                                                                                     <  1:2333206/149‑1 (MQ=255)
                     ggCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCTATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTgca                                                                                                                                <  2:2582649/149‑1 (MQ=255)
                      gCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTgcag                                                                                                                               >  1:963408/1‑149 (MQ=255)
                       cTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTgcagc                                                                                                                              >  2:2199967/1‑149 (MQ=255)
                         aCGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCTATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGcagcag                                                                                                                            <  2:1504032/149‑1 (MQ=255)
                         aCGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGcagcag                                                                                                                            >  1:2618832/1‑149 (MQ=255)
                           gTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGt                                                                                                                          >  1:4677262/1‑149 (MQ=255)
                           gTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGt                                                                                                                          <  2:1492875/149‑1 (MQ=255)
                           gTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGt                                                                                                                          >  1:1067443/1‑149 (MQ=255)
                            tCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGtt                                                                                                                         <  2:4608106/149‑1 (MQ=255)
                               tCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGt                                                                                                                      >  1:856235/1‑149 (MQ=255)
                               tCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGt                                                                                                                      >  1:3245653/1‑149 (MQ=255)
                                      gATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGc                                                                                                               <  2:3240631/149‑1 (MQ=255)
                                      gATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGc                                                                                                               <  2:2810045/149‑1 (MQ=255)
                                       aTGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCTATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGcc                                                                                                              >  1:4457547/1‑149 (MQ=255)
                                         ggTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCcgc                                                                                                            <  1:1580108/149‑1 (MQ=255)
                                            gatgaACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCTATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTc                                                                                                         <  2:386125/149‑1 (MQ=255)
                                              tgaACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCg                                                                                                       <  2:1822418/149‑1 (MQ=255)
                                              tgaACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCg                                                                                                       >  1:1646970/1‑149 (MQ=255)
                                              tgaACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGGGTGATGGGCTGCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGAGTATGACCCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCg                                                                                                       >  2:1412061/1‑149 (MQ=255)
                                                      gcgAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCt                                                                                               >  2:3407841/1‑149 (MQ=255)
                                                      gcgAGCTGGCACCTGGTGCGCAGCGTACCGCGGGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGAGAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATCAGCCGCGTCCGAAAACGCt                                                                                               >  2:4104786/1‑149 (MQ=255)
                                                         aGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGtt                                                                                            <  1:4751622/149‑1 (MQ=255)
                                                         aGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGtt                                                                                            >  1:4115658/1‑149 (MQ=255)
                                                             ggCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGaa                                                                                        >  2:1142577/1‑149 (MQ=255)
                                                             ggCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGaa                                                                                        <  2:693370/149‑1 (MQ=255)
                                                                  cTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCggg                                                                                   <  2:3231015/149‑1 (MQ=255)
                                                                  cTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCggg                                                                                   <  1:1382724/149‑1 (MQ=255)
                                                                  cTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCggg                                                                                   <  1:3771973/149‑1 (MQ=255)
                                                                    ggTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTg                                                                                 <  1:4298110/149‑1 (MQ=255)
                                                                       gcgcAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGaaa                                                                              <  1:2675670/149‑1 (MQ=255)
                                                                       gcgcAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGaaa                                                                              <  2:1127945/149‑1 (MQ=255)
                                                                        cgcAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAAt                                                                             <  1:1952459/149‑1 (MQ=255)
                                                                         gcAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATg                                                                            <  1:4557964/149‑1 (MQ=255)
                                                                         gcAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATg                                                                            >  1:192342/1‑149 (MQ=255)
                                                                         gcAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATg                                                                            <  1:1766381/149‑1 (MQ=255)
                                                                            gCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCTATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTc                                                                         <  1:575438/149‑1 (MQ=255)
                                                                               tACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCgt                                                                      <  1:3183085/149‑1 (MQ=255)
                                                                                aCCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCgtg                                                                     >  2:841208/1‑149 (MQ=255)
                                                                                 ccGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCgtgt                                                                    <  2:1596610/149‑1 (MQ=255)
                                                                                  cgcgTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGtt                                                                   <  2:4720055/149‑1 (MQ=255)
                                                                                      tgtgATGGGCTTCATCGGCGGTACTTCCTATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAatg                                                               >  1:1455262/1‑149 (MQ=255)
                                                                                      tgtgATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAatg                                                               >  2:2082642/1‑149 (MQ=255)
                                                                                      tgtgATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAatg                                                               <  2:4557020/149‑1 (MQ=255)
                                                                                       gtgATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAatga                                                              >  1:635718/1‑149 (MQ=255)
                                                                                        tgATGGGCTTCATCGGCGGTACTTCCTATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAatgat                                                             <  2:3095101/149‑1 (MQ=255)
                                                                                        tgATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAatgat                                                             <  2:4361095/149‑1 (MQ=255)
                                                                                        tgATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAatgat                                                             >  1:4094890/1‑149 (MQ=255)
                                                                                        tgATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAatgat                                                             <  2:4361003/149‑1 (MQ=255)
                                                                                         gATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAatgatg                                                            >  1:2157505/1‑149 (MQ=255)
                                                                                            gggCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTc                                                         >  2:120634/1‑149 (MQ=255)
                                                                                            gggCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTc                                                         <  1:1184719/149‑1 (MQ=255)
                                                                                            gggCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTc                                                         <  1:3052577/149‑1 (MQ=255)
                                                                                              gCTTCATCGGCGGTACTTCCTATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCGGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCg                                                       >  2:839214/1‑149 (MQ=255)
                                                                                              gCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGACATGGTCCGTGTTAATGATGGTCCg                                                       >  2:1334954/1‑149 (MQ=255)
                                                                                              gCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCg                                                       >  2:221808/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:9065/149‑1 (MQ=255)
                                                                                                         ggTACTTCCTATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACtt                                            <  1:1229423/149‑1 (MQ=255)
                                                                                                         ggTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACtt                                            <  2:784121/149‑1 (MQ=255)
                                                                                                          gTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTc                                           >  1:2543477/1‑149 (MQ=255)
                                                                                                            aCTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCaa                                         <  2:341171/149‑1 (MQ=255)
                                                                                                                  taTCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGgtgt                                   <  1:2437939/148‑1 (MQ=255)
                                                                                                                      gTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTgttgtt                               <  1:1330711/149‑1 (MQ=255)
                                                                                                                       tCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGttgttg                              <  1:2170943/149‑1 (MQ=255)
                                                                                                                         cTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTgaa                            <  1:2911054/149‑1 (MQ=255)
                                                                                                                         cTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTgaa                            >  2:1097510/1‑149 (MQ=255)
                                                                                                                         cTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTgaa                            <  2:1009033/149‑1 (MQ=255)
                                                                                                                         cTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTgaa                            <  1:3682830/149‑1 (MQ=255)
                                                                                                                          tGCGCCAATCAGCGAGAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTgaag                           <  1:1459445/149‑1 (MQ=255)
                                                                                                                                   cAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTGAAGAAGTGGAtt                  >  2:3592398/1‑149 (MQ=255)
                                                                                                                                   cAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTGAAGAAGTGGAtt                  >  2:661644/1‑149 (MQ=255)
                                                                                                                                     gCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTGAAGTAGTGGATTAc                <  1:610728/149‑1 (MQ=255)
                                                                                                                                     gCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTGAAGAAGTGGATTAc                >  1:3586672/1‑149 (MQ=255)
                                                                                                                                     gCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTGAAGAAGTGGATTAc                <  1:276315/149‑1 (MQ=255)
                                                                                                                                     gCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTGAAGAAGTGGATTAc                >  2:2962436/1‑149 (MQ=255)
                                                                                                                                       gATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTGAAGAAGTGGATTACga              <  1:3103565/149‑1 (MQ=255)
                                                                                                                                        aTAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTGAAGAAGTGGATTACgag             >  1:528798/1‑149 (MQ=255)
                                                                                                                                           aagaagTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTGAAGAAGTGGATTACGAGaaa          <  1:438080/149‑1 (MQ=255)
                                                                                                                                             gaagTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTGAAGAAGTGGATTACGAGAAAtc        >  2:886993/1‑149 (MQ=255)
                                                                                                                                               agTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTGAAGAAGTGGATTACGAGAAAtctc      <  1:2522353/149‑1 (MQ=255)
                                                                                                                                                  cGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTGAAGAAGTGGATTACGAGAAATCTCGTc   <  1:3493310/149‑1 (MQ=255)
                                                                                                                                                  cGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTGAAGAAGTGGATTACGAGAAATCTCGTc   <  2:4072705/149‑1 (MQ=255)
                                                                                                                                                  cGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTGAAGAAGTGGATTACGAGAAATCTCGTc   <  1:1730400/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTGAAGAAGTGGATTACGAGAAATCTCGTCT  >  CP009273/4167848‑4168143

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: