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.6% 886.6 / ‑0.4 282S57R (AGC→AGAdamDNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base T (136/142);  total (138/144)
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.92e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 286 total aligned reads displayed.

ACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTT  >  CP009273/3508955‑3509248
                                                                                                                                                   |                                                                                                                                                  
 cggaacggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTa                                                                                                                                                 >  2:797926/1‑150 (MQ=255)
 cggaacggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTa                                                                                                                                                 >  1:201738/1‑150 (MQ=255)
  ggaacggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTAt                                                                                                                                                >  2:3750609/1‑150 (MQ=255)
  ggaacggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTAt                                                                                                                                                >  2:3823579/1‑150 (MQ=255)
      cggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGat                                                                                                                                            <  1:3104091/150‑1 (MQ=255)
        gaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGatat                                                                                                                                          <  2:3247129/150‑1 (MQ=255)
          tCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCg                                                                                                                                        <  2:4096435/150‑1 (MQ=255)
           ccTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCgg                                                                                                                                       <  2:4984467/150‑1 (MQ=255)
              ggCTTTTGTTGAACTCTTCGCTGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCaa                                                                                                                                    <  1:1044323/150‑1 (MQ=255)
               gCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAg                                                                                                                                   >  2:4134908/1‑150 (MQ=255)
                cTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAgg                                                                                                                                  <  2:2576184/150‑1 (MQ=255)
                  tttGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGat                                                                                                                                <  1:889753/150‑1 (MQ=255)
                   ttgttgAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGata                                                                                                                               <  1:386535/150‑1 (MQ=255)
                    tgttgAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGatat                                                                                                                              <  2:2246279/150‑1 (MQ=255)
                    tgttgAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGatat                                                                                                                              <  2:1090234/150‑1 (MQ=255)
                    tgttgAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGatat                                                                                                                              <  2:4899432/150‑1 (MQ=255)
                     gttgAACTCGTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGAtata                                                                                                                             <  1:2451213/150‑1 (MQ=255)
                      ttgAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATaa                                                                                                                            <  2:1707101/150‑1 (MQ=255)
                      ttgAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATaa                                                                                                                            <  2:4385660/150‑1 (MQ=255)
                       tgAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAAc                                                                                                                           >  2:4590188/1‑150 (MQ=255)
                       tgAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAAc                                                                                                                           <  2:4752415/150‑1 (MQ=255)
                         aaCTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACga                                                                                                                         >  2:514328/1‑150 (MQ=255)
                          aCTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGCTATAACgag                                                                                                                        >  1:1513266/1‑150 (MQ=255)
                          aCTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACgag                                                                                                                        >  2:4409701/1‑150 (MQ=255)
                             ctTCGCGGAACTGATAGTAAACCTCGGCTCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGTCCTGTAAGTACTCATCAGTACGAATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAATGATATAACGAGaaa                                                                                                                     <  2:3058501/150‑1 (MQ=255)
                                  cgGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCg                                                                                                                >  1:3044333/1‑150 (MQ=255)
                                      aCTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGgtgt                                                                                                            >  1:417947/1‑150 (MQ=255)
                                       cTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGtt                                                                                                           <  2:2863622/150‑1 (MQ=255)
                                          aTAGTATACCTCGTCGCATTTTGTTTCGTGAACAAACAGCTAGCGTGCGGGCTGTTCTTACTCATCTGTACGAATCTTCACAATGTTATAGAGACTGTTCAGGTCTCTATTGATATCGGCAAGGTTATAACGAGATAAGTCGGTGTTgag                                                                                                        <  1:384289/150‑1 (MQ=255)
                                          aTAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTgag                                                                                                        >  2:2474449/1‑150 (MQ=255)
                                          aTAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTgag                                                                                                        <  2:261544/150‑1 (MQ=255)
                                          aTAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTgag                                                                                                        <  1:5151014/150‑1 (MQ=255)
                                            aGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGaa                                                                                                      >  2:2087568/1‑150 (MQ=255)
                                             gTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGCCTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGgaa                                                                                                     >  1:4578445/1‑150 (MQ=255)
                                             gTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGaaa                                                                                                     <  1:54108/150‑1 (MQ=255)
                                                  ccTCGGCTCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTTCTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCg                                                                                                <  1:2603939/150‑1 (MQ=255)
                                                    tCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGAc                                                                                              <  1:2082714/150‑1 (MQ=255)
                                                      ggCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGAccc                                                                                            <  2:1379518/150‑1 (MQ=255)
                                                       gcgcAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCg                                                                                           <  2:1792540/150‑1 (MQ=255)
                                                         gcAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGc                                                                                         <  2:3776125/150‑1 (MQ=255)
                                                           aaTTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCACGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCAc                                                                                       >  1:132303/1‑150 (MQ=255)
                                                           aaTTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCAc                                                                                       >  2:491841/1‑150 (MQ=255)
                                                                     gggAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAgg                                                                             >  1:1405639/1‑150 (MQ=255)
                                                                           aaacaGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAAc                                                                       >  2:3168767/1‑150 (MQ=255)
                                                                            aacaGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAAcc                                                                      >  1:4045630/1‑150 (MQ=255)
                                                                                  tCGCGTGCGGCCTGTACTTACTAATCATTACGCATATTCACAATGTTATAGAGACTGATCTGGTCTCTATTTATATCGGAAAGTATATAACGATAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACAt                                                                <  1:4838913/150‑1 (MQ=255)
                                                                                         cGGCATGTACGTACTAATAAGTACGAATATTCAAAATTTTATATAGAATGATATGTTATCTATTGATATCGTAAAGTATATAACTATAAAAGTGGGTGTTTAGAAACACCGACCCGGCACCTACAAAATGCTCAACCAGACATTCGCCCt                                                         <  2:3599871/150‑1 (MQ=255)
                                                                                          ggCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCtt                                                        >  1:2394540/1‑150 (MQ=255)
                                                                                            cgtttacttactAATCTTTACTAATATTAACATTTTTATTTAGACTTATCTGGTATATTTTTATATCGGAACGTATATAACTATAAATTTGGTTGTTGATAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTgg                                                      <  1:4219019/146‑1 (MQ=255)
                                                                                            ccTGTACGTACTCAACAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGATAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTgg                                                      <  2:588645/150‑1 (MQ=255)
                                                                                             ctttacgtaaTCATCAGTACGAATCTTCACAATGTTATAGAGACTGATCAGGTATCTATTGATATCGCCAAGGATATAACGATAAATGTGGGTGTTGATAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTggg                                                     <  1:5236793/147‑1 (MQ=255)
                                                                                             cTGTACTTACTCATCTTTACGAATCTTCACAATTTTATAGAGACTGATCAGGTCTCTTTTGATATCGGCACGGATATAACGAGAAAAGTCGGTGTTGATAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTggg                                                     <  2:5073718/150‑1 (MQ=255)
                                                                                              tGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAACAGGCTCAACCAGACATTCGCCCTTGGGc                                                    >  2:4299267/1‑150 (MQ=255)
                                                                                               gtacgtacTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCa                                                   >  1:59105/1‑150 (MQ=255)
                                                                                                  cgtacTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAAt                                                >  2:4672722/1‑150 (MQ=255)
                                                                                                   ttacTCATATTTACGAAACTTCAAATTGTTATCGGTACTGATCTGGTCTCTTTTTATAGCGCCAAGTATATAACGTGAAATGTCGGTGTTGACAAACACCCACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATg                                               <  2:2323647/149‑1 (MQ=255)
                                                                                                       tcatcaGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGt                                           >  1:1037610/1‑150 (MQ=255)
                                                                                                        catcaGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGtt                                          >  2:2895751/1‑150 (MQ=255)
                                                                                                         atcaGTCCGCATCTACACAATGTGATAGAGACTGATCCGGTCTCTATTGATATCGGCAAGGATATACCGCGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGttt                                         <  2:1159769/150‑1 (MQ=255)
                                                                                                         atcaGTACGCATCTTCACAATGTTATAGAGACTGATCTGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGcc                                             <  2:1615755/146‑1 (MQ=255)
                                                                                                         atcaGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGcc                                             >  1:1615755/1‑146 (MQ=255)
                                                                                                         atcaGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGttt                                         >  2:1488144/1‑150 (MQ=255)
                                                                                                            aGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTCTTGATATCGGCAAGGATATCACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAat                                      >  1:2347826/1‑150 (MQ=255)
                                                                                                               aCGCATCTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATatc                                   >  2:4815001/1‑150 (MQ=255)
                                                                                                                cGCATCTCCACAATTTTATCGAGACTGAACTGGTCTCTATTGATATCGGCAAGCATATAACGATAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATatca                                  <  1:2948039/150‑1 (MQ=255)
                                                                                                                cGCAACTTCACATTGTTATAGAGACTGATCTGGTCTCTATTGATATCGGCAAGTATATAACGATAAATGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATatca                                  <  1:2179228/150‑1 (MQ=255)
                                                                                                                 gCATCTTCACCTGGTTATAGTGACTGATCCGGTCTATTTTGATATGGGCAAGGTTTTACCGAGAAATGTGGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATatcat                                 <  1:3229042/150‑1 (MQ=255)
                                                                                                                   aTCTTCACAATGTTATTGAGACTGATCTGGTCTCTTTTGATATCGGCAAGTATATACCGAGAAAAGTGGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATAtcatca                               <  2:1930377/150‑1 (MQ=255)
                                                                                                                     cTTCACAATGTTATAGAGACTTATCAGGTCTCTATTGATATATGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGTCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAg                             <  2:701953/150‑1 (MQ=255)
                                                                                                                     cTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAATTCGGTGTTGAGAAACACCGACACGGCACCTACAAAAGGCTCAACCATACATTCGCCCTTGGGCAAATGCAGTTTATTATCATCAAg                             >  1:749614/1‑150 (MQ=255)
                                                                                                                     cTTCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAg                             >  1:2959828/1‑150 (MQ=255)
                                                                                                                      ttaaaattGTTATAGAGACTGAACTGGTATATTTTTATTTCGGAAAGTATATAACTATAAATGTCGGTTTTGATAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGc                            <  1:3750609/143‑1 (MQ=255)
                                                                                                                      ttCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGc                            >  2:111233/1‑150 (MQ=255)
                                                                                                                      ttCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGc                            >  2:248010/1‑150 (MQ=255)
                                                                                                                       tCACAATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCACATGCCGTTTAATATCATCAAGCa                           >  2:2779227/1‑150 (MQ=255)
                                                                                                                       tCACAATGTTATAGAGACTGAACAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCa                           <  1:2689517/150‑1 (MQ=255)
                                                                                                                        cacaATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAg                          >  2:2368470/1‑150 (MQ=255)
                                                                                                                         acaATGTTATAGAGACTGAACCGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAgg                         <  1:4887239/150‑1 (MQ=255)
                                                                                                                          caATGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAggg                        >  1:1951389/1‑150 (MQ=255)
                                                                                                                           aaTGTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCACATGCCGTTTAATATCAGCAAGCAgggg                       >  1:2321065/1‑150 (MQ=255)
                                                                                                                             ttttATCGAGACTTATCTGTTCTATTTTTATATCGGCATTTATATAACTATAAATGTGGGTTTTTACAAACCCCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGAt                     <  1:510876/148‑1 (MQ=255)
                                                                                                                              tttATAGAGACTGATCAGTTCTATTTTTATATCGGAATGTATATAACGATAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATa                    <  1:5193028/149‑1 (MQ=255)
                                                                                                                              gTTTTAGAGACTGATCTGGTCTCTATTGTTATCGGCATGGATATAACGAGAAATGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATa                    <  1:3511925/150‑1 (MQ=255)
                                                                                                                              gTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGATAAATGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATa                    <  1:3440704/150‑1 (MQ=255)
                                                                                                                              gTTATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATa                    >  2:2938664/1‑150 (MQ=255)
                                                                                                                               ttATAGAGACTGATCAGGTCTCTATTGATATCGGCAAGTATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATAc                   <  1:362563/150‑1 (MQ=255)
                                                                                                                                 atatacaCTGATCAGGTATCTATTTATATCGCCATGTATATAACGATAAAAGTCGTTTTTGATAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACtt                 <  1:663235/144‑1 (MQ=255)
                                                                                                                                 atatAGACTTATCTGGTCTCTATTGATATCGGCAAGTATATAACGATAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACtt                 <  1:2126637/150‑1 (MQ=255)
                                                                                                                                 ataGATACTGATCAGGTCTATTTTTATATCGGCAAGTATATAACGATAAAAGTCGGTGTTGAGAAACACCGCCCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACtt                 <  1:1145365/150‑1 (MQ=255)
                                                                                                                                  taGAGACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTg                >  1:467993/1‑150 (MQ=255)
                                                                                                                                   agagACTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGc               >  1:1803042/1‑150 (MQ=255)
                                                                                                                                     agaCTGATCAGGTCTCTATTGATATCGGCAAGGATCTAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCCAATGCCGTGTAATATCATCAAGCAGGGGCTACTTGccc             >  2:867310/1‑150 (MQ=255)
                                                                                                                                     agaCTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGccc             >  1:4466742/1‑150 (MQ=255)
                                                                                                                                     agaCTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCACAACCACACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTAGccc             >  1:2497273/1‑150 (MQ=255)
                                                                                                                                        cTGATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCt          <  2:1248894/150‑1 (MQ=255)
                                                                                                                                          gATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGc        >  2:2012226/1‑150 (MQ=255)
                                                                                                                                          gATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGc        >  2:3435574/1‑150 (MQ=255)
                                                                                                                                          gATCAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTGTAATATCATCAAGCAGGGGATACTTGCCCCCTGc        >  1:1920392/1‑150 (MQ=255)
                                                                                                                                             cAGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCa     >  2:4923704/1‑150 (MQ=255)
                                                                                                                                              aGGTCTCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCAc    >  2:3634444/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTT  >  CP009273/3508955‑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: