Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,509,102 G→T 100% S57R (AGC→AGA)  dam ← DNA adenine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,509,1020GT98.1% 667.0 / 2.0 211S57R (AGC→AGA) damDNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base T (106/101);  total (108/103)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 212 total aligned reads displayed.

CGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCA  >  CP009273/3508956‑3509245
                                                                                                                                                  |                                                                                                                                               
cggaacggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACNAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTctct                                                                                                                                               <  1:2081946/149‑1 (MQ=255)
 ggaacggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTa                                                                                                                                              >  2:1102903/1‑149 (MQ=255)
  gaacggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTAt                                                                                                                                             >  1:366073/1‑149 (MQ=255)
    acggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTAt                                                                                                                                             <  1:3399530/147‑1 (MQ=255)
    acggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTAt                                                                                                                                             >  2:3399530/1‑147 (MQ=255)
    acggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTAt                                                                                                                                             <  1:1040483/147‑1 (MQ=255)
     tggaTCATGGATTTTTTTGAACTCTTCGAGGTAGTGATAGTAATCCTCGTCGCAATTTGTTTCTGTAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGGATCTTCACAATGTTATAGAGACTGATCAGGTGTCTATTGa                                                                                                                                          <  2:2206184/148‑1 (MQ=255)
      ggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGat                                                                                                                                         >  2:3972060/1‑149 (MQ=255)
      ggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGat                                                                                                                                         <  2:1621695/149‑1 (MQ=255)
        aTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGatat                                                                                                                                       >  2:866399/1‑149 (MQ=255)
        aTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGatat                                                                                                                                       <  1:3267485/149‑1 (MQ=255)
         tCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATc                                                                                                                                      <  1:415270/149‑1 (MQ=255)
         tCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATc                                                                                                                                      <  1:1698418/149‑1 (MQ=255)
            tGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGc                                                                                                                                   >  1:4088576/1‑149 (MQ=255)
              gCTTTTGTTGAACTCTTCGCGGAACTGATATTAAACCTCGGCGCAATTTGTTTCGTGAACAAACAGCTAGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGTTCTCTATTGATATCGGCaa                                                                                                                                 >  1:3618293/1‑149 (MQ=255)
              gCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTTATATCGGCaa                                                                                                                                 <  1:3288346/149‑1 (MQ=255)
                 tttGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGa                                                                                                                              <  1:3167533/149‑1 (MQ=255)
                 tttGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGa                                                                                                                              <  2:1941665/149‑1 (MQ=255)
                 tttGTTGAACTCTTAGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGa                                                                                                                              <  1:440853/149‑1 (MQ=255)
                 tgtgttgAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAATGa                                                                                                                              <  2:624480/147‑1 (MQ=255)
                  ttgttgAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGat                                                                                                                             <  1:2218330/149‑1 (MQ=255)
                      tgAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATaa                                                                                                                         <  1:2788917/149‑1 (MQ=255)
                      tgAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATaa                                                                                                                         <  1:2067083/149‑1 (MQ=255)
                        aaCTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACg                                                                                                                       >  2:401031/1‑149 (MQ=255)
                         aCTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACga                                                                                                                      >  2:3470427/1‑149 (MQ=255)
                           tctTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACgaga                                                                                                                    >  1:1172232/1‑149 (MQ=255)
                               cgcgGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAg                                                                                                                <  1:3450976/149‑1 (MQ=255)
                               cgcgGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAg                                                                                                                <  2:3738726/149‑1 (MQ=255)
                                gcgGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGt                                                                                                               <  1:668447/149‑1 (MQ=255)
                                gcgGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGt                                                                                                               <  1:1918262/149‑1 (MQ=255)
                                 cgGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTc                                                                                                              <  2:366073/149‑1 (MQ=255)
                                 cgGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTc                                                                                                              <  1:2841054/149‑1 (MQ=255)
                                 cgGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTc                                                                                                              <  1:2840070/149‑1 (MQ=255)
                                   gAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCgg                                                                                                            <  2:2473357/149‑1 (MQ=255)
                                    aaCTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGgt                                                                                                           >  2:704786/1‑149 (MQ=255)
                                      cTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGgtgt                                                                                                         <  1:1074447/149‑1 (MQ=255)
                                      cTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGgtgt                                                                                                         <  1:1185513/149‑1 (MQ=255)
                                      cTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGgtgt                                                                                                         >  1:1977528/1‑149 (MQ=255)
                                        gATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTg                                                                                                       <  1:624544/149‑1 (MQ=255)
                                           aGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTgaga                                                                                                    <  2:2319636/149‑1 (MQ=255)
                                           aGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTgaga                                                                                                    <  1:4013793/149‑1 (MQ=255)
                                           aGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTgaga                                                                                                    >  2:1705045/1‑149 (MQ=255)
                                             tAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGaaa                                                                                                  >  2:4036469/1‑149 (MQ=255)
                                                aCCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAacac                                                                                               >  1:1146235/1‑149 (MQ=255)
                                                    cGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGGACGTACTCATCAGTACGCATCTTCACAAGGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGAc                                                                                           <  1:3931107/149‑1 (MQ=255)
                                                    cGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCATGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGAc                                                                                           <  1:1596954/149‑1 (MQ=255)
                                                      gcgcAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGAccc                                                                                         >  1:834210/1‑149 (MQ=255)
                                                      gcgcAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGAccc                                                                                         <  1:2039881/149‑1 (MQ=255)
                                                          aaTTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCa                                                                                     >  1:3105499/1‑149 (MQ=255)
                                                            tttGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACTTACTCATCAGTACGCATCTTCACAATTTTATAGAGACTGATAAGGTCTCTATTGATATCGGCAATGATATAACGAGAAAAGTCGGTTTTGAGAAACACCGACCCGGCAcc                                                                                   <  2:2623382/149‑1 (MQ=255)
                                                             ttGTTTCGGTAACAAACAGCTCGCGTGCGGCCTGTAGGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTTATCAGGTCTCTTTTGATATCGGCAAGTATATAACGATAAAAGTCGGTGTTGATAAACACCGACCCGGCACCt                                                                                  <  2:2186650/149‑1 (MQ=255)
                                                             ttGTTTCGGTAAAAAACAGCTCGCGTGGGGCATGTACGTACTCATCTTTACGAATCTTCAAAATTTTATATAGAGTTATCTGGTCTCTTTTTATATCGGCATGTATATAACGATAAAAGTGGGTGTTGAGAAACACCGACCCGGCACCt                                                                                  <  2:2044536/149‑1 (MQ=255)
                                                             ttGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCt                                                                                  >  2:901564/1‑149 (MQ=255)
                                                                tttCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACa                                                                               >  2:2725586/1‑149 (MQ=255)
                                                                tttCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACa                                                                               >  2:1847389/1‑149 (MQ=255)
                                                                     ggAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAgg                                                                          >  1:262/1‑149 (MQ=255)
                                                                        acaaagaGCTCGCTTGCTGCCTGTAGGTACTCATCAGTACGCATCTTCAATATGTTTTTGATAATGTTCTGTTGTATATTGATATGGGATTGTATTTAACGATAAATGTGTGTGTTTATAAACACCGACCCGGCACCTACAAAAGGCTc                                                                       <  2:3543236/149‑1 (MQ=255)
                                                                               gCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGa                                                                >  1:1530123/1‑149 (MQ=255)
                                                                                  ggcgTGCGGCATGTACGTACACATCAGTACGCATCTTCACAAGGTTATCGAGACTGATCTGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGATAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACAt                                                             <  2:121331/148‑1 (MQ=255)
                                                                                  cgcgTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACAt                                                             <  2:2130926/149‑1 (MQ=255)
                                                                                           ccTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTg                                                    <  2:3740602/149‑1 (MQ=255)
                                                                                              gtacgtacTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGc                                                 >  2:3687561/1‑149 (MQ=255)
                                                                                                      tcatctGTACGAATCTCCACATTGTTATAGCCACTGATCTGGTCGCTATTGATATCGGAAAGTATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCg                                         <  1:3061587/149‑1 (MQ=255)
                                                                                                      tcatcaGTACGCATCTTCACAATGTTATAGAGACTGATCAGTTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGctg                                         <  2:2291573/149‑3 (MQ=255)
                                                                                                      tcatcaGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCg                                         >  1:516285/1‑149 (MQ=255)
                                                                                                        atcaGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGtt                                       <  2:3251625/149‑1 (MQ=255)
                                                                                                        atcaGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGtt                                       >  2:1338892/1‑149 (MQ=255)
                                                                                                        atcaGTACGAATCTTCACAATTTTATAGATACTGATCAGGTATATATTTATATCGGCAAGGATATAACGATAAAAGTCGGTGTTGATAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGtt                                       <  2:1172232/149‑1 (MQ=255)
                                                                                                         tcaGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGTATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGttt                                      <  2:1187330/149‑1 (MQ=255)
                                                                                                           aGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTaa                                    <  2:1686740/149‑1 (MQ=255)
                                                                                                            gTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAat                                   >  1:2404675/1‑149 (MQ=255)
                                                                                                             tACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAata                                  >  2:3244039/1‑149 (MQ=255)
                                                                                                               cGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATatc                                <  2:2075684/149‑1 (MQ=255)
                                                                                                                gCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATatca                               >  1:3821632/1‑149 (MQ=255)
                                                                                                                  aTCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATatcatc                             >  1:933822/1‑149 (MQ=255)
                                                                                                                   tCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATAtcatca                            >  2:3526443/1‑149 (MQ=255)
                                                                                                                    cTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCaa                           >  2:2900210/1‑149 (MQ=255)
                                                                                                                    cTTCAAATTTTTATATAGACTGATCTGGTATTTATTTATATCGTAAAGGATAAAACGATAAAAGTGGGTGTTGACAAAAACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCaa                           <  1:3960028/149‑1 (MQ=255)
                                                                                                                     ttCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAg                          >  2:3307989/1‑149 (MQ=255)
                                                                                                                     ttCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAg                          >  2:586038/1‑149 (MQ=255)
                                                                                                                       cacaATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCa                        >  2:3805206/1‑149 (MQ=255)
                                                                                                                        acaATGTTATAGAGACTGATCCGGTCTCTTTTGATATCGGCATGTATATAACGATAAAAGTGGGTTTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAg                       <  1:1256863/149‑1 (MQ=255)
                                                                                                                        acaATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAg                       >  1:1813404/1‑149 (MQ=255)
                                                                                                                           aTGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAgggg                    >  1:3418892/1‑149 (MQ=255)
                                                                                                                              ttATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGTATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATa                 <  2:3424143/149‑1 (MQ=255)
                                                                                                                              ttATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATa                 >  1:126189/1‑149 (MQ=255)
                                                                                                                               tataGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATAc                >  2:1762167/1‑149 (MQ=255)
                                                                                                                                ataGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACt               >  2:3908067/1‑149 (MQ=255)
                                                                                                                                 taGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACtt              >  1:3249154/1‑149 (MQ=255)
                                                                                                                                 taGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACtt              >  2:1470622/1‑149 (MQ=255)
                                                                                                                                  agagCCTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTCATATCTTCAAGCAGGGGCTCCTTg             >  2:3596735/1‑149 (MQ=255)
                                                                                                                                   gagaCTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGc            <  2:3591999/149‑1 (MQ=255)
                                                                                                                                   gagaCTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGc            <  2:1546697/149‑1 (MQ=255)
                                                                                                                                   gagaCTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCGACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAGTAGCATCAAGCAGGGGATACTTga            >  2:942539/1‑148 (MQ=255)
                                                                                                                                    agaCTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGcc           >  2:591366/1‑149 (MQ=255)
                                                                                                                                     gaCTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTCCAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGCTCCTTGccc          >  2:2844083/1‑149 (MQ=255)
                                                                                                                                       cTGATCAGGTCTCTATTTATATCGGCAAGTATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGccccc        <  2:538932/149‑1 (MQ=255)
                                                                                                                                            cAGGTCTATATTGATCTCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGccc   <  2:3117235/149‑1 (MQ=255)
                                                                                                                                            cAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGccc   >  2:3213800/1‑149 (MQ=255)
                                                                                                                                             aGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCa  >  1:2820372/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
CGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCA  >  CP009273/3508956‑3509245

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