Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,084,798 C→A 90.0% L102L (CTC→CTA rsmE → 16S rRNA m(3)U1498 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,084,7980CA90.0% 891.4 / 95.7 341L102L (CTC→CTArsmE16S rRNA m(3)U1498 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (14/20);  new base A (147/160);  total (161/180)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

Only 100 of 345 total aligned reads displayed.

AAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCGGAAATC  >  CP009273/3084650‑3084942
                                                                                                                                                    |                                                                                                                                                
aaaTTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTc                                                                                                                                                  <  1:4363957/149‑1 (MQ=255)
  aTTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAgg                                                                                                                                                >  2:1952012/1‑149 (MQ=255)
  aTTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAgg                                                                                                                                                <  2:3451773/149‑1 (MQ=255)
  aTTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAgg                                                                                                                                                >  2:1140892/1‑149 (MQ=255)
        aGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAg                                                                                                                                          <  1:213415/149‑1 (MQ=255)
         gcgcCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGc                                                                                                                                         >  2:2989379/1‑149 (MQ=255)
         gcgcCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCAAGAAATCGATCGAACTAGGTGTAAGc                                                                                                                                         >  2:3840818/1‑149 (MQ=255)
          cgcCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGcc                                                                                                                                        <  1:3048617/149‑1 (MQ=255)
             cAGCAAAAATAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCa                                                                                                                                     <  1:1159842/149‑1 (MQ=255)
             cAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCa                                                                                                                                     <  2:3613311/149‑1 (MQ=255)
              aGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATATACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCAt                                                                                                                                    >  1:3810540/1‑149 (MQ=255)
               gaaaaaaaaGCGTGGAAGTGAATGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTTATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCAtt                                                                                                                                   <  2:3268662/147‑1 (MQ=255)
               gCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCAtt                                                                                                                                   <  2:1597757/149‑1 (MQ=255)
                 aaaaaaaGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTAc                                                                                                                                 <  1:1372305/149‑1 (MQ=255)
                  aaaaaaGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACg                                                                                                                                >  2:3608813/1‑149 (MQ=255)
                       aGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACt                                                                                                                           >  2:751385/1‑149 (MQ=255)
                        gCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACtt                                                                                                                          <  1:1284177/149‑1 (MQ=255)
                          gTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACtttt                                                                                                                        >  2:1055303/1‑149 (MQ=255)
                             gAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTc                                                                                                                     >  2:3670563/1‑149 (MQ=255)
                              aaGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCt                                                                                                                    <  1:3274960/149‑1 (MQ=255)
                                gTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGACATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGa                                                                                                                  >  1:1730949/1‑149 (MQ=255)
                                gTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGa                                                                                                                  <  2:2873467/149‑1 (MQ=255)
                                gTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGa                                                                                                                  <  1:3537054/149‑1 (MQ=255)
                                  gAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAgc                                                                                                                >  1:4255523/1‑149 (MQ=255)
                                    aGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAgcgc                                                                                                              <  1:3191396/149‑1 (MQ=255)
                                           ggAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTgcggcg                                                                                                       >  2:2136794/1‑149 (MQ=255)
                                            gAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGt                                                                                                      <  1:3238676/149‑1 (MQ=255)
                                              aGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTa                                                                                                    >  2:3259050/1‑149 (MQ=255)
                                                  cAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACt                                                                                                >  2:2356942/1‑149 (MQ=255)
                                                  cAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACt                                                                                                >  2:2711299/1‑149 (MQ=255)
                                                    gATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTgg                                                                                              >  1:1552833/1‑149 (MQ=255)
                                                      tCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGAt                                                                                            <  1:523102/149‑1 (MQ=255)
                                                      tCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGAt                                                                                            <  2:2281226/149‑1 (MQ=255)
                                                             tCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAAc                                                                                     >  2:1718359/1‑149 (MQ=255)
                                                               gcgAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGt                                                                                   <  1:4181456/149‑1 (MQ=255)
                                                                cgAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTc                                                                                  <  2:2094695/149‑1 (MQ=255)
                                                                     ctcCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAAc                                                                             <  1:4776763/149‑1 (MQ=255)
                                                                     ctcCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAAc                                                                             >  2:415387/1‑149 (MQ=255)
                                                                      tcCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACa                                                                            >  1:1006143/1‑149 (MQ=255)
                                                                       ccGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACaa                                                                           <  2:4301564/149‑1 (MQ=255)
                                                                       ccGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACaa                                                                           <  1:1413587/149‑1 (MQ=255)
                                                                       ccGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACaa                                                                           <  1:4591486/149‑1 (MQ=255)
                                                                        cGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACaag                                                                          <  1:1906361/149‑1 (MQ=255)
                                                                          cTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACaagaa                                                                        <  1:1531311/149‑1 (MQ=255)
                                                                          cTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACaagaa                                                                        <  1:3884016/149‑1 (MQ=255)
                                                                               tatTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGCTAGTGAACGTCTGAACAAGAAGCTTc                                                                   >  1:3751463/1‑149 (MQ=255)
                                                                               tatTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTc                                                                   <  2:2580537/149‑1 (MQ=255)
                                                                               tatTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTc                                                                   >  2:2553017/1‑149 (MQ=255)
                                                                                atTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTca                                                                  <  1:4510233/149‑1 (MQ=255)
                                                                                atTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTca                                                                  <  2:634307/149‑1 (MQ=255)
                                                                                 ttCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTcag                                                                 <  2:3550138/149‑1 (MQ=255)
                                                                                 ttCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTcag                                                                 <  2:1330487/149‑1 (MQ=255)
                                                                                     ccTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGt                                                             <  2:2943098/149‑1 (MQ=255)
                                                                                     ccTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGt                                                             >  1:1993705/1‑149 (MQ=255)
                                                                                     ccTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGt                                                             <  1:60887/149‑1 (MQ=255)
                                                                                     ccTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGt                                                             >  1:2221190/1‑149 (MQ=255)
                                                                                       tCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTgg                                                           <  2:4448776/149‑1 (MQ=255)
                                                                                        cGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGc                                                          <  2:3810965/149‑1 (MQ=255)
                                                                                            cAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCagaa                                                      <  2:4552919/149‑1 (MQ=255)
                                                                                             aGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCagaag                                                     >  1:843192/1‑149 (MQ=255)
                                                                                              ggTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCagaaga                                                    <  2:2950211/149‑1 (MQ=255)
                                                                                              ggTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCagaaga                                                    <  2:714265/149‑1 (MQ=255)
                                                                                              ggTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCagaaga                                                    >  2:434936/1‑149 (MQ=255)
                                                                                              ggTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCagaaga                                                    <  1:488665/149‑1 (MQ=255)
                                                                                                tGATGTCGCGTGGTGAAAAAATGGAATGTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGAtt                                                  <  2:533313/149‑1 (MQ=255)
                                                                                                 gATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTg                                                 <  2:2057329/149‑1 (MQ=255)
                                                                                                  aTGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGc                                                >  1:2882007/1‑149 (MQ=255)
                                                                                                  aTGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGc                                                >  2:484529/1‑149 (MQ=255)
                                                                                                   tGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCa                                               >  2:1958473/1‑149 (MQ=255)
                                                                                                    gTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCaa                                              <  1:4043846/149‑1 (MQ=255)
                                                                                                         gtggtgAAAAAATGGAATTGACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATtgct                                         <  1:2076692/149‑1 (MQ=255)
                                                                                                           ggtgAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATtgctgc                                       >  2:3972812/1‑149 (MQ=255)
                                                                                                             tgAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCt                                     <  1:1384332/149‑1 (MQ=255)
                                                                                                             tgAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCt                                     <  2:4688564/149‑1 (MQ=255)
                                                                                                               aaaaaaTGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCtgt                                   <  1:1100956/149‑1 (MQ=255)
                                                                                                                aaaaaTGGAATTTACTATCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCtgtg                                  <  2:3376752/149‑1 (MQ=255)
                                                                                                                aaaaaTGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCtgtg                                  >  1:3639414/1‑149 (MQ=255)
                                                                                                                   aaTGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGc                               >  1:4412069/1‑149 (MQ=255)
                                                                                                                    aTGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCa                              <  1:532103/149‑1 (MQ=255)
                                                                                                                     tGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAg                             >  2:1584357/1‑149 (MQ=255)
                                                                                                                     tGGAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAg                             >  1:3376388/1‑149 (MQ=255)
                                                                                                                       gAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAgtg                           >  1:242348/1‑149 (MQ=255)
                                                                                                                       gAATTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAgtg                           >  1:2810273/1‑149 (MQ=255)
                                                                                                                        aaTTTACTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAgtgt                          >  2:368990/1‑149 (MQ=255)
                                                                                                                              cTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGt                    >  1:3264242/1‑149 (MQ=255)
                                                                                                                              cTATCCAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGt                    >  1:874573/1‑149 (MQ=255)
                                                                                                                                   cAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCg               >  1:705857/1‑149 (MQ=255)
                                                                                                                                   cAGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCg               <  2:826803/149‑1 (MQ=255)
                                                                                                                                    aGAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCgg              >  2:3221253/1‑149 (MQ=255)
                                                                                                                                     gAAATCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCggg             >  2:4391349/1‑149 (MQ=255)
                                                                                                                                      aaaTCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGt            >  2:3109385/1‑149 (MQ=255)
                                                                                                                                      aaaTCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGt            >  1:1805656/1‑149 (MQ=255)
                                                                                                                                       aaTCGATCGAACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTg           >  2:4007693/1‑149 (MQ=255)
                                                                                                                                         tcgatcgaACTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGcc         <  2:146639/149‑1 (MQ=255)
                                                                                                                                         tcgatcgaACTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGcc         >  2:1735318/1‑149 (MQ=255)
                                                                                                                                           gatcgaACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCgg       >  1:3941686/1‑149 (MQ=255)
                                                                                                                                            atcgaACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCGGa      >  1:4576723/1‑149 (MQ=255)
                                                                                                                                            atcgaACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCGGa      >  2:1195713/1‑149 (MQ=255)
                                                                                                                                              cgaACTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCGGaaa    >  1:2039672/1‑149 (MQ=255)
                                                                                                                                                aaCTAGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCGGAAATc  >  1:4729304/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
AAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGCATATTCACCTCGGTCAGGTGATGTCGCGTGGTGAAAAAATGGAATTTACTATCCAGAAATCGATCGAACTCGGTGTAAGCCTCATTACGCCACTTTTTTCTGAGCGCTGCGGCGTTAAACTGGATAGTGAACGTCTGAACAAGAAGCTTCAGCAGTGGCAGAAGATTGCAATTGCTGCCTGTGAGCAGTGTGGTCGTAACCGGGTGCCGGAAATC  >  CP009273/3084650‑3084942

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: