Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,024,703 G→T 40.6% 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,7030GT40.6% 396.1 / inf 349G241G (GGC→GGArlmI23S rRNA m(5)C1962 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base G (111/96);  new base T (60/82);  total (171/178)
Fisher's exact test for biased strand distribution p-value = 3.90e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 351 total aligned reads displayed.

CGCGATCGCGATAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAG  >  CP009273/1024555‑1024852
                                                                                                                                                    |                                                                                                                                                     
   gatcgcgatAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCa                                                                                                                                                   >  2:691574/1‑150 (MQ=255)
   gatcgcgatAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCa                                                                                                                                                   >  2:1350650/1‑150 (MQ=255)
    atcgcgatAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCAt                                                                                                                                                  <  1:677270/150‑1 (MQ=255)
      cgcgatAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCa                                                                                                                                                <  1:604506/150‑1 (MQ=255)
       gcgatAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAg                                                                                                                                               >  1:2871106/1‑150 (MQ=255)
         gatAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTg                                                                                                                                             <  1:1538717/150‑1 (MQ=255)
         gatAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTg                                                                                                                                             <  2:631860/150‑1 (MQ=255)
          atAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGc                                                                                                                                            >  1:29134/1‑150 (MQ=255)
          atAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGc                                                                                                                                            <  1:4637809/150‑1 (MQ=255)
             aGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGa                                                                                                                                         <  2:869830/150‑1 (MQ=255)
              gTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGAt                                                                                                                                        >  2:715294/1‑150 (MQ=255)
              gTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGAt                                                                                                                                        <  2:4443799/150‑1 (MQ=255)
              gTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGAt                                                                                                                                        >  2:947806/1‑150 (MQ=255)
               tACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATa                                                                                                                                       >  2:427532/1‑150 (MQ=255)
                aCGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATAc                                                                                                                                      <  1:3762712/150‑1 (MQ=255)
                    agcaATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCg                                                                                                                                  >  2:2968065/1‑150 (MQ=255)
                      caATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGaa                                                                                                                                >  1:2541855/1‑150 (MQ=255)
                         tttAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAacc                                                                                                                             <  1:3197590/150‑1 (MQ=255)
                          ttAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAacca                                                                                                                            >  1:705940/1‑150 (MQ=255)
                           tAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAaccac                                                                                                                           >  2:4236635/1‑150 (MQ=255)
                              aGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCgg                                                                                                                        <  2:3576879/150‑1 (MQ=255)
                               gACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGt                                                                                                                       >  2:707213/1‑150 (MQ=255)
                                 catcatCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGtat                                                                                                                     >  2:1976390/1‑150 (MQ=255)
                                  atcatcACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGtata                                                                                                                    >  2:1210742/1‑150 (MQ=255)
                                  atcatcACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGtata                                                                                                                    <  1:4086242/150‑1 (MQ=255)
                                    catcaCGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAgg                                                                                                                  >  1:3890627/1‑150 (MQ=255)
                                     atcaCGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGga                                                                                                                 >  2:5073854/1‑150 (MQ=255)
                                     atcaCGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGga                                                                                                                 >  1:1944982/1‑150 (MQ=255)
                                      tcaCGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGgag                                                                                                                >  1:1440349/1‑150 (MQ=255)
                                       caCGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGgaga                                                                                                               >  1:1517023/1‑150 (MQ=255)
                                        aCGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGaa                                                                                                              >  1:3722173/1‑150 (MQ=255)
                                         cGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACACCCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGaaa                                                                                                             >  2:1134771/1‑150 (MQ=255)
                                         cGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGaaa                                                                                                             >  2:3706102/1‑150 (MQ=255)
                                          ggACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAGCGCTGACAACCTGGCTGCAACCGCCCCTCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAAc                                                                                                            >  2:4518005/1‑150 (MQ=255)
                                          ggACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAAc                                                                                                            >  1:2169162/1‑150 (MQ=255)
                                           gACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACa                                                                                                           <  2:1530953/150‑1 (MQ=255)
                                              aaaCTCAGCCTTGCTCAGATCCAGTTTGTTCATCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGCAACCACCGGTATCGGAGAAACAGtt                                                                                                        >  2:4786564/1‑150 (MQ=255)
                                              aaaCTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGtt                                                                                                        >  2:2657864/1‑150 (MQ=255)
                                              aaaCTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGtt                                                                                                        <  1:1779257/150‑1 (MQ=255)
                                                  tCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGc                                                                                                    >  1:3156534/1‑150 (MQ=255)
                                                  tCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGc                                                                                                    >  1:2666387/1‑150 (MQ=255)
                                                  tCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGc                                                                                                    >  1:5007962/1‑150 (MQ=255)
                                                      ccTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCacac                                                                                                <  1:3245357/150‑1 (MQ=255)
                                                       cTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACg                                                                                               >  2:2773646/1‑150 (MQ=255)
                                                        ttGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGt                                                                                              <  2:653096/150‑1 (MQ=255)
                                                        ttGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGt                                                                                              <  1:1804095/150‑1 (MQ=255)
                                                          gCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCACCCTCCCATCAGTGCCGATACGGCGAAACCACCGGTCTAGGAGAAACAGATCAGCACACGttt                                                                                            >  2:3495406/1‑150 (MQ=255)
                                                          gCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGttt                                                                                            <  1:4444704/150‑1 (MQ=255)
                                                            tCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTAt                                                                                          <  2:1601339/150‑1 (MQ=255)
                                                              aGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTAttt                                                                                        >  2:407859/1‑150 (MQ=255)
                                                               gATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTAtttt                                                                                       <  1:1609958/150‑1 (MQ=255)
                                                               gATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTAtttt                                                                                       >  2:2862415/1‑150 (MQ=255)
                                                                aTCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTc                                                                                      <  1:852680/150‑1 (MQ=255)
                                                                 tCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCa                                                                                     <  1:2334810/150‑1 (MQ=255)
                                                                 tCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCa                                                                                     >  2:385318/1‑150 (MQ=255)
                                                                  ccAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTACAGCACACGTTTATTTTCaa                                                                                    >  2:2447840/1‑150 (MQ=255)
                                                                   cAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAAc                                                                                   <  2:4588490/150‑1 (MQ=255)
                                                                    aGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACg                                                                                  <  1:4552445/150‑1 (MQ=255)
                                                                      ttagTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAg                                                                               <  2:2131321/147‑1 (MQ=255)
                                                                         gTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAgcg                                                                             <  1:3523712/150‑1 (MQ=255)
                                                                         gTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAgcg                                                                             >  1:2996595/1‑150 (MQ=255)
                                                                             aGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGcggcgg                                                                         <  2:524221/150‑1 (MQ=255)
                                                                                   aCGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGcc                                                                   <  1:1196182/150‑1 (MQ=255)
                                                                                     gTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAg                                                                 >  1:3944576/1‑150 (MQ=255)
                                                                                       tCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAggc                                                               >  2:3214989/1‑150 (MQ=255)
                                                                                        cTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAggcg                                                              >  1:1662297/1‑150 (MQ=255)
                                                                                            cGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTa                                                          <  2:2528218/150‑1 (MQ=255)
                                                                                            cGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTa                                                          >  1:4295297/1‑150 (MQ=255)
                                                                                            cGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTa                                                          >  1:4046216/1‑150 (MQ=255)
                                                                                             gTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTAt                                                         <  1:3095622/150‑1 (MQ=255)
                                                                                              tGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATc                                                        <  1:1858023/150‑1 (MQ=255)
                                                                                              tGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATc                                                        <  1:411722/150‑1 (MQ=255)
                                                                                               gCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCa                                                       >  1:2259298/1‑150 (MQ=255)
                                                                                               gCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCa                                                       <  2:835093/150‑1 (MQ=255)
                                                                                                cAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCAc                                                      <  2:3819055/150‑1 (MQ=255)
                                                                                                 aaTATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACg                                                     <  1:3489111/150‑1 (MQ=255)
                                                                                                    atCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTg                                                  <  1:4880242/150‑1 (MQ=255)
                                                                                                    atCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTg                                                  >  1:3330716/1‑150 (MQ=255)
                                                                                                       ccgcgcGTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTc                                               <  2:4883291/148‑1 (MQ=255)
                                                                                                        aGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTcc                                              >  1:1006782/1‑150 (MQ=255)
                                                                                                            cTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTCGCCAGGCGGCTATCACGCTGGTCCAGgt                                          >  2:621446/1‑150 (MQ=255)
                                                                                                            cTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGgt                                          <  1:4771981/150‑1 (MQ=255)
                                                                                                                cTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGgtagta                                      <  2:147062/150‑1 (MQ=255)
                                                                                                                cTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGgtagta                                      <  1:503570/150‑1 (MQ=255)
                                                                                                                 tGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGtagtag                                     >  2:1886603/1‑150 (MQ=255)
                                                                                                                    gAGGTATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGccc                                  <  2:1444629/150‑1 (MQ=255)
                                                                                                                      ggTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGt                                >  2:417791/1‑150 (MQ=255)
                                                                                                                        tATCAACGCTGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTGTATTTTCAACGTATCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGtgt                              >  2:3536404/1‑148 (MQ=255)
                                                                                                                        tATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGttt                              <  2:1025434/150‑1 (MQ=255)
                                                                                                                        tATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGttt                              >  2:1302563/1‑150 (MQ=255)
                                                                                                                          tCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTtg                            >  2:4129716/1‑150 (MQ=255)
                                                                                                                           cAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTtgt                           <  2:4844822/150‑1 (MQ=255)
                                                                                                                                cTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCg                      >  1:119559/1‑150 (MQ=255)
                                                                                                                                 tGACAACCTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCCGGCGGCTCTCACGCTGGTCCAGGTAGTAGCCCGTTGTGTGTCCGt                     >  2:3218722/1‑150 (MQ=255)
                                                                                                                                     aaCCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTg                 <  1:4609822/150‑1 (MQ=255)
                                                                                                                                     aaCCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGGTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTCGTAGCCCGTGTTGTGTCCGTGCTg                 >  2:4292708/1‑150 (MQ=255)
                                                                                                                                       ccTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGaa               <  1:1392650/150‑1 (MQ=255)
                                                                                                                                        cTGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAat              <  1:249137/150‑1 (MQ=255)
                                                                                                                                         tGGCTGCAACCTCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAata             <  1:529413/150‑1 (MQ=255)
                                                                                                                                            cTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATcc          >  1:1000284/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGCGATCGCGATAAGTACGCAGCAATTTAAAGACATCATCACGGACAAACTCAGCCTTGCTCAGATCCAGTTTGTTCAGCTCAACGTTCTGCCGTGCAATATCCAGCGCTTCCTGGGAGGTATCAACGCTGACAACCTGGCTGCAACCGCCCATCAGTGCCGATACGGCGAAACCACCGGTATAGGAGAAACAGTTCAGCACACGTTTATTTTCAACGTAGCGGCGGGTAGCCAGGCGGCTATCACGCTGGTCCAGGTAGTAGCCCGTTTTGTGTCCGTGCTGAATATCCACCAGCAG  >  CP009273/1024555‑1024852

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: