Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,509,102 G→T 89.5% 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,1020GT89.5% 651.2 / 67.6 250S57R (AGC→AGAdamDNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (15/12);  new base T (108/115);  total (123/127)
Fisher's exact test for biased strand distribution p-value = 5.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 252 total aligned reads displayed.

CGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTT  >  CP009273/3508956‑3509248
                                                                                                                                                  |                                                                                                                                                  
  gaacggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTAt                                                                                                                                                <  2:959183/149‑1 (MQ=255)
  gaacggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTAt                                                                                                                                                <  1:4130771/149‑1 (MQ=255)
  gaacggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTAt                                                                                                                                                <  1:1740655/149‑1 (MQ=255)
   aacggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTAtt                                                                                                                                               <  1:121080/149‑1 (MQ=255)
     cggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGa                                                                                                                                             >  2:2800814/1‑149 (MQ=255)
      ggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGat                                                                                                                                            <  2:1330117/149‑1 (MQ=255)
       gaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGata                                                                                                                                           <  1:596182/149‑1 (MQ=255)
        aTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGatat                                                                                                                                          >  2:4311750/1‑149 (MQ=255)
        aTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGatat                                                                                                                                          <  2:1735004/149‑1 (MQ=255)
        aTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGatat                                                                                                                                          <  2:382678/149‑1 (MQ=255)
        aTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGatat                                                                                                                                          <  2:732815/149‑1 (MQ=255)
        aTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGatat                                                                                                                                          >  2:1782887/1‑149 (MQ=255)
        aTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATAGatat                                                                                                                                          >  2:4642581/1‑149 (MQ=255)
         tCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATc                                                                                                                                         <  2:2849793/149‑1 (MQ=255)
         tCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATc                                                                                                                                         >  2:3497315/1‑149 (MQ=255)
         tCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATc                                                                                                                                         >  1:2056202/1‑149 (MQ=255)
         tCCTGGCTTGTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATc                                                                                                                                         <  2:2849852/149‑1 (MQ=255)
          ccTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCg                                                                                                                                        <  1:1472493/149‑1 (MQ=255)
             ggCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCa                                                                                                                                     >  1:2915161/1‑149 (MQ=255)
             ggCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCa                                                                                                                                     >  2:1200457/1‑149 (MQ=255)
                  gtgttgAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGCCTGATCAGGTCTCTATTGATATCGGCAAGGat                                                                                                                                <  2:3391045/148‑1 (MQ=255)
                    gttgAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGatat                                                                                                                              <  1:883742/149‑1 (MQ=255)
                     ttgAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGAtata                                                                                                                             <  2:570703/149‑1 (MQ=255)
                           tctTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACgaga                                                                                                                       >  1:3000089/1‑149 (MQ=255)
                           tctTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACgaga                                                                                                                       <  1:4625151/149‑1 (MQ=255)
                             ttCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGaaa                                                                                                                     <  2:4106288/149‑1 (MQ=255)
                                gcgGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGt                                                                                                                  >  2:2214642/1‑149 (MQ=255)
                                 cgGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTc                                                                                                                 <  2:523309/149‑1 (MQ=255)
                                 cgGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTAAGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTc                                                                                                                 <  2:2594664/149‑1 (MQ=255)
                                  ggAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGCTATAACGAGAAAAGTCg                                                                                                                >  1:3229799/1‑149 (MQ=255)
                                   gAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCgg                                                                                                               <  2:4204255/149‑1 (MQ=255)
                                   gAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCgg                                                                                                               <  2:1336245/149‑1 (MQ=255)
                                   gAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCgg                                                                                                               >  2:522364/1‑149 (MQ=255)
                                    aaCTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCACGGATATAACGAGAAAAGTCGgt                                                                                                              >  1:4466984/1‑149 (MQ=255)
                                    aaCTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGgt                                                                                                              <  1:4380760/149‑1 (MQ=255)
                                     aCTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGgtg                                                                                                             >  1:442884/1‑149 (MQ=255)
                                       tGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGtt                                                                                                           >  1:4500102/1‑149 (MQ=255)
                                       tGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGtt                                                                                                           <  1:4651206/149‑1 (MQ=255)
                                        gATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTg                                                                                                          >  1:357253/1‑149 (MQ=255)
                                         aTAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTga                                                                                                         <  1:3144138/149‑1 (MQ=255)
                                           aGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTgaga                                                                                                       <  1:3959808/149‑1 (MQ=255)
                                            gTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGCTATAACGAGAAAAGTCGGTGTTGAGaa                                                                                                      >  1:1402310/1‑149 (MQ=255)
                                            gTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGaa                                                                                                      >  1:554772/1‑149 (MQ=255)
                                              aaaCCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAac                                                                                                    <  2:1510063/149‑1 (MQ=255)
                                              aaaCCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAac                                                                                                    <  1:2394452/149‑1 (MQ=255)
                                               aaCCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGCCTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAaca                                                                                                   >  1:3930266/1‑149 (MQ=255)
                                               aaCCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAaca                                                                                                   <  1:2465635/149‑1 (MQ=255)
                                               aaCCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAaca                                                                                                   >  1:2261753/1‑149 (MQ=255)
                                                   taggCGCTTTTTGTTTCGGGAACAAACAGCTCTCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATTTTATAGAGTCTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACAGa                                                                                               <  2:4775591/147‑1 (MQ=255)
                                                   tCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGa                                                                                               <  1:1601183/149‑1 (MQ=255)
                                                    cGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGAc                                                                                              <  2:1104483/149‑1 (MQ=255)
                                                          aaTTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCa                                                                                        >  2:3391766/1‑149 (MQ=255)
                                                             ttGTTTCGGTAACAAACAGCTCGCGTGCGGCCTGTACGTAATCATCATTACGAATCTTAACATTGTTATTGAGACTGATCAGGTCGCTTTTGATATCGGCATGGATATACCGTGAAATGTGGGTGTTGAGAAACACCGCCCCGGCACCt                                                                                     <  2:667467/149‑1 (MQ=255)
                                                             ttGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCt                                                                                     >  1:1229289/1‑149 (MQ=255)
                                                                    gggAAAAAACAGGTCGGTTGCGGCGTGTAGGTACTAATCTGTACGCATCTTCAAATTGTTATACAGAGTGATCTGGTGTCTATTGTTTTCGGCAAGGATATAACGATAAACGTCGGTGTTTATAAACACCGACCCGGCACCTACAAAAg                                                                              <  2:1446001/149‑1 (MQ=255)
                                                                     ggAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAgg                                                                             >  1:2675054/1‑149 (MQ=255)
                                                                       aacaaacaGCTCGCGTGCGGCCTGTACGTACTAATAAGTACGAATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCt                                                                           <  2:1949070/149‑1 (MQ=255)
                                                                        acaaacaGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTc                                                                          >  1:3339537/1‑149 (MQ=255)
                                                                            acaGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAATGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAAcc                                                                      <  2:2946735/149‑1 (MQ=255)
                                                                            acaGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAAcc                                                                      >  2:1864002/1‑149 (MQ=255)
                                                                               gCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGa                                                                   >  1:3434370/1‑149 (MQ=255)
                                                                                 tCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACa                                                                 >  2:884115/1‑149 (MQ=255)
                                                                                  cgcgTGCGGCCTGTACGTACTAATCAGTACGCATCTTCACAATGTTATAGATACTTATCAGGTCTCTATTGATATCGGCAAGTATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACAt                                                                <  2:272156/149‑1 (MQ=255)
                                                                                     gTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCg                                                             >  1:2867588/1‑149 (MQ=255)
                                                                                            cTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTgg                                                      >  2:2251383/1‑149 (MQ=255)
                                                                                              gtacgtacTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGc                                                    >  2:4359273/1‑149 (MQ=255)
                                                                                              gtacgtacTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGc                                                    >  2:4358300/1‑149 (MQ=255)
                                                                                               tacttactCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCa                                                   <  2:4154631/149‑1 (MQ=255)
                                                                                                 cgtacTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATAAGGTCTCTATTGATATCGGCAAGGATATAACTAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCaaa                                                 <  2:4690237/149‑1 (MQ=255)
                                                                                                 cgtacTCATCAGTACGCATCTTCACAATGTTATAGAGACTGAACAGTTCTCTATTGATATCGGCAAGTATATAACGAGAAAAGTCGGTGTTTATAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCaaa                                                 <  2:4076393/149‑1 (MQ=255)
                                                                                                    acTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTTATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGc                                              >  1:4414681/1‑149 (MQ=255)
                                                                                                    acTCATCAGTACGCAACTACAAAATGTTATCGAGACGGATCAGGCCGCTATTTATATCGGCAAGTATATAACGAGAAAAGTGGGTTTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGc                                              <  2:590409/149‑1 (MQ=255)
                                                                                                        atcaGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGtt                                          >  1:1241773/1‑149 (MQ=255)
                                                                                                            gTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAat                                      <  2:2085644/149‑1 (MQ=255)
                                                                                                            gTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAat                                      <  2:4655070/149‑1 (MQ=255)
                                                                                                             tACTCAACTTCACAATTTTATAGTGACTGAGCAGGTGGTTTTTGTTAGCGGCACGGATATAACGAGAAACGTCGGTGTTGATAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAata                                     <  2:4407805/149‑1 (MQ=255)
                                                                                                             tACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCACATGCCGTTTAata                                     >  2:2088046/1‑149 (MQ=255)
                                                                                                             tACGAATCTTCACAAGGTTATAGAGACTTATCAGGTCTCTTTTGATATCGGCACGTATATAACGATAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAata                                     <  1:3493410/149‑1 (MQ=255)
                                                                                                                gCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATatca                                  <  2:3896284/149‑1 (MQ=255)
                                                                                                                gCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATatca                                  >  1:3937597/1‑149 (MQ=255)
                                                                                                                 cATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGCGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACAATCGCCCTTGGGCCAATGCCGGTTAATatcat                                 >  1:2165416/1‑149 (MQ=255)
                                                                                                                 aaTCTTCAAAATTTTATAGAGACTTAACTGGACTCTATTTATATCGGAAAGTATATAACGATAAATGTCGGTTTTGACAAACACCGCCCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATatcat                                 <  1:995312/148‑1 (MQ=255)
                                                                                                                 aaTCTACACAATGTTATACACACTGATCCGGTCTATATTGATATCGGAAAGTATATACCGACAAAAGTGGGTTTTGACAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATatcat                                 <  1:3769534/148‑1 (MQ=255)
                                                                                                                        acaATGTTTTAGAGAATTAACAGTTCTGTATTGATATCGGCATTTATATAACGTTAAAATTGGTTGTTTTTAAAAACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATGAAGCAg                          <  2:2841244/149‑1 (MQ=255)
                                                                                                                            tGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGa                      >  1:100774/1‑149 (MQ=255)
                                                                                                                             gTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGAt                     <  2:1954603/149‑1 (MQ=255)
                                                                                                                             gTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGAt                     >  2:639172/1‑149 (MQ=255)
                                                                                                                                 taGAGACTGATCAGGACTCTATTGATATCGGCAAGGATATAACGAGAAAATACGGAGTAGAGAAACACCTACACGGCACCTACAAAAGTCTCAACCAGACATTCGCCCTAGGGCAAATGCCGTTTAATATCATCAAGCAGCGGATACta                 >  1:2005985/1‑148 (MQ=255)
                                                                                                                                  agagACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTg                >  2:3767012/1‑149 (MQ=255)
                                                                                                                                   gagaCTTATCTGGTCTCTATTGATATCGGCATGTATATAACGAGAAAAGGCGGTGTTTAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGc               <  1:3012649/149‑1 (MQ=255)
                                                                                                                                   gagaCTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGc               >  2:1423024/1‑149 (MQ=255)
                                                                                                                                   gagaCTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGc               >  1:964496/1‑149 (MQ=255)
                                                                                                                                       cTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGccccc           >  1:2498910/1‑149 (MQ=255)
                                                                                                                                        tGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCt          >  1:1448789/1‑149 (MQ=255)
                                                                                                                                         gATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTg         >  1:1515736/1‑149 (MQ=255)
                                                                                                                                          aTCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGc        >  1:4244395/1‑149 (MQ=255)
                                                                                                                                          aTCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGc        >  2:771312/1‑149 (MQ=255)
                                                                                                                                            cAGGTCTCTATTGATATCGGCACGTCAATAACGAGAAACGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGccc      <  2:4448907/149‑1 (MQ=255)
                                                                                                                                               gTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACt   >  2:860066/1‑149 (MQ=255)
                                                                                                                                                tCGCTATTGATAGCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACtt  <  2:2961581/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTT  >  CP009273/3508956‑3509248

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: