Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,596,243 T→A 45.6% intergenic (‑25/+243) leuQ ← / ← rsmC tRNA‑Leu/16S rRNA m(2)G1207 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,596,2430TA45.6% 117.4 / inf 329intergenic (‑25/+243)leuQ/rsmCtRNA‑Leu/16S rRNA m(2)G1207 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (99/80);  new base A (90/60);  total (189/140)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

Only 100 of 333 total aligned reads displayed.

CACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGC  >  CP009273/4596095‑4596392
                                                                                                                                                    |                                                                                                                                                     
     tCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTAcc                                                                                                                                                 >  2:562874/1‑150 (MQ=255)
      cGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCa                                                                                                                                                >  1:3832127/1‑150 (MQ=255)
        cccATTGCTAGCAATATCATCGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACg                                                                                                                                              <  2:1868824/150‑1 (MQ=255)
        cccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACg                                                                                                                                              >  2:2103169/1‑150 (MQ=255)
        cccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACg                                                                                                                                              >  2:718037/1‑150 (MQ=255)
        cccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACg                                                                                                                                              >  1:3189610/1‑150 (MQ=255)
         ccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACgg                                                                                                                                             >  1:762090/1‑150 (MQ=255)
         ccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACgg                                                                                                                                             >  2:3064043/1‑150 (MQ=255)
           agagCGAGCAATATCAACTTGGTTTTTGGTGCGTGGGGGTGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAg                                                                                                                                           <  1:1271710/147‑1 (MQ=255)
               cGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGgcgc                                                                                                                                       >  2:2444035/1‑150 (MQ=255)
                tagaactatCAAAGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCa                                                                                                                                      <  1:5117770/144‑1 (MQ=255)
                  gCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTAcc                                                                                                                                                 <  2:3242458/137‑1 (MQ=255)
                  gCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTAcc                                                                                                                                                 >  1:3242458/1‑137 (MQ=255)
                    aaTATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCt                                                                                                                                  <  1:3547867/150‑1 (MQ=255)
                      tatCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAg                                                                                                                                >  1:3003907/1‑150 (MQ=255)
                         cAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTgg                                                                                                                             >  1:3635683/1‑150 (MQ=255)
                          atcggggTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGt                                                                                                                            <  2:4900372/145‑1 (MQ=255)
                          aaCGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGt                                                                                                                            >  2:266634/1‑150 (MQ=255)
                            cgctgcgtttGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGttt                                                                                                                          <  1:2845562/143‑1 (MQ=255)
                             tagtttattGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGtttt                                                                                                                         <  2:4773999/143‑1 (MQ=255)
                              cgtgttttGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTc                                                                                                                        <  2:2174923/146‑1 (MQ=255)
                               gatatttGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCa                                                                                                                       <  1:2264509/146‑1 (MQ=255)
                               gataattGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCa                                                                                                                       <  1:1815479/145‑1 (MQ=255)
                                 tttttGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGc                                                                                                                     >  1:2809685/1‑150 (MQ=255)
                                    ttGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCttt                                                                                                                  <  1:687423/150‑1 (MQ=255)
                                    ttGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCttt                                                                                                                  <  1:1708576/150‑1 (MQ=255)
                                    ttGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCtgt                                                                                                                  >  1:3860435/1‑148 (MQ=255)
                                     tGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCtttt                                                                                                                 <  2:2338263/150‑1 (MQ=255)
                                     tGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCtttt                                                                                                                 <  1:1773545/150‑1 (MQ=255)
                                          cGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTg                                                                                                                                                      >  2:1364290/1‑108 (MQ=255)
                                          cGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCa                                                                                                            <  1:11250/150‑1 (MQ=255)
                                            aGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAAt                                                                                                          <  1:2977776/150‑1 (MQ=255)
                                             ggggggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATa                                                                                                         <  1:1383632/150‑1 (MQ=255)
                                             ggggggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATa                                                                                                         >  2:2075817/1‑150 (MQ=255)
                                               ggggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGt                                                                                                       >  1:3982065/1‑150 (MQ=255)
                                                  gggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTaa                                                                                                    <  1:453758/150‑1 (MQ=255)
                                                  gggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTaa                                                                                                    >  1:4981293/1‑150 (MQ=255)
                                                       ttGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATc                                                                                               >  1:2561663/1‑150 (MQ=255)
                                                       ttGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATc                                                                                               >  1:1591947/1‑150 (MQ=255)
                                                               ccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAgg                                                                                       <  1:100627/150‑1 (MQ=255)
                                                                     ccGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAAt                                                                                 >  2:4435649/1‑150 (MQ=255)
                                                                      cGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAAtt                                                                                <  1:4896710/150‑1 (MQ=255)
                                                                      cGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAAtt                                                                                >  1:750961/1‑150 (MQ=255)
                                                                         aaGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCt                                                                             <  1:1986161/150‑1 (MQ=255)
                                                                         aaGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCt                                                                             <  1:4613783/150‑1 (MQ=255)
                                                                            gACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGa                                                                          <  1:4880309/150‑1 (MQ=255)
                                                                             acacTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGaa                                                                         <  1:1653537/150‑1 (MQ=255)
                                                                              cacTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGaaa                                                                        <  1:4804817/150‑1 (MQ=255)
                                                                                cTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGaaaaa                                                                      >  2:3006428/1‑150 (MQ=255)
                                                                                  aaCACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATg                                                                    <  1:2088312/150‑1 (MQ=255)
                                                                                  aaCACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATg                                                                    >  1:2390416/1‑150 (MQ=255)
                                                                                   acacCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGt                                                                   >  1:1949440/1‑150 (MQ=255)
                                                                                      ccTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCa                                                                >  2:3983150/1‑150 (MQ=255)
                                                                                       cTGAAGCTAGCGCGTCTACCAAGTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCAt                                                               <  2:3877256/150‑1 (MQ=255)
                                                                                         gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATACTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCa                                                             >  1:1714555/1‑150 (MQ=255)
                                                                                          aaGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAg                                                            <  1:3375937/150‑1 (MQ=255)
                                                                                          aaGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAg                                                            <  1:2622133/150‑1 (MQ=255)
                                                                                          aaGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAg                                                            >  1:3479365/1‑150 (MQ=255)
                                                                                           aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTGTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAgg                                                           >  2:1081358/1‑150 (MQ=255)
                                                                                           aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAgg                                                           >  1:3112998/1‑150 (MQ=255)
                                                                                            gCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGa                                                          >  2:1994891/1‑150 (MQ=255)
                                                                                              tAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGaaa                                                        <  2:4528538/150‑1 (MQ=255)
                                                                                                gcgcGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACt                                                      >  1:932087/1‑150 (MQ=255)
                                                                                                   cgTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGc                                                   >  2:1137296/1‑150 (MQ=255)
                                                                                                   cgTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCATTTCGTCAATAGTTAATTATCGACAGAGGTGTAAATGCTGGAAAAATGTCCATCAGGAAACTAGc                                                   >  2:2987551/1‑150 (MQ=255)
                                                                                                      cTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTg                                                >  2:650527/1‑150 (MQ=255)
                                                                                                        aCCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCa                                              >  2:3950839/1‑150 (MQ=255)
                                                                                                        aCCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCa                                              <  1:4180819/150‑1 (MQ=255)
                                                                                                         ccAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAg                                             <  1:4336293/150‑1 (MQ=255)
                                                                                                           aaTTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGt                                           >  2:2709522/1‑150 (MQ=255)
                                                                                                           aaTTCCGCCACCGTCGCATACCATCAATTCTTAAAAAGAATAGCTCCCACGGAGGCGCATTCGAGTGGTTTTCAGCTTTTCGTCCATAGTTCATTATCGACAGAGGTGTAATTGCTGGAAAAAGGTCCATCAGGAAACTAGCGTGCAGGt                                           >  1:581157/1‑150 (MQ=255)
                                                                                                             ttCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGttt                                         >  2:4657075/1‑150 (MQ=255)
                                                                                                              tCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTg                                        <  2:1003325/150‑1 (MQ=255)
                                                                                                                 gCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTa                                     <  2:885192/150‑1 (MQ=255)
                                                                                                                  ccACCTTAGCATACCATCAATTCTTAAAAAGTATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATCGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTat                                    <  2:1963650/150‑1 (MQ=255)
                                                                                                                   cACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTatg                                   >  2:3061072/1‑150 (MQ=255)
                                                                                                                       ttCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTatgcatg                               <  2:933264/150‑1 (MQ=255)
                                                                                                                          gCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTCATTATCGACCGAGGTGTAATTGCTGGACAAATGTCCATCAGGAACCTAGCGTGCAGGGGGGGGCTGCATGCgg                            >  2:2749360/1‑150 (MQ=255)
                                                                                                                          gCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCgg                            <  1:1562574/150‑1 (MQ=255)
                                                                                                                           cATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCggg                           >  1:834573/1‑150 (MQ=255)
                                                                                                                           cATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCggg                           >  2:154832/1‑150 (MQ=255)
                                                                                                                              aCCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGc                        <  2:3159375/150‑1 (MQ=255)
                                                                                                                               ccATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGca                       <  2:5096636/150‑1 (MQ=255)
                                                                                                                               ccATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGca                       >  1:1540045/1‑150 (MQ=255)
                                                                                                                                cATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGGGTGGTATGCATGCGGGGGcag                      >  2:4261149/1‑150 (MQ=255)
                                                                                                                                  tCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcagat                    <  2:1332999/150‑1 (MQ=255)
                                                                                                                                   cAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcagatg                   >  2:165583/1‑150 (MQ=255)
                                                                                                                                     aTTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcagatgcc                 >  2:5042976/1‑150 (MQ=255)
                                                                                                                                      ttCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTAAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcagatgcca                >  1:4757057/1‑150 (MQ=255)
                                                                                                                                      ttCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcagatgcca                <  2:806998/150‑1 (MQ=255)
                                                                                                                                         ttAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcagatgccagat             <  1:4237536/150‑1 (MQ=255)
                                                                                                                                             aaaGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGa         >  2:2253019/1‑150 (MQ=255)
                                                                                                                                             aaaGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGa         >  1:1455496/1‑150 (MQ=255)
                                                                                                                                              aaGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGAc        >  2:2057421/1‑150 (MQ=255)
                                                                                                                                               aGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACg       >  1:4799669/1‑150 (MQ=255)
                                                                                                                                               aGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACg       <  2:750961/150‑1 (MQ=255)
                                                                                                                                                 aaTTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGAGGCGACGCt     >  2:3805232/1‑150 (MQ=255)
                                                                                                                                                 aaTAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCt     <  2:3585990/150‑1 (MQ=255)
                                                                                                                                                  aTTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTg    >  1:412166/1‑150 (MQ=255)
                                                                                                                                                   ttGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTgg   >  1:874164/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGC  >  CP009273/4596095‑4596392

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: