Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,024,703 G→T 100% 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,7030GT97.7% 678.9 / 6.0 217G241G (GGC→GGArlmI23S rRNA m(5)C1962 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base T (100/112);  total (103/114)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 217 total aligned reads displayed.

GATCGCGATAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCA  >  CP009273/1024558‑1024848
                                                                                                                                                 |                                                                                                                                                 
 atcgcgatAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTccc                                                                                                                                                <  2:3541433/148‑1 (MQ=255)
  tcgcgatAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCa                                                                                                                                               >  2:280319/1‑148 (MQ=255)
  tcgcgatAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCACTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTcccc                                                                                                                                               >  2:2197344/1‑147 (MQ=255)
   cgcgatAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCAt                                                                                                                                              <  1:3158035/148‑1 (MQ=255)
     cgatAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCa                                                                                                                                            <  2:1880060/148‑1 (MQ=255)
     cgatAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCa                                                                                                                                            <  2:545667/148‑1 (MQ=255)
        tAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTg                                                                                                                                         >  2:424871/1‑148 (MQ=255)
         aaGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGc                                                                                                                                        >  2:202731/1‑148 (MQ=255)
          aGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGcc                                                                                                                                       >  2:159378/1‑148 (MQ=255)
             aCGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGAt                                                                                                                                    >  1:3445478/1‑148 (MQ=255)
             aCGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGAt                                                                                                                                    <  1:1714424/148‑1 (MQ=255)
                cagcaATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACg                                                                                                                                 >  1:487805/1‑148 (MQ=255)
                  gcaATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGc                                                                                                                               <  1:2607505/148‑1 (MQ=255)
                          aaGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAaccac                                                                                                                       >  1:841100/1‑148 (MQ=255)
                          aaGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAaccac                                                                                                                       <  1:69229/148‑1 (MQ=255)
                          aaGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAaccac                                                                                                                       <  1:71685/148‑1 (MQ=255)
                          aaGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAaccac                                                                                                                       >  1:841083/1‑148 (MQ=255)
                          aaGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAaccac                                                                                                                       <  1:1300220/148‑1 (MQ=255)
                              catcatCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGt                                                                                                                   <  2:1704844/148‑1 (MQ=255)
                              catcatCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGt                                                                                                                   <  2:358129/148‑1 (MQ=255)
                                tcatcaCGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGtat                                                                                                                 <  1:2001542/148‑1 (MQ=255)
                                    caCGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGga                                                                                                             <  2:884885/148‑1 (MQ=255)
                                         aCAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAAc                                                                                                        <  1:2285025/148‑1 (MQ=255)
                                          cAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACa                                                                                                       >  2:126924/1‑148 (MQ=255)
                                                cAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAg                                                                                                 <  2:2920628/148‑1 (MQ=255)
                                                     ttGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCacac                                                                                            <  1:2527594/148‑1 (MQ=255)
                                                     ttGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCacac                                                                                            >  1:1005621/1‑148 (MQ=255)
                                                      tGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACg                                                                                           >  2:292198/1‑148 (MQ=255)
                                                      tGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACg                                                                                           <  1:2933435/148‑1 (MQ=255)
                                                       gCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGt                                                                                          <  2:2157659/148‑1 (MQ=255)
                                                         tCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGttt                                                                                        <  2:1613784/148‑1 (MQ=255)
                                                         tCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGttt                                                                                        <  2:193276/148‑1 (MQ=255)
                                                         tCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGttt                                                                                        <  2:2679838/148‑1 (MQ=255)
                                                          cAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTa                                                                                       <  1:697579/148‑1 (MQ=255)
                                                          cAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTa                                                                                       <  1:116981/148‑1 (MQ=255)
                                                           aGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTAt                                                                                      <  1:2162463/148‑1 (MQ=255)
                                                            gATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTAtt                                                                                     <  2:894651/148‑1 (MQ=255)
                                                            gATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTAtt                                                                                     >  1:368022/1‑148 (MQ=255)
                                                             aTCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTAttt                                                                                    >  2:13206/1‑148 (MQ=255)
                                                               ccAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTc                                                                                  >  2:1114400/1‑148 (MQ=255)
                                                               ccAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTc                                                                                  >  2:3525619/1‑148 (MQ=255)
                                                               ccAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTc                                                                                  >  2:687346/1‑148 (MQ=255)
                                                                cAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCa                                                                                 <  1:611227/148‑1 (MQ=255)
                                                                 aGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCaa                                                                                <  2:2122074/148‑1 (MQ=255)
                                                                   tttGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACg                                                                              <  2:2153367/148‑1 (MQ=255)
                                                                   tttGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACg                                                                              >  1:1837048/1‑148 (MQ=255)
                                                                    ttGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGt                                                                             >  2:3333276/1‑148 (MQ=255)
                                                                    ttGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGt                                                                             <  1:278721/148‑1 (MQ=255)
                                                                    ttGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCGTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGt                                                                             <  2:1809094/148‑1 (MQ=255)
                                                                         cAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAgcgg                                                                        <  2:255055/148‑1 (MQ=255)
                                                                              cAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGt                                                                   >  1:3231629/1‑148 (MQ=255)
                                                                                aCGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATTAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAg                                                                 >  2:2462529/1‑148 (MQ=255)
                                                                                aCGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAg                                                                 >  2:1773324/1‑148 (MQ=255)
                                                                                 cGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGc                                                                >  2:1234484/1‑148 (MQ=255)
                                                                                 cGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGc                                                                >  1:1261229/1‑148 (MQ=255)
                                                                                   ttCTGCCGTGCAATCTCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCa                                                              <  2:153378/148‑1 (MQ=255)
                                                                                   ttCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCa                                                              <  2:800336/148‑1 (MQ=255)
                                                                                        ccGTGCAATATCCCGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAggcgg                                                         <  2:1403262/148‑1 (MQ=255)
                                                                                        ccGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAggcgg                                                         <  2:1508423/148‑1 (MQ=255)
                                                                                          gTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCt                                                       <  2:651762/148‑1 (MQ=255)
                                                                                           tGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTa                                                      <  1:3125/148‑1 (MQ=255)
                                                                                               atatCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCAc                                                  >  2:2743417/1‑148 (MQ=255)
                                                                                                tatCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACg                                                 <  1:3519145/148‑1 (MQ=255)
                                                                                                 atCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGc                                                >  1:1608790/1‑148 (MQ=255)
                                                                                                   ccAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTg                                              >  2:2338175/1‑148 (MQ=255)
                                                                                                    cAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTgg                                             <  1:29290/148‑1 (MQ=255)
                                                                                                      gcgcTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTc                                           <  2:710180/148‑1 (MQ=255)
                                                                                                       cgcTTCCTGGGAGGTATCAACGCTGACAACCTGGCTCCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTcc                                          >  2:2471719/1‑148 (MQ=255)
                                                                                                           tCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGgt                                      <  1:1777674/148‑1 (MQ=255)
                                                                                                              tGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGgtagt                                   >  1:1195157/1‑148 (MQ=255)
                                                                                                              tGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGgtagt                                   <  1:980494/148‑1 (MQ=255)
                                                                                                                 gAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGc                                >  2:631330/1‑148 (MQ=255)
                                                                                                                 gAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGc                                >  1:1941566/1‑148 (MQ=255)
                                                                                                                 gAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGc                                >  1:357677/1‑148 (MQ=255)
                                                                                                                   ggTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGccc                              >  2:1867784/1‑148 (MQ=255)
                                                                                                                    gTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAGCCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCg                             >  1:3334849/1‑148 (MQ=255)
                                                                                                                    gTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCg                             <  2:582412/148‑1 (MQ=255)
                                                                                                                    gTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCg                             >  2:3489366/1‑148 (MQ=255)
                                                                                                                    gTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCg                             <  1:629698/148‑1 (MQ=255)
                                                                                                                     tATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGt                            >  1:2032493/1‑148 (MQ=255)
                                                                                                                      aTCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGtt                           <  1:277565/148‑1 (MQ=255)
                                                                                                                      aTCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCGATCACGCTGGTCCAGGTCGTAGCCCGtt                           >  1:932993/1‑148 (MQ=255)
                                                                                                                           cGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTtgtg                      >  2:3151972/1‑148 (MQ=255)
                                                                                                                            gCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTgtgt                     <  2:841083/148‑1 (MQ=255)
                                                                                                                            gCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTgtgt                     <  2:466158/148‑1 (MQ=255)
                                                                                                                            gCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTgtgt                     <  2:841100/148‑1 (MQ=255)
                                                                                                                            gCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTgtgt                     >  2:2867918/1‑148 (MQ=255)
                                                                                                                               gACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCg                  <  1:1642690/148‑1 (MQ=255)
                                                                                                                                aCAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGt                 <  2:1124445/148‑1 (MQ=255)
                                                                                                                                     cTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGa            <  1:1675149/148‑1 (MQ=255)
                                                                                                                                     cTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGa            <  1:2144575/148‑1 (MQ=255)
                                                                                                                                       ggCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAat          <  1:473236/148‑1 (MQ=255)
                                                                                                                                       ggCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAat          <  1:2437154/148‑1 (MQ=255)
                                                                                                                                        gCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAata         <  2:102763/148‑1 (MQ=255)
                                                                                                                                         cTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAatat        <  2:3155615/148‑1 (MQ=255)
                                                                                                                                         cTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAatat        <  1:3078280/148‑1 (MQ=255)
                                                                                                                                            cAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATcca     >  1:2897028/1‑148 (MQ=255)
                                                                                                                                              aCCTCCCATCAGTGCCGATACGGTGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATccacc   >  1:3443099/1‑148 (MQ=255)
                                                                                                                                              aCCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATccacc   <  1:1306752/148‑1 (MQ=255)
                                                                                                                                              aCCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATccacc   <  1:753647/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GATCGCGATAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCA  >  CP009273/1024558‑1024848

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: