Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,884,978 C→A 5.4% D183Y (GAT→TAT)  tsaB ← tRNA(ANN) t(6)A37 threonylcarbamoyladenosine modification protein; binding partner and protease for TsaD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,884,9780CA 5.4% 867.1 / 15.8 229D183Y (GAT→TAT) tsaBtRNA(ANN) t(6)A37 threonylcarbamoyladenosine modification protein; binding partner and protease for TsaD
Reads supporting (aligned to +/- strand):  ref base C (113/103);  new base A (7/6);  total (120/109)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 236 total aligned reads displayed.

CATTCTTTGCCCGGAAGTTTCTTCCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATA  >  CP009273/1884830‑1885126
                                                                                                                                                    |                                                                                                                                                    
cattcTTTGCCCGGAAGTTTCTTCCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATc                                                                                                                                                      >  2:3853644/1‑149 (MQ=255)
  ttcTTTGCCCGGAAGTTTCTTCAATTCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGAAACGCAATCGGCAGCATATCTTAAGAAGCAGGAAGTAACACTTCGCCATcgc                                                                                                                                                    >  2:362224/1‑149 (MQ=255)
    cTTTGCCCGGAAGTTTCTTCCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCgcgc                                                                                                                                                  <  1:1553175/149‑1 (MQ=255)
    cTTTGCCCGGAAGTTTCTTCCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCgcgc                                                                                                                                                  <  1:1684533/149‑1 (MQ=255)
      ttGCCCGGAAGTTTCTTCCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCaa                                                                                                                                                <  1:2813947/149‑1 (MQ=255)
         cccGGAAGTTTCTTCCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAAc                                                                                                                                             <  2:4061640/149‑1 (MQ=255)
         cccGGAAGTTTCTTCCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAAc                                                                                                                                             >  2:345361/1‑149 (MQ=255)
           cGGAAGTTTCTTCCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCa                                                                                                                                           >  2:3773989/1‑149 (MQ=255)
           cGGAAGTTTCTTCCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCa                                                                                                                                           >  2:667486/1‑149 (MQ=255)
           cGGAAGTTTCTTCCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCa                                                                                                                                           >  2:2035950/1‑149 (MQ=255)
           cGGAAGTTTCTTCCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCAACGCGCAAAACCa                                                                                                                                           >  2:1179722/1‑149 (MQ=255)
             gAAGTTTCTTCCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGc                                                                                                                                         >  2:3953741/1‑149 (MQ=255)
               aGTTTCTTCCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGccc                                                                                                                                       >  2:2566586/1‑149 (MQ=255)
                 tttCTTCCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGc                                                                                                                                     <  1:1955672/149‑1 (MQ=255)
                   tcttcCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATAGCGCAAAACCAGCCCGctc                                                                                                                                   >  1:1442481/1‑149 (MQ=255)
                        cATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCttt                                                                                                                              <  2:796081/149‑1 (MQ=255)
                        cATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCttt                                                                                                                              <  2:2854755/149‑1 (MQ=255)
                          tGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTAc                                                                                                                            >  1:4026672/1‑149 (MQ=255)
                           gCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTAcc                                                                                                                           >  1:540448/1‑149 (MQ=255)
                             gACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCga                                                                                                                         >  1:977899/1‑149 (MQ=255)
                             gACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCga                                                                                                                         <  2:24609/149‑1 (MQ=255)
                             gACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCga                                                                                                                         >  1:635672/1‑149 (MQ=255)
                               cGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATAGCGCAAAACCAGCCCGCTCTCTTTACCgaga                                                                                                                       <  2:2633761/149‑1 (MQ=255)
                                 ttgttACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATc                                                                                                                     <  1:2269195/149‑1 (MQ=255)
                                 ttgttACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATc                                                                                                                     >  2:2838444/1‑149 (MQ=255)
                                 ttgttACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATc                                                                                                                     <  1:179786/149‑1 (MQ=255)
                                    ttACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCgg                                                                                                                  >  2:321624/1‑149 (MQ=255)
                                      aCGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGcc                                                                                                                >  2:2690804/1‑149 (MQ=255)
                                           aataaaCCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCt                                                                                                           >  1:2643252/1‑149 (MQ=255)
                                              aaaCCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGc                                                                                                        >  2:3523292/1‑149 (MQ=255)
                                              aaaCCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGc                                                                                                        >  2:2889554/1‑149 (MQ=255)
                                               aaCCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGcc                                                                                                       <  1:1723595/149‑1 (MQ=255)
                                                   ggTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGcc                                                                                                   >  2:556915/1‑149 (MQ=255)
                                                   ggTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGcc                                                                                                   <  2:1677324/149‑1 (MQ=255)
                                                        ggCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATAGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTcc                                                                                              >  2:668626/1‑149 (MQ=255)
                                                         gCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTccc                                                                                             <  1:1332074/149‑1 (MQ=255)
                                                         gCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTccc                                                                                             >  1:2220127/1‑149 (MQ=255)
                                                               tCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGt                                                                                       <  2:3688393/149‑1 (MQ=255)
                                                                  accgccaccgTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCAc                                                                                    >  2:2369052/1‑149 (MQ=255)
                                                                        accgTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTc                                                                              >  1:316660/1‑149 (MQ=255)
                                                                        accgTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTc                                                                              >  1:1633724/1‑149 (MQ=255)
                                                                        accgTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTc                                                                              <  1:1019708/149‑1 (MQ=255)
                                                                               tACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGaa                                                                       <  1:159722/149‑1 (MQ=255)
                                                                               tACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGaa                                                                       >  1:1541920/1‑149 (MQ=255)
                                                                                aCCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGaaa                                                                      <  2:3810701/149‑1 (MQ=255)
                                                                                  ccTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGc                                                                    <  2:1218518/149‑1 (MQ=255)
                                                                                         aaaCATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCa                                                             >  1:3142832/1‑149 (MQ=255)
                                                                                         aaaCATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCa                                                             >  1:2138317/1‑149 (MQ=255)
                                                                                           aCATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTATCACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGGAtt                                                           <  2:3897915/149‑1 (MQ=255)
                                                                                            cATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATAGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTc                                                          <  2:1100229/149‑1 (MQ=255)
                                                                                              tCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGt                                                        <  2:2084102/149‑1 (MQ=255)
                                                                                               cTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGtt                                                       <  1:2269707/149‑1 (MQ=255)
                                                                                                 gACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCa                                                     <  2:830047/149‑1 (MQ=255)
                                                                                                 gACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCa                                                     >  1:3739412/1‑149 (MQ=255)
                                                                                                   cacGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATg                                                   <  1:3906299/149‑1 (MQ=255)
                                                                                                    acGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATgg                                                  >  2:3748265/1‑149 (MQ=255)
                                                                                                    acGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATgg                                                  >  1:3239092/1‑149 (MQ=255)
                                                                                                      gCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGAc                                                <  1:929689/149‑1 (MQ=255)
                                                                                                       cAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACg                                               >  1:3667931/1‑149 (MQ=255)
                                                                                                       cAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACg                                               >  1:3668107/1‑149 (MQ=255)
                                                                                                        aaTCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGa                                              >  2:3036818/1‑149 (MQ=255)
                                                                                                          tCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGAtt                                            <  1:345361/149‑1 (MQ=255)
                                                                                                          tCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGAtt                                            >  2:3825314/1‑149 (MQ=255)
                                                                                                          tCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGAtt                                            >  1:3657865/1‑149 (MQ=255)
                                                                                                           cGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGAttt                                           <  2:1592734/149‑1 (MQ=255)
                                                                                                             gcagcaTATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCg                                         >  1:3822431/1‑149 (MQ=255)
                                                                                                                gcaTATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGt                                      >  1:3367374/1‑149 (MQ=255)
                                                                                                                  atatCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGttt                                    >  1:1141238/1‑149 (MQ=255)
                                                                                                                   tatCTTCAGCAGCAGGCCGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTg                                   >  2:1727391/1‑149 (MQ=255)
                                                                                                                    atCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGa                                  <  2:1531059/149‑1 (MQ=255)
                                                                                                                     tCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAg                                 <  1:3746301/149‑1 (MQ=255)
                                                                                                                     tCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAg                                 <  1:3746239/149‑1 (MQ=255)
                                                                                                                     tCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAg                                 >  1:3548493/1‑149 (MQ=255)
                                                                                                                       ttCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTa                               >  1:3950230/1‑149 (MQ=255)
                                                                                                                        tCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTAc                              <  2:1927438/149‑1 (MQ=255)
                                                                                                                        tCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTAc                              <  2:3461032/149‑1 (MQ=255)
                                                                                                                            cagcagGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCt                          <  1:1355980/149‑1 (MQ=255)
                                                                                                                             agcagGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCtt                         <  1:862363/149‑1 (MQ=255)
                                                                                                                             agcagGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCtt                         <  1:2116517/149‑1 (MQ=255)
                                                                                                                             agcagGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCtt                         >  2:1728690/1‑149 (MQ=255)
                                                                                                                             agcagGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCGCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCtt                         <  2:1467413/149‑1 (MQ=255)
                                                                                                                              gcaggcagTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTc                        <  2:3304220/149‑1 (MQ=255)
                                                                                                                                  gcagTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGtt                    <  2:3386317/149‑1 (MQ=255)
                                                                                                                                  gcagTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGtt                    <  2:2490610/149‑1 (MQ=255)
                                                                                                                                    agTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTttc                  >  1:1102273/1‑149 (MQ=255)
                                                                                                                                    agTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTttc                  >  2:1331099/1‑149 (MQ=255)
                                                                                                                                    agTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTttc                  >  2:1772451/1‑149 (MQ=255)
                                                                                                                                     gTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTttct                 >  2:3152141/1‑149 (MQ=255)
                                                                                                                                      tAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTttctt                >  1:3837571/1‑149 (MQ=255)
                                                                                                                                      tAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTttctt                >  1:2032153/1‑149 (MQ=255)
                                                                                                                                        acacTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGACAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCa              >  2:3598952/1‑149 (MQ=255)
                                                                                                                                        acacTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCa              >  2:500415/1‑149 (MQ=255)
                                                                                                                                        acacTTCGCCATAGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCa              <  2:122020/149‑1 (MQ=255)
                                                                                                                                         cacTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCAc             <  2:635672/149‑1 (MQ=255)
                                                                                                                                            ttCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGt          >  2:1774891/1‑149 (MQ=255)
                                                                                                                                             tCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTg         <  1:2847042/149‑1 (MQ=255)
                                                                                                                                               gCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGcc       >  1:2224321/1‑149 (MQ=255)
                                                                                                                                               gCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATACGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCAGCGTGcc       >  2:3584012/1‑149 (MQ=255)
                                                                                                                                                  aTCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGa    <  1:2939292/149‑1 (MQ=255)
                                                                                                                                                    cgcgCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATa  >  1:4111928/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CATTCTTTGCCCGGAAGTTTCTTCCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATA  >  CP009273/1884830‑1885126

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: