Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,003,597 G→T 92.3% D100Y (GAC→TAC)  rlmL → fused 23S rRNA m(2)G2445 and m(7)G2069 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,003,5970GT92.3% 427.9 / 28.8 156D100Y (GAC→TAC) rlmLfused 23S rRNA m(2)G2445 and m(7)G2069 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base G (4/8);  new base T (79/65);  total (83/73)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 156 total aligned reads displayed.

CCGCCTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGAAACCG  >  CP009273/1003452‑1003741
                                                                                                                                                 |                                                                                                                                                
  gccTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTacac                                                                                                                                              >  1:502687/1‑148 (MQ=255)
  gccTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTacac                                                                                                                                              <  1:585530/148‑1 (MQ=255)
  gccTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTacac                                                                                                                                              <  2:222130/148‑1 (MQ=255)
      ggCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATc                                                                                                                                          >  1:2084135/1‑148 (MQ=255)
            gcgTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAAc                                                                                                                                    <  2:319361/148‑1 (MQ=255)
              gTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAAcag                                                                                                                                  <  2:932950/148‑1 (MQ=255)
                attatGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAAcagtc                                                                                                                                >  1:3043487/1‑148 (MQ=255)
                   atGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAAcagtcagt                                                                                                                             >  2:1531600/1‑148 (MQ=255)
                    tGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTa                                                                                                                            >  2:978538/1‑148 (MQ=255)
                    tGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTa                                                                                                                            <  2:2978198/148‑1 (MQ=255)
                            cTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGa                                                                                                                    <  1:918336/148‑1 (MQ=255)
                                gCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtgaa                                                                                                                >  2:2625723/1‑148 (MQ=255)
                                  gAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtgaaag                                                                                                              <  1:3266801/148‑1 (MQ=255)
                                   aGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtgaaagt                                                                                                             >  1:3080799/1‑148 (MQ=255)
                                   aGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtgaaagt                                                                                                             >  1:2053217/1‑148 (MQ=255)
                                   aGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtgaaagt                                                                                                             <  2:2896697/148‑1 (MQ=255)
                                     tgtAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtgaaagtga                                                                                                           >  1:434423/1‑148 (MQ=255)
                                      gtAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtgaaagtgaa                                                                                                          <  1:2422097/148‑1 (MQ=255)
                                        aaGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtgaaagtgaaag                                                                                                        >  1:116864/1‑148 (MQ=255)
                                          ggTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGAc                                                                                                      >  1:1758869/1‑148 (MQ=255)
                                          ggTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGAc                                                                                                      >  1:1256945/1‑148 (MQ=255)
                                           gTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGAcg                                                                                                     >  2:2798878/1‑148 (MQ=255)
                                           gTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGAcg                                                                                                     >  2:302768/1‑148 (MQ=255)
                                            tttACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGAcgc                                                                                                    >  1:2906782/1‑148 (MQ=255)
                                              tACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGa                                                                                                  <  1:2936207/148‑1 (MQ=255)
                                               aCAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGAt                                                                                                 >  2:2703410/1‑148 (MQ=255)
                                               aCAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCGTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGAt                                                                                                 <  2:2368892/148‑1 (MQ=255)
                                                cAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGAtc                                                                                                >  1:1142900/1‑148 (MQ=255)
                                                 aGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGAtcg                                                                                               <  2:2125655/148‑1 (MQ=255)
                                                      tttAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGAt                                                                                          >  2:1887807/1‑148 (MQ=255)
                                                      tttAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGAt                                                                                          >  2:1888340/1‑148 (MQ=255)
                                                       ttAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATg                                                                                         >  1:1880739/1‑148 (MQ=255)
                                                        tAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGc                                                                                        >  1:1827575/1‑148 (MQ=255)
                                                         aGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCt                                                                                       <  2:34538/148‑1 (MQ=255)
                                                         aGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCt                                                                                       <  2:1763101/148‑1 (MQ=255)
                                                           aCCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCttt                                                                                     >  2:843565/1‑148 (MQ=255)
                                                            ccTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTc                                                                                    >  2:2449578/1‑148 (MQ=255)
                                                            ccTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTc                                                                                    >  1:1464880/1‑148 (MQ=255)
                                                            ccTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTc                                                                                    >  1:394807/1‑148 (MQ=255)
                                                              tctATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCAc                                                                                  >  1:965117/1‑148 (MQ=255)
                                                                  tctcGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCgg                                                                              <  1:1347466/148‑1 (MQ=255)
                                                                  tctcGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCgg                                                                              <  1:2181835/148‑1 (MQ=255)
                                                                      ggTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGaaaa                                                                          >  2:3084133/1‑148 (MQ=255)
                                                                      ggTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGaaaa                                                                          <  2:992407/148‑1 (MQ=255)
                                                                      ggTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGaaaa                                                                          <  1:771975/148‑1 (MQ=255)
                                                                      ggTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGaaaa                                                                          <  1:1018497/148‑1 (MQ=255)
                                                                       gtgtTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGaaaaa                                                                         >  1:155571/1‑148 (MQ=255)
                                                                         gtTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATc                                                                       <  2:2660541/148‑1 (MQ=255)
                                                                          ttCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCt                                                                      >  2:471909/1‑148 (MQ=255)
                                                                          ttCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCt                                                                      <  1:3195678/148‑1 (MQ=255)
                                                                           tCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTg                                                                     <  1:751578/148‑1 (MQ=255)
                                                                             aGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGcc                                                                   <  2:1801927/148‑1 (MQ=255)
                                                                              ggCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCcg                                                                  >  2:1446373/1‑148 (MQ=255)
                                                                                   tCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTcc                                                             >  1:3153010/1‑148 (MQ=255)
                                                                                   tCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTcc                                                             >  2:2057501/1‑148 (MQ=255)
                                                                                     aaCTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCaa                                                           <  1:1873316/148‑1 (MQ=255)
                                                                                     aaCTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCaa                                                           <  2:575603/148‑1 (MQ=255)
                                                                                        tGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATg                                                        <  1:1385625/148‑1 (MQ=255)
                                                                                        tGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATg                                                        <  1:3032603/148‑1 (MQ=255)
                                                                                         ggACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCTTTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGt                                                       <  2:2601345/148‑1 (MQ=255)
                                                                                          gACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGtt                                                      >  1:2276430/1‑148 (MQ=255)
                                                                                          gACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGtt                                                      >  1:23136/1‑148 (MQ=255)
                                                                                           aCAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGACAAATCTGCCGCGTCCAAATGTTg                                                     >  2:3032083/1‑148 (MQ=255)
                                                                                           aCAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTg                                                     <  1:2703410/148‑1 (MQ=255)
                                                                                                 aTGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAGAATCTGCCGCGTCCAAATGTTGATcgcg                                               >  1:3032887/1‑148 (MQ=255)
                                                                                                 aTGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATcgcg                                               >  2:3289069/1‑148 (MQ=255)
                                                                                                  tGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGa                                              >  1:1993306/1‑148 (MQ=255)
                                                                                                  tGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGa                                              <  2:1144550/148‑1 (MQ=255)
                                                                                                   gTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGAt                                             >  1:1992092/1‑148 (MQ=255)
                                                                                                   gTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGAt                                             >  1:1992917/1‑148 (MQ=255)
                                                                                                    tttAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATg                                            <  1:1776077/148‑1 (MQ=255)
                                                                                                    tttAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATg                                            <  1:431805/148‑1 (MQ=255)
                                                                                                    tttAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATg                                            <  1:1013511/148‑1 (MQ=255)
                                                                                                     ttAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATgc                                           >  1:26785/1‑148 (MQ=255)
                                                                                                      tAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATgcg                                          <  1:3024492/148‑1 (MQ=255)
                                                                                                           cTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGa                                     <  2:2179006/148‑1 (MQ=255)
                                                                                                             ggCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGata                                   <  2:709126/148‑1 (MQ=255)
                                                                                                                 cgACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCcg                               >  1:3245451/1‑148 (MQ=255)
                                                                                                                   aCCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCcgcg                             <  2:3297071/148‑1 (MQ=255)
                                                                                                                    ccTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGt                            <  1:496584/148‑1 (MQ=255)
                                                                                                                      ttCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTa                          <  1:1720557/148‑1 (MQ=255)
                                                                                                                      ttCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTa                          <  2:1799653/148‑1 (MQ=255)
                                                                                                                        cGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAATTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAAc                        <  2:320783/148‑1 (MQ=255)
                                                                                                                          cTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGt                      >  2:1364235/1‑148 (MQ=255)
                                                                                                                           tGTCCACTTCAGTGGTTTGACTTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTc                     <  2:1849966/148‑1 (MQ=255)
                                                                                                                           tGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTc                     <  1:1132373/148‑1 (MQ=255)
                                                                                                                           tGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTc                     <  2:1849990/148‑1 (MQ=255)
                                                                                                                              ccACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTgg                  >  2:16625/1‑148 (MQ=255)
                                                                                                                              ccACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTgg                  >  1:530508/1‑148 (MQ=255)
                                                                                                                                aCTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCt                <  1:1505223/148‑1 (MQ=255)
                                                                                                                                 cTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTg               >  2:825868/1‑148 (MQ=255)
                                                                                                                                      gTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATaa          <  1:1524723/148‑1 (MQ=255)
                                                                                                                                       tGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATaaa         >  2:1933499/1‑148 (MQ=255)
                                                                                                                                         gTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGa       <  2:1427667/148‑1 (MQ=255)
                                                                                                                                          tttGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGaa      <  2:394807/148‑1 (MQ=255)
                                                                                                                                          tttGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGaa      >  2:2069864/1‑148 (MQ=255)
                                                                                                                                          tttGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGaa      <  2:1464880/148‑1 (MQ=255)
                                                                                                                                           ttGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGaaa     >  1:1401093/1‑148 (MQ=255)
                                                                                                                                           ttGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGaaa     >  1:2324819/1‑148 (MQ=255)
                                                                                                                                              aaTTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGAAACCg  >  2:3082891/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CCGCCTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGAAACCG  >  CP009273/1003452‑1003741

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: