Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,721,551 G→T 10.6% A238E (GCA→GAA)  glnQ ← glutamine transport ATP‑binding protein GlnQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,721,5510GT10.6% 195.0 / 12.7 66A238E (GCA→GAA) glnQglutamine transport ATP‑binding protein GlnQ
Reads supporting (aligned to +/- strand):  ref base G (39/20);  new base T (4/3);  total (43/23)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.37e-01

GCGGTTCAGCATCTGACTTCGCCTCGCCCTCTTCGCGGGCAAGCCCGCTCACACAGGTACATCACCCCCATCAAGGCAGTGGTGGACATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGAACAGCACGCGGTCAGCCACTTCGCGGGCAAAGCCCATTTCATGGGTGACGCAGACCATGGTCAT  >  NC_002947/5721404‑5721693
                                                                                                                                                   |                                                                                                                                              
gCGGTTCAGCATCTGACTTCGCCTCGCCCTCTTCGCGGGCAAGCCCGCTCACACAGGTACATCACCCCCATCAAGGCAGTGGTGGACATGCAGGAGAGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGc                                                                                                                                               >  4:69949/1‑149 (MQ=255)
    ttCAGCATCTGACTTCGCCTCGCCCTCTTCGCGGGCAAGCCCGCTCACACAGGTACATCACCCCCATCAAGGCAGTGGTGGACATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTgc                                                                                                                                           >  2:31841/1‑149 (MQ=255)
    ttCAGCATCTGACTTCGCCTCGCCCTCTTCGCGGGCAAGCCCGCTCACACAGGTACATCACCCCCATCAAGGCAGTGGTGGACATGCAGGAGCGGGCTTGCCCGCGAACAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTgc                                                                                                                                           >  1:149826/1‑149 (MQ=255)
     tCAGCATCTGACTTCGCCTCGCCCTCTTCGCGGGCAAGCCCGCTCACACAGGTACATCACCCCCATCAAGGCAGTGGTGGACATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTgcg                                                                                                                                          >  1:308641/1‑149 (MQ=255)
           tCTGACTTCGCCTCGCCCTCTTCGCGGGCAAGCCCGCTCACACAGGTACATCACCCCCATCAAGGCAGTGGTGGACATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCAcgcg                                                                                                                                    >  2:38478/1‑149 (MQ=255)
              gACTTCGCCTCGCCCTCTTCGCGGGCAAGCCCGCTCACACAGGTACATCACCCCCATCAAGGCAGTGGTGGACATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGCATTCCTgcgc                                                                                                                                         <  3:354294/141‑1 (MQ=255)
              gACTTCGCCTCGCCCTCTTCGCGGGCAAGCCCGCTCACACAGGTACATCACCCCCATCAAGGCAGTGGTGGACATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGCATTCCTgcgc                                                                                                                                         >  4:354294/1‑141 (MQ=255)
                      ctcgccCTCTTCGCGGGCAAGCCCGCTCACACAGGTACATCACCCCCATCAAGGCAGTGGTGGACATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGgc                                                                                                                         <  2:31090/149‑1 (MQ=255)
                       tcgccCTCTTCGCGGGCAAGCCCGCTCTCACAGGTACATCACCCCCATCAAGGCAGTGGTGGACAGGCAGGAGCGGGCTTGCCCGCGAAGACGCCCGCGCGGGTACAGCACCAGGCGCAGGAATGCCTGCGCACGCGGGGCGGTCGgcg                                                                                                                        >  8:257138/1‑149 (MQ=255)
                                   cgGGCAAGCCCGCTCACACAGGTACATCACCCCCATCAAGGCAGTGGTGGACATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATg                                                                                                            >  8:366736/1‑149 (MQ=255)
                                     ggCAAGCCCGCTCACACAGGTACATCACCCCCATCAAGGCAGTGGTGGACATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCg                                                                                                          <  7:48135/149‑1 (MQ=255)
                                      gCAAGCCCGCTCACACAGGTACATCACCCCCATCAAGGCAGTGGTGGACATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCgg                                                                                                         >  4:15954/1‑149 (MQ=255)
                                          gCCCGCTCACACAGGTACATCACCCCCATCAAGGCAGTGGTGGACATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGAGGTTACAGCACCTGGCGCAGGAATGCCTGAGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCgg                                                                                                     >  5:105159/1‑149 (MQ=255)
                                              gCTCACACAGGTACATCACCCCCATCAAGGCAGTGGTGGACATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGcc                                                                                                 >  2:385394/1‑149 (MQ=255)
                                                  acacaGGTACATCACCCCCATCAAGGCAGTGGTGGACATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAAt                                                                                             >  5:54548/1‑149 (MQ=255)
                                                  acacaGGTACATCACCCCCATCAAGGCAGTGGTGGACATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAAt                                                                                             >  1:18961/1‑149 (MQ=255)
                                                  acacaGGTACATCACCCCCATCAAGGCAGTGGTGGACATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCCGCGCAGCGAAGAATGCGGCCGGGGCCGAAt                                                                                             <  1:31841/149‑1 (MQ=255)
                                                   cacaGGTACATCACCCCCATCAAGGCAGTGGTGGACATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATc                                                                                            >  3:7281/1‑149 (MQ=255)
                                                    acaGGTACATCACCCCCATCAAGGCAGTGGTGGACATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCt                                                                                           <  5:91178/149‑1 (MQ=255)
                                                      aGGTACATCACCCCCATCAAGGCAGTGGTGGACATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTc                                                                                         <  1:385394/149‑1 (MQ=255)
                                                      aGGTACATCACCCCCATCAAGGCAGTGGTGGACATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTc                                                                                         <  2:308641/149‑1 (MQ=255)
                                                        gTACATCACCCCCATCAAGGCAGTGGTGGACATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCca                                                                                       <  7:65358/149‑1 (MQ=255)
                                                            aTCACCCCCATCAAGGCAGTGGTGGACATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCcagcag                                                                                   >  8:226192/1‑149 (MQ=255)
                                                             tCACCCCCATCAAGGCAGTGGTGGACATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCagcagc                                                                                  >  4:203263/1‑149 (MQ=255)
                                                             tCACCCCCATCAAGGCAGTGGTGGACATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCagcagc                                                                                  >  8:76378/1‑149 (MQ=255)
                                                                cccccATCAAGGCAGTGGTGGACATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTg                                                                               <  6:54548/149‑1 (MQ=255)
                                                                   ccATCAAGCCAGTGGTGGACATGCTGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCa                                                                            <  2:149826/149‑1 (MQ=255)
                                                                       cAAGGCAGTGGTGGACATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGCATCTTCCAGCAGCGTGCCATGGt                                                                        >  8:132433/1‑149 (MQ=255)
                                                                        aaGGCAGTGGTGGACATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTc                                                                       >  1:366015/1‑149 (MQ=255)
                                                                          ggCAGTGGTGGACATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCg                                                                      >  6:317733/1‑148 (MQ=255)
                                                                           gCAGTGGTGGACATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCg                                                                      >  3:32295/1‑147 (MQ=255)
                                                                                 gtggACATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGACGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGAACAg                                                              >  6:146592/1‑149 (MQ=255)
                                                                                    gACATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGAAc                                                                >  8:326332/1‑144 (MQ=255)
                                                                                    gACATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGAAc                                                                <  7:326332/144‑1 (MQ=255)
                                                                                    gACATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGAACAGCAc                                                           >  3:309168/1‑149 (MQ=255)
                                                                                      cATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGAACAGCAcgc                                                         <  7:76378/149‑1 (MQ=255)
                                                                                          cAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGAACAGCACGCGGt                                                      >  6:62719/1‑148 (MQ=255)
                                                                                                  ggCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGAACAGCCCGCGGTCAGGCACtt                                             >  7:317116/1‑149 (MQ=255)
                                                                                                   gCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGGTTGCCATGGTCGAAGAACAGCACGCGGTCAGCCACTTc                                            >  1:62067/1‑149 (MQ=255)
                                                                                                     ttGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGAACAGCACGCGGTCAGCCACTTcgc                                          >  8:53693/1‑149 (MQ=255)
                                                                                                     ttGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGAACAGCACGCGGTCAGCCACTTcgc                                          >  2:333102/1‑149 (MQ=255)
                                                                                                       gCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGAACAGCACGCGGTCAGCCACTTcgc                                          <  5:62719/147‑1 (MQ=255)
                                                                                                        cccGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGAACAGCACGCGGTCAGCCACTTCGCgg                                        <  3:15954/148‑1 (MQ=255)
                                                                                                         ccGCGAAGAAGCCCGCGCGGTTACCGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGCACCGCCCGCGGTCAGCCACTTCGCggg                                       >  5:210649/1‑148 (MQ=255)
                                                                                                         ccGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGAACAGCACGCGGTCAGCCACTTCGCggg                                       >  1:273517/1‑148 (MQ=255)
                                                                                                             gaagaaGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGAACAGCACGCGGTCAGCCACTTCGCgg                                        <  8:282794/143‑1 (MQ=255)
                                                                                                             gaagaaGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGAACAGCACGCGGTCAGCCACTTCGCgg                                        >  7:282794/1‑143 (MQ=255)
                                                                                                                   gCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCGGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGCTCGAAGAACAGCCCGCGGTCAGCCACGTCGCGGGCAAAGCCCAtt                            >  8:311176/1‑149 (MQ=255)
                                                                                                                       gcgcgGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGAACAGCACGCGGTCAGCCACTTCGCGGGCAAAGCCCATTTCAt                        <  5:146592/149‑1 (MQ=255)
                                                                                                                       gcgcgGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGAACAGCACGCGGTCAGCCACTTCGCGGGCAAAGCCCATTTCAt                        <  7:366736/149‑1 (MQ=255)
                                                                                                                        cgcgGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGAACAGCACGCGGTCAGCCACTTCGCGGGCAAAGCCCATTTCATg                       <  2:273517/149‑1 (MQ=255)
                                                                                                                         gcgGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGAACAGCACGCGGTCAGCCACTTCGCGGGCAAAGCCCATTTCATgg                      <  4:32295/149‑1 (MQ=255)
                                                                                                                           ggTTACAGCACCTGGCGCAGGCATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGACGAACAGCACGCGGTCAGCCACTTCGCGGGCAAAGCCCATTTCATGGGt                    <  7:257138/149‑1 (MQ=255)
                                                                                                                                 aGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGAACAGCACGCGGTCAGCCACTTCGCGGGCAAAGCCCATTTCATGGGTGACg                >  6:44855/1‑147 (MQ=255)
                                                                                                                                 aGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGAACAGCACGCGGTCAGCCACTTCGCGGGCAAAGCCCATTTCATGGGTGACg                >  3:253125/1‑147 (MQ=255)
                                                                                                                                 aGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGAACAGCACGCGGTCAGCCACTTCGCGGGCAAAGCCCATTTCATGGGTGACg                >  5:164566/1‑147 (MQ=255)
                                                                                                                                 aGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGAACAGAACGCGGTCAGCCACTTCGCGGGCAAAGCCCATTTCATGGGTGACGCa              >  7:22872/1‑149 (MQ=255)
                                                                                                                                   cACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGAACAGCACGCGGTCAGCCACTTCGCGGGCAAAGCCCATTTCATGGGTGACGCAGa            >  3:291626/1‑149 (MQ=255)
                                                                                                                                     ccTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGAACAGCACGCGGTCAGCCACTTCGCGGGCAAAGCCCATTTCATGGGTGACGCAGAc           >  4:42054/1‑148 (MQ=255)
                                                                                                                                       gtgcgcAGGCATTCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGAACAGCACGCGGTCAGCCACTTCGCGGGCAAAGCCCATTTCATGGGTGa                  >  5:298252/3‑139 (MQ=255)
                                                                                                                                       gtgcgcAGGCATTCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGAACAGCACGCGGTCAGCCACTTCGCGGGCAAAGCCCATTTCATGGGTGa                  <  6:298252/137‑1 (MQ=255)
                                                                                                                                        tgcgcAGGCATTCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGAACAGCACGCGGTCAGCCACTTCGCGGGCAAAGCCCATTTCATGGGTGACGCAGAcc          >  2:113724/2‑146 (MQ=255)
                                                                                                                                        tgcgcAGGCATTCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGAACAGCACGCGGTCAGCCACTTCGCGGGCAAAGCCCATTTCATGGGTGACGCAGAcc          <  1:113724/145‑1 (MQ=255)
                                                                                                                                        tgcgcAGGCATTCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGAACAGCACGCGGTCAGCCACTTCGCGGGCAAAGCCCATTTCATGGGTGACGCAGACCATg       >  7:124745/2‑149 (MQ=255)
                                                                                                                                         gcgcAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGAACAGCACGCGGTCAGCCACTTCGCGGGCAAAGCCCATTTCATGGGTGACGCAGACCATgg      >  7:124159/1‑149 (MQ=255)
                                                                                                                                             aGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGAACAGCACGCGGTCAGCCACTTCGCGGGCAAAGCCCATTTCATGGGTGACGCAGACCATGGTCAt  <  2:18961/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
GCGGTTCAGCATCTGACTTCGCCTCGCCCTCTTCGCGGGCAAGCCCGCTCACACAGGTACATCACCCCCATCAAGGCAGTGGTGGACATGCAGGAGCGGGCTTGCCCGCGAAGAAGCCCGCGCGGTTACAGCACCTGGCGCAGGAATGCCTGCGCACGCGGGTCTTTCGGCGCAGCGAAGAATGCGGCCGGGGCCGAATCTTCCAGCAGCTTGCCATGGTCGAAGAACAGCACGCGGTCAGCCACTTCGCGGGCAAAGCCCATTTCATGGGTGACGCAGACCATGGTCAT  >  NC_002947/5721404‑5721693

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: