Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,024,703 G→T 91.0% G241G (GGC→GGA rlmI ← 23S rRNA m(5)C1962 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,024,7030GT91.0% 531.0 / 45.3 201G241G (GGC→GGArlmI23S rRNA m(5)C1962 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base G (10/8);  new base T (84/99);  total (94/107)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 201 total aligned reads displayed.

GATCGCGATAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACC  >  CP009273/1024558‑1024847
                                                                                                                                                 |                                                                                                                                                
gatcgcgatAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTcc                                                                                                                                                <  2:536453/148‑1 (MQ=255)
   cgcgatAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCAt                                                                                                                                             <  1:1912701/148‑1 (MQ=255)
     cgatAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCa                                                                                                                                           <  2:1466537/148‑1 (MQ=255)
     cgatAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCa                                                                                                                                           >  2:3163118/1‑148 (MQ=255)
        tAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTg                                                                                                                                        <  2:372213/148‑1 (MQ=255)
         aaGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGc                                                                                                                                       >  1:2614288/1‑148 (MQ=255)
           gTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCg                                                                                                                                     <  1:1497435/148‑1 (MQ=255)
             aCGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGAt                                                                                                                                   >  1:842951/1‑148 (MQ=255)
             aCGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGAt                                                                                                                                   >  1:527052/1‑148 (MQ=255)
              cGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATa                                                                                                                                  >  2:1886744/1‑148 (MQ=255)
              cGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATa                                                                                                                                  <  1:533066/148‑1 (MQ=255)
                  gcaATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGc                                                                                                                              <  2:1749021/148‑1 (MQ=255)
                  gcaATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGc                                                                                                                              >  2:2105961/1‑148 (MQ=255)
                     aTTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGaa                                                                                                                           >  2:905124/1‑148 (MQ=255)
                     aTTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGaa                                                                                                                           <  2:522688/148‑1 (MQ=255)
                        tAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAacc                                                                                                                        >  2:2938058/1‑148 (MQ=255)
                        tAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAacc                                                                                                                        >  2:2035755/1‑148 (MQ=255)
                           aGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAaccacc                                                                                                                     <  2:1166377/148‑1 (MQ=255)
                           aGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAaccacc                                                                                                                     <  1:1394938/148‑1 (MQ=255)
                             aCATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCgg                                                                                                                   >  1:285534/1‑148 (MQ=255)
                             aCATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCgg                                                                                                                   >  1:281078/1‑148 (MQ=255)
                              catcatCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGt                                                                                                                  >  1:762231/1‑148 (MQ=255)
                                 catcaCGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGtata                                                                                                               >  1:965282/1‑148 (MQ=255)
                                 catcaCGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGtata                                                                                                               <  2:2661480/148‑1 (MQ=255)
                                   tcaCGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAgg                                                                                                             <  2:1934993/148‑1 (MQ=255)
                                   tcaCGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAgg                                                                                                             <  2:1935002/148‑1 (MQ=255)
                                   tcaCGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAgg                                                                                                             <  1:1969731/148‑1 (MQ=255)
                                    caCGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATGAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGga                                                                                                            >  2:3128177/1‑148 (MQ=255)
                                    caCGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGga                                                                                                            >  1:2197093/1‑148 (MQ=255)
                                     aCGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGgag                                                                                                           >  2:2720492/1‑148 (MQ=255)
                                       ggACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGaa                                                                                                         <  1:951253/148‑1 (MQ=255)
                                         aCAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAAc                                                                                                       <  1:2683810/148‑1 (MQ=255)
                                          cAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACa                                                                                                      >  1:1049178/1‑148 (MQ=255)
                                          cAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACa                                                                                                      <  2:2549403/148‑1 (MQ=255)
                                               tCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCa                                                                                                 >  1:1944281/1‑148 (MQ=255)
                                                 aGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGc                                                                                               >  1:1287235/1‑148 (MQ=255)
                                                 aGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGc                                                                                               <  2:3030584/148‑1 (MQ=255)
                                                   ccTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGcac                                                                                             >  1:648496/1‑148 (MQ=255)
                                                      tGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACg                                                                                          >  2:2046046/1‑148 (MQ=255)
                                                        cccAGATCCACTTGCGTCAGCACAACGGTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACACCCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGtt                                                                                        <  2:1896186/146‑1 (MQ=255)
                                                          cAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTa                                                                                      <  1:2018522/148‑1 (MQ=255)
                                                          cAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTa                                                                                      >  1:1842775/1‑148 (MQ=255)
                                                           aGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTAt                                                                                     <  1:130721/148‑1 (MQ=255)
                                                           aGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTAt                                                                                     >  2:1596562/1‑148 (MQ=255)
                                                            gATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTAtt                                                                                    <  1:3002825/148‑1 (MQ=255)
                                                            gATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAAGATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACNTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTAtt                                                                                    >  1:1678066/1‑148 (MQ=255)
                                                             aTCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTAttt                                                                                   <  2:1824332/148‑1 (MQ=255)
                                                              tCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTAtttt                                                                                  >  2:2706896/1‑148 (MQ=255)
                                                                cAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCa                                                                                >  2:973544/1‑148 (MQ=255)
                                                                   tttGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACg                                                                             >  2:56106/1‑148 (MQ=255)
                                                                    ttGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGt                                                                            <  1:244371/148‑1 (MQ=255)
                                                                           gCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAgcggcg                                                                     >  1:1597623/1‑148 (MQ=255)
                                                                            cTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGcggcgg                                                                    >  1:40007/1‑148 (MQ=255)
                                                                                   ttCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCCGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCa                                                             <  1:3077141/148‑1 (MQ=255)
                                                                                   ttCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCa                                                             <  1:2330539/148‑1 (MQ=255)
                                                                                    tCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAg                                                            <  2:1377097/148‑1 (MQ=255)
                                                                                        ccGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAggcgg                                                        <  2:600257/148‑1 (MQ=255)
                                                                                         cGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAggcagc                                                       <  1:2081763/148‑1 (MQ=255)
                                                                                            gCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTAt                                                    <  1:2035755/148‑1 (MQ=255)
                                                                                             cAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATc                                                   <  2:2535855/148‑1 (MQ=255)
                                                                                             cAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATc                                                   >  1:2068848/1‑148 (MQ=255)
                                                                                             cAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATc                                                   >  2:2576215/1‑148 (MQ=255)
                                                                                              aaTATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCa                                                  >  2:2441518/1‑148 (MQ=255)
                                                                                                tatCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTCTCCCg                                                >  2:804238/1‑148 (MQ=255)
                                                                                                 atCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGc                                               >  2:433239/1‑148 (MQ=255)
                                                                                                   ccAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTg                                             <  1:1901359/148‑1 (MQ=255)
                                                                                                       cgcTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTcc                                         >  2:1375122/1‑148 (MQ=255)
                                                                                                       cgcTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGGATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTcc                                         <  2:2466075/148‑1 (MQ=255)
                                                                                                        gcTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCa                                        <  1:3120608/148‑1 (MQ=255)
                                                                                                         cTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAg                                       <  1:1461223/148‑1 (MQ=255)
                                                                                                         cTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAg                                       <  1:1463318/148‑1 (MQ=255)
                                                                                                           tCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGgt                                     >  1:449845/1‑148 (MQ=255)
                                                                                                           tCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGgt                                     >  1:1981959/1‑148 (MQ=255)
                                                                                                              tGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGgtagt                                  <  2:1971529/148‑1 (MQ=255)
                                                                                                              tGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGgtagt                                  <  2:2383405/148‑1 (MQ=255)
                                                                                                                  aGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGcc                              <  1:3057631/148‑1 (MQ=255)
                                                                                                                  aGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGcc                              <  2:842951/148‑1 (MQ=255)
                                                                                                                     tATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACTGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCgg                           >  2:1850265/1‑147 (MQ=255)
                                                                                                                     tATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGt                           <  1:2554358/148‑1 (MQ=255)
                                                                                                                     tATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGt                           >  1:2589070/1‑148 (MQ=255)
                                                                                                                      aTCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTCGTAGCCCGtt                          >  2:2817107/1‑148 (MQ=255)
                                                                                                                      aTCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGtt                          <  2:3241486/148‑1 (MQ=255)
                                                                                                                          aCGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTtgt                      >  1:2539952/1‑148 (MQ=255)
                                                                                                                           cGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTtgtg                     >  2:2882058/1‑148 (MQ=255)
                                                                                                                           cGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTtgtg                     >  1:2771961/1‑148 (MQ=255)
                                                                                                                           cGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTGTtgtg                     >  2:2881921/1‑148 (MQ=255)
                                                                                                                             cTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTc                   <  2:2546711/148‑1 (MQ=255)
                                                                                                                              tGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTcc                  >  1:432713/1‑148 (MQ=255)
                                                                                                                               gACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCg                 <  1:2720492/148‑1 (MQ=255)
                                                                                                                               gACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCg                 >  2:979200/1‑148 (MQ=255)
                                                                                                                                 cAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTg               <  2:1821140/148‑1 (MQ=255)
                                                                                                                                      tGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGaa          <  2:527052/148‑1 (MQ=255)
                                                                                                                                      tGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGaa          >  1:3104258/1‑148 (MQ=255)
                                                                                                                                       ggCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAat         <  2:3250008/148‑1 (MQ=255)
                                                                                                                                        gCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAata        >  2:2248085/1‑148 (MQ=255)
                                                                                                                                        gCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTACAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAata        <  2:1536844/148‑1 (MQ=255)
                                                                                                                                         cTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAatat       >  2:1694136/1‑148 (MQ=255)
                                                                                                                                          tGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATc      >  2:1339785/1‑148 (MQ=255)
                                                                                                                                           gCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATcc     <  2:762231/148‑1 (MQ=255)
                                                                                                                                           gCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATcc     <  1:1214687/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GATCGCGATAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACC  >  CP009273/1024558‑1024847

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: