Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,958,333 C→A 18.2% E220* (GAA→TAA)  rhlB ← ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,958,3330CA18.2% 845.9 / 142.4 295E220* (GAA→TAA) rhlBATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base C (134/107);  new base A (27/27);  total (161/134)
Fisher's exact test for biased strand distribution p-value = 5.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 296 total aligned reads displayed.

ACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTAC  >  CP009273/3958184‑3958481
                                                                                                                                                     |                                                                                                                                                    
 caACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCa                                                                                                                                                     <  1:3361053/150‑1 (MQ=255)
   aCGGTTTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAat                                                                                                                                                   <  1:4357258/150‑1 (MQ=255)
    cGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAata                                                                                                                                                  <  2:3829976/150‑1 (MQ=255)
     ggTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAatat                                                                                                                                                 >  2:932498/1‑150 (MQ=255)
     ggTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAatat                                                                                                                                                 >  1:1800283/1‑150 (MQ=255)
     ggTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAatat                                                                                                                                                 >  1:2633125/1‑150 (MQ=255)
      gtgtTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAatata                                                                                                                                                >  2:3210706/1‑150 (MQ=255)
       tgtTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAatatat                                                                                                                                               >  2:1023754/1‑150 (MQ=255)
        gtTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTAAATATAtt                                                                                                                                              <  2:1634705/150‑1 (MQ=255)
            ggTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGc                                                                                                                                          <  2:3514029/150‑1 (MQ=255)
             gtgtTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTAAATATATTCGGCa                                                                                                                                         <  2:2331448/150‑1 (MQ=255)
              tgtTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCAt                                                                                                                                        >  1:647627/1‑150 (MQ=255)
              tgtTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCAt                                                                                                                                        <  1:1945989/150‑1 (MQ=255)
                 tGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGt                                                                                                                                     <  1:2465450/150‑1 (MQ=255)
                  ggCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGtt                                                                                                                                    <  2:2853826/150‑1 (MQ=255)
                  ggCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGtt                                                                                                                                    <  1:3916881/150‑1 (MQ=255)
                    cGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCa                                                                                                                                  >  1:1503036/1‑150 (MQ=255)
                      aaaaTAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATc                                                                                                                                >  2:1747149/1‑150 (MQ=255)
                       aaaTAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCt                                                                                                                               >  1:829805/1‑150 (MQ=255)
                       aaaTAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCt                                                                                                                               >  1:1680961/1‑150 (MQ=255)
                       aaaTAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCt                                                                                                                               >  1:3807505/1‑150 (MQ=255)
                        aataatCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTg                                                                                                                              <  1:2313783/150‑1 (MQ=255)
                            atCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCg                                                                                                                          <  2:3041436/150‑1 (MQ=255)
                            atCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCg                                                                                                                          <  2:4430236/150‑1 (MQ=255)
                               gCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAAc                                                                                                                       >  2:4312104/1‑150 (MQ=255)
                               gCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTGCAATATATTCGGCATTGTTCATCTGCTCGAAc                                                                                                                       >  2:1297753/1‑150 (MQ=255)
                                  ctgtctgGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGcc                                                                                                                    >  2:3284483/1‑150 (MQ=255)
                                   tgtctgGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCa                                                                                                                   <  2:3162824/150‑1 (MQ=255)
                                   tgtctgGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTAAATATATTCGGCATTGTTCATCTGCTCGAACGCCa                                                                                                                   <  2:1722019/150‑1 (MQ=255)
                                    gtctgGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAg                                                                                                                  <  2:2987046/150‑1 (MQ=255)
                                         gCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCAc                                                                                                             >  1:1248417/1‑150 (MQ=255)
                                            aCTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCTGCATTGTTCATCTGCTCGAACGCCAGTTCACGTa                                                                                                          >  1:4850411/1‑150 (MQ=255)
                                               cttcttCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTAccc                                                                                                       <  1:3507439/150‑1 (MQ=255)
                                                ttcttcGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCg                                                                                                      >  1:1592335/1‑150 (MQ=255)
                                                ttcttcGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTAAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCg                                                                                                      <  1:4408567/150‑1 (MQ=255)
                                                            cGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGt                                                                                          >  2:2438312/1‑150 (MQ=255)
                                                            cGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGt                                                                                          <  1:3246018/150‑1 (MQ=255)
                                                            cGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTGTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGt                                                                                          >  2:3490978/1‑150 (MQ=255)
                                                             gTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCGCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTg                                                                                         >  2:4252104/1‑150 (MQ=255)
                                                                 gcagcaAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCgg                                                                                     >  1:4263725/1‑150 (MQ=255)
                                                                  cagcaAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGga                                                                                    <  1:408204/150‑1 (MQ=255)
                                                                   agcaAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGgag                                                                                   <  2:3384831/150‑1 (MQ=255)
                                                                   agcaAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGgag                                                                                   >  1:377219/1‑150 (MQ=255)
                                                                    gcaAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGgaga                                                                                  >  2:1065328/1‑150 (MQ=255)
                                                                      aaaCGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAAc                                                                                <  2:293590/150‑1 (MQ=255)
                                                                       aaCGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACa                                                                               >  2:2566345/1‑150 (MQ=255)
                                                                            tttttttCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATg                                                                          >  1:2294100/2‑150 (MQ=255)
                                                                              tttttCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGtt                                                                        <  1:3975000/150‑1 (MQ=255)
                                                                                tttCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGa                                                                      >  1:2529828/1‑150 (MQ=255)
                                                                                 ttcttcGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCTCTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAg                                                                     <  1:2177020/150‑1 (MQ=255)
                                                                                  tcttcGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTAAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAgg                                                                    >  2:1806780/1‑150 (MQ=255)
                                                                                   cttcGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGgc                                                                   <  2:1692119/150‑1 (MQ=255)
                                                                                      cGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCt                                                                <  2:4003125/150‑1 (MQ=255)
                                                                                       gTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTAAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTg                                                               >  1:4019919/1‑150 (MQ=255)
                                                                                          aGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGtt                                                            <  1:4118502/150‑1 (MQ=255)
                                                                                            aaGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTAAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTg                                                          <  1:1008115/150‑1 (MQ=255)
                                                                                             aGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCTGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGGTGGTTTGc                                                         >  2:2483797/1‑150 (MQ=255)
                                                                                             aGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGc                                                         >  1:2825992/1‑150 (MQ=255)
                                                                                                   gAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAggcgg                                                   >  1:4621511/1‑150 (MQ=255)
                                                                                                    aaaaGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTAAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAggcggc                                                  >  1:92247/1‑150 (MQ=255)
                                                                                                     aaaGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCa                                                 <  2:1045634/150‑1 (MQ=255)
                                                                                                     aaaGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCa                                                 >  1:3284875/1‑150 (MQ=255)
                                                                                                       aGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATa                                               >  1:4528443/1‑150 (MQ=255)
                                                                                                        gCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATAc                                              <  2:1509872/150‑1 (MQ=255)
                                                                                                         cgcgtcttTAAGACGGTGGCCCGTTTTCTGGTCCGGGTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATAcg                                             <  1:3741979/146‑1 (MQ=255)
                                                                                                          tcttctTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATAcgg                                            >  2:202437/1‑150 (MQ=255)
                                                                                                          tcttctTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATAcgg                                            <  1:3453512/150‑1 (MQ=255)
                                                                                                          tcttctTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATAcgg                                            >  2:2991891/1‑150 (MQ=255)
                                                                                                          tcttctTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTAAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATAcgg                                            >  1:2130415/1‑150 (MQ=255)
                                                                                                           cttcttTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATAcggc                                           >  2:714414/1‑150 (MQ=255)
                                                                                                             tcttTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATAcggcgg                                         >  1:1568852/1‑150 (MQ=255)
                                                                                                             tcttTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTAAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATAcggcgg                                         <  2:2843675/150‑1 (MQ=255)
                                                                                                              cttTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGa                                        >  2:4778865/1‑150 (MQ=255)
                                                                                                                ttAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAAc                                      >  1:1989633/1‑150 (MQ=255)
                                                                                                                 tAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAAca                                     >  1:3159520/1‑150 (MQ=255)
                                                                                                                 tAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTAAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAAca                                     >  2:4674317/1‑150 (MQ=255)
                                                                                                                  aaTACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTAAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAAcag                                    >  1:2939045/1‑150 (MQ=255)
                                                                                                                  aaTACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTAAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAAcag                                    >  1:1402182/1‑150 (MQ=255)
                                                                                                                    tACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTAAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAAcagcc                                  >  2:1418686/1‑150 (MQ=255)
                                                                                                                     aCGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAAcagcca                                 <  2:1157128/150‑1 (MQ=255)
                                                                                                                      cGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAAcagccag                                >  2:4811972/1‑150 (MQ=255)
                                                                                                                        gTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCg                              <  1:1897652/150‑1 (MQ=255)
                                                                                                                             ccGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGatat                         <  2:1592042/150‑1 (MQ=255)
                                                                                                                             ccGTTTTCTGTTCCGGTTCCACTTAAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGatat                         <  2:461118/150‑1 (MQ=255)
                                                                                                                               gTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCt                       <  2:3589625/150‑1 (MQ=255)
                                                                                                                                ttttCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCtt                      <  2:1373296/150‑1 (MQ=255)
                                                                                                                                ttttCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCtt                      <  1:3052124/150‑1 (MQ=255)
                                                                                                                                ttttCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCtt                      <  2:2895893/150‑1 (MQ=255)
                                                                                                                                ttttCTGTTCCGGTTCCACTTAAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCtt                      >  1:1519437/1‑150 (MQ=255)
                                                                                                                                 tttCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCttt                     <  1:1747149/150‑1 (MQ=255)
                                                                                                                                   tCTGTTCCGGTTCCACTTAAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTtaa                   <  2:2116155/150‑1 (MQ=255)
                                                                                                                                   ccTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTtaa                   <  1:2974747/149‑1 (MQ=255)
                                                                                                                                         ccGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGc             >  1:1838713/1‑150 (MQ=255)
                                                                                                                                           ggTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGccc           >  2:1207175/1‑150 (MQ=255)
                                                                                                                                            gTTCCACTTCAATATATTCGGCAGTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCa          <  2:4408932/150‑1 (MQ=255)
                                                                                                                                               ccACTTAAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGAt       <  2:1912517/150‑1 (MQ=255)
                                                                                                                                                cACTTAAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATc      <  2:2057870/150‑1 (MQ=255)
                                                                                                                                                  cTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGt    >  1:3634946/1‑150 (MQ=255)
                                                                                                                                                   ttCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTa   <  2:2389829/150‑1 (MQ=255)
                                                                                                                                                    tCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTAc  >  2:4546763/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTAC  >  CP009273/3958184‑3958481

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