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,7030GT95.9% 741.0 / 17.5 244G241G (GGC→GGArlmI23S rRNA m(5)C1962 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base G (7/3);  new base T (119/114);  total (127/117)
Fisher's exact test for biased strand distribution p-value = 3.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 244 total aligned reads displayed.

GATCGCGATAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACC  >  CP009273/1024558‑1024847
                                                                                                                                                 |                                                                                                                                                
 atcgcgatAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTccc                                                                                                                                               >  2:65181/1‑148 (MQ=255)
  tcgcgatAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCa                                                                                                                                              >  1:98063/1‑148 (MQ=255)
       atAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCACCCTCCCATCAGt                                                                                                                                         >  1:2720156/1‑148 (MQ=255)
       atAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGt                                                                                                                                         >  2:1467031/1‑148 (MQ=255)
       atAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGt                                                                                                                                         >  1:50927/1‑148 (MQ=255)
        tAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTg                                                                                                                                        >  2:2607344/1‑148 (MQ=255)
           gTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCg                                                                                                                                     >  2:2422820/1‑148 (MQ=255)
             aCGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTCTCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGAt                                                                                                                                   >  2:3329228/1‑148 (MQ=255)
              cGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATa                                                                                                                                  <  2:2276136/148‑1 (MQ=255)
                 agcaATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACgg                                                                                                                               <  1:3010872/148‑1 (MQ=255)
                   caATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCg                                                                                                                             >  1:3718913/1‑148 (MQ=255)
                     aTTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGaa                                                                                                                           >  2:2095099/1‑148 (MQ=255)
                     aTTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGaa                                                                                                                           >  2:2158249/1‑148 (MQ=255)
                     aTTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGaa                                                                                                                           >  2:2158784/1‑148 (MQ=255)
                     aTTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGaa                                                                                                                           <  1:134886/148‑1 (MQ=255)
                     aTTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGaa                                                                                                                           <  1:2089401/148‑1 (MQ=255)
                      tttAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGaaa                                                                                                                          <  1:2534997/148‑1 (MQ=255)
                       ttAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCCTCAGTGCCGATACGGCGAAac                                                                                                                         >  1:1591319/1‑148 (MQ=255)
                       ttAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACTGCGAAac                                                                                                                         >  2:160829/1‑148 (MQ=255)
                       ttAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAac                                                                                                                         >  2:3501817/1‑148 (MQ=255)
                        tAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTTTTCAGCTCAAAGATCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATAAACGATGACAACCTGGCTCCAACCTCCCATCAGAGCCGATACTGCGAAacc                                                                                                                        >  1:3198493/1‑148 (MQ=255)
                            gACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCg                                                                                                                    <  2:880127/148‑1 (MQ=255)
                            gACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCg                                                                                                                    >  2:452718/1‑148 (MQ=255)
                              catcatCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGt                                                                                                                  >  2:2478744/1‑148 (MQ=255)
                               atcatcACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGta                                                                                                                 <  1:1984749/148‑1 (MQ=255)
                                tcatcaCGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGtat                                                                                                                >  1:1225200/1‑148 (MQ=255)
                                tcatcaCGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGtat                                                                                                                >  1:3129463/1‑148 (MQ=255)
                                 catcaCGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGtata                                                                                                               <  2:779229/148‑1 (MQ=255)
                                   tcaCGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAgg                                                                                                             >  1:545495/1‑148 (MQ=255)
                                   tcaCGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAgg                                                                                                             <  1:3048353/148‑1 (MQ=255)
                                    caCGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGga                                                                                                            >  1:2744884/1‑148 (MQ=255)
                                        gACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGaaa                                                                                                        <  1:2496690/148‑1 (MQ=255)
                                          cAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACa                                                                                                      <  1:1660178/148‑1 (MQ=255)
                                          cAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACa                                                                                                      >  2:398553/1‑148 (MQ=255)
                                             aCTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGtt                                                                                                   <  1:1298480/148‑1 (MQ=255)
                                             aCTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGtt                                                                                                   >  1:165552/1‑148 (MQ=255)
                                             aCTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGtt                                                                                                   >  2:1664684/1‑148 (MQ=255)
                                                cAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAg                                                                                                >  1:3210207/1‑148 (MQ=255)
                                                cAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAg                                                                                                >  1:3210491/1‑148 (MQ=255)
                                                           aGATCCAGTTTGTTCAGCTCAACGTTCGGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTAt                                                                                     <  2:3686195/148‑1 (MQ=255)
                                                              tCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTAtttt                                                                                  <  1:3065371/148‑1 (MQ=255)
                                                              tCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTAtttt                                                                                  >  1:2011245/1‑148 (MQ=255)
                                                              tCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTAtttt                                                                                  >  2:67229/1‑148 (MQ=255)
                                                                 aGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCaa                                                                               >  2:340574/1‑148 (MQ=255)
                                                                 aGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCaa                                                                               >  2:2032334/1‑148 (MQ=255)
                                                                   tttGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACg                                                                             >  2:1169438/1‑148 (MQ=255)
                                                                   tttGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACg                                                                             >  2:126422/1‑148 (MQ=255)
                                                                   tttGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACg                                                                             <  1:1337289/148‑1 (MQ=255)
                                                                   tttGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGATTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTTCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACg                                                                             >  1:2471819/1‑148 (MQ=255)
                                                                    ttGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGt                                                                            >  2:2051247/1‑148 (MQ=255)
                                                                    ttGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGt                                                                            <  2:17463/148‑1 (MQ=255)
                                                                        tCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAgcg                                                                        <  1:65181/148‑1 (MQ=255)
                                                                            cTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGcggcgg                                                                    >  2:2516927/1‑148 (MQ=255)
                                                                            cTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGcggcgg                                                                    >  1:2963866/1‑148 (MQ=255)
                                                                             tCAACGTTCTGCCGTGCAATATCGAGCGCTTCCTGGGATGTATCAACGCTGACATACTGGCGGGAACCTCCCATCAGTGCCGATACTGCGACACCACCGTTATTGGATTAACAGTTCTGCACACTTTTTTTGTCAACGTAGCGGCggg                                                                   <  2:3198493/148‑1 (MQ=255)
                                                                              cAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGt                                                                  <  1:292035/148‑1 (MQ=255)
                                                                              cAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGt                                                                  >  2:2204079/1‑148 (MQ=255)
                                                                               aaCGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGtc                                                                 >  2:2678281/1‑147 (MQ=255)
                                                                                aCGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAg                                                                <  1:1716479/148‑1 (MQ=255)
                                                                                 cGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGc                                                               <  2:1396559/148‑1 (MQ=255)
                                                                                   ttCTGCCGTGCAATAGCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCa                                                             <  2:2748085/148‑1 (MQ=255)
                                                                                    tCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAg                                                            >  2:3439420/1‑148 (MQ=255)
                                                                                        ccGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAggcgg                                                        <  1:2406408/148‑1 (MQ=255)
                                                                                         cGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAggcggc                                                       >  1:1488329/1‑148 (MQ=255)
                                                                                            gCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTAt                                                    >  2:208802/1‑148 (MQ=255)
                                                                                            gCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTAt                                                    >  1:1241197/1‑148 (MQ=255)
                                                                                            gCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTAt                                                    >  2:231534/1‑148 (MQ=255)
                                                                                             cAATATCCAGCGCTTCCTGGGATGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATc                                                   >  2:2999663/1‑148 (MQ=255)
                                                                                               atatCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCAc                                                 >  2:2813428/1‑148 (MQ=255)
                                                                                               atatCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCAc                                                 >  1:1115489/1‑148 (MQ=255)
                                                                                               atatCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCAc                                                 >  2:381587/1‑148 (MQ=255)
                                                                                                tatCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACg                                                >  1:757169/1‑148 (MQ=255)
                                                                                                 atCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGc                                               <  2:1225200/148‑1 (MQ=255)
                                                                                                 atCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAATGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGc                                               <  2:1574224/148‑1 (MQ=255)
                                                                                                      gcgcTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTc                                          <  2:1429441/148‑1 (MQ=255)
                                                                                                       cgcTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTcc                                         <  2:2744884/148‑1 (MQ=255)
                                                                                                       cgcTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTcc                                         >  1:1824710/1‑148 (MQ=255)
                                                                                                        gcTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCa                                        <  2:1488546/148‑1 (MQ=255)
                                                                                                         cTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAg                                       <  1:2051247/148‑1 (MQ=255)
                                                                                                          ttCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAgg                                      <  2:314966/148‑1 (MQ=255)
                                                                                                             cTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGgtag                                   <  2:1315197/148‑1 (MQ=255)
                                                                                                              tGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGgtagt                                  <  2:3231582/148‑1 (MQ=255)
                                                                                                                ggAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGtagtag                                <  2:2708652/148‑1 (MQ=255)
                                                                                                                    gTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCg                            >  1:428202/1‑148 (MQ=255)
                                                                                                                    gTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCg                            <  2:2888005/148‑1 (MQ=255)
                                                                                                                     tATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGt                           >  2:2980918/1‑148 (MQ=255)
                                                                                                                      aTCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGtt                          >  1:1002664/1‑148 (MQ=255)
                                                                                                                      aTCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGtt                          <  1:1664684/148‑1 (MQ=255)
                                                                                                                        cAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGtttt                        >  1:3563870/1‑148 (MQ=255)
                                                                                                                         aaCGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTtg                       >  1:1015064/1‑148 (MQ=255)
                                                                                                                           cGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTtgtg                     >  2:326212/1‑148 (MQ=255)
                                                                                                                              tGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGAGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTcc                  <  1:712381/148‑1 (MQ=255)
                                                                                                                               gACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCg                 >  1:1374190/1‑148 (MQ=255)
                                                                                                                               gACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCg                 >  1:827685/1‑148 (MQ=255)
                                                                                                                                aCAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTGTTGTGTCCGt                >  1:2361634/1‑148 (MQ=255)
                                                                                                                                 cAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTg               <  2:485421/148‑1 (MQ=255)
                                                                                                                                       ggCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAat         <  1:1491378/148‑1 (MQ=255)
                                                                                                                                          tGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATc      <  1:3501817/148‑1 (MQ=255)
                                                                                                                                              aCCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATccacc  >  2:3602971/1‑148 (MQ=255)
                                                                                                                                              aCCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATccacc  >  2:2949198/1‑148 (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: