Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,751,719 G→T 6.9% intergenic (‑157/‑314) thiI ← / → glnA tRNA sulfurtransferase/glutamine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,751,7190GT 6.9% 162.4 / 3.9 58intergenic (‑157/‑314)thiI/glnAtRNA sulfurtransferase/glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base G (27/27);  new base T (2/2);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.78e-01

GCGCGCGCAGGGCCTGCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCAAC  >  NC_002947/5751572‑5751864
                                                                                                                                                   |                                                                                                                                                 
gcgcgtgcAGGGCCTGCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCgcgc                                                                                                                                                  >  1:239557/1‑149 (MQ=255)
    gtgcAGGGCCTGCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTAc                                                                                                                                              <  2:239557/147‑1 (MQ=255)
        aGGGCCTGCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGc                                                                                                                                          <  7:157120/149‑1 (MQ=255)
           gCCTGCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGt                                                                                                                                       <  5:352165/149‑1 (MQ=255)
           gCCTGCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGt                                                                                                                                       <  7:1474/149‑1 (MQ=255)
            ccTGCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTg                                                                                                                                      <  6:59210/149‑1 (MQ=255)
             cTGCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTg                                                                                                                                      >  7:119685/1‑148 (MQ=255)
                ccGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCa                                                                                                                                    >  6:264734/1‑147 (MQ=255)
                    aaaTGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGccc                                                                                                                              <  5:264734/149‑1 (MQ=255)
                      aTGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCACGt                                                                                                                            >  7:66484/1‑149 (MQ=255)
                       tGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTa                                                                                                                           <  6:226250/149‑1 (MQ=255)
                        tagcgggtgaaagtgcacccttacgatgcgttttttgatcgtgtcgcaccaatagggtgcgtggcgtgagtggtgtGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAg                                                                                                                          <  8:59538‑M1/73‑1 (MQ=255)
                              cgcgAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCgg                                                                                                                    <  6:330746/149‑1 (MQ=255)
                              cgcgAATTATATCGGAAAATGCTCAGGATTTGACCAACTTTTGATCAGAAGCTGTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCgg                                                                                                                    <  3:315122/149‑1 (MQ=255)
                                gaaagtgcacccttacgatgcgttttttgatcgtgtcgcaccaatagggtgcgtggcgtgagtggtgtGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGTTGCAAAGCCCGTAGGAGCGGGt                                                                                                                  >  7:59538‑M1/69‑149 (MQ=255)
                                 gAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCTTTGCACCGCCTGTAGGCGCGGGtt                                                                                                                 <  2:318112/149‑1 (MQ=25)
                                    gggcacccagacgatgcgttttttgatcgtgtcgcaccaatagggtgcgtggcgtgagtggtctGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTAc                                                                                                              >  8:360798‑M1/65‑149 (MQ=255)
                                         tCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCcgc                                                                                                          >  3:119440/1‑148 (MQ=255)
                                            gAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACCGGCGGTGCAAAGCCCGTATGAGCGGGTTTACCcgcg                                                                                                         >  1:23914/1‑146 (MQ=255)
                                               aTTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAgg                                                                                                   <  2:119216/149‑1 (MQ=255)
                                                ttccgcAGGATTTGACCAACTTTTGATCAGCAGCTTTCAGATAAATCCAGAGGCCAGGTTTTGCACCGCCTGTCCCGGCCCTTTCGCGGGTAACCCCGCGCCTACCGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                                  <  8:119685/144‑1 (MQ=37)
                                                ttGCTCAGGATTTGCCAAACTTTAGATCAGACTCGTTGAGATAAATGCTCAGGCCAGGTCTTGCCGCCCCTGTGCCGTCCCTTGCGCGGGTAAACCCGCGCCTACAGGCGGTGCAATGCCCGTAGCAGCGGGTTTACCCGCGAAGAGGc                                                                                                  <  8:66484/149‑1 (MQ=25)
                                                    tCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTAccc                                                                                                            >  3:286769/1‑135 (MQ=255)
                                                    tCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTAccc                                                                                                            <  4:286769/135‑1 (MQ=255)
                                                     cAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGccc                                                                                                >  2:26638/1‑146 (MQ=255)
                                                     cAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGccc                                                                                                >  5:311226/1‑146 (MQ=255)
                                                             gACCAAGTTTTGATCATAAGATTTGAGAAAATTGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCTTTTCGCGGTTAACCCCGGGCCTACTGGCGGTGAAAAGCCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTgg                                                                                      <  2:376513/148‑1 (MQ=37)
                                                             gACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTgg                                                                                      <  1:26638/148‑1 (MQ=255)
                                                             gACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTAAAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCAAGAACCTgg                                                                                      >  1:376513/1‑148 (MQ=255)
                                                              aCCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTTAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGa                                                                                                                        >  1:123267/1‑113 (MQ=255)
                                                              aCCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGa                                                                                                                        <  2:123267/113‑1 (MQ=255)
                                                                cAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGcaca                                                                                  <  5:26072/149‑1 (MQ=255)
                                                                            aGAAGCTTTTATATAAATGCAGAGGCCAGGTGTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGGGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTcacgcc                                                                       <  8:112419/148‑1 (MQ=255)
                                                                             gAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTcacgccac                                                                     >  4:183561/1‑149 (MQ=255)
                                                                              aaGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCAcacca                                                                               >  8:319514/1‑138 (MQ=255)
                                                                              aaGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCAcacca                                                                               <  7:319514/138‑1 (MQ=255)
                                                                                   ttGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCTTTGCACCGCCTGTAGGCGCggg                                                                                                                   <  5:275870/97‑1 (MQ=18)
                                                                                   ttGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCTTTGCACCGCCTGTAGGCGCggg                                                                                                                   >  6:275870/1‑97 (MQ=18)
                                                                                   ttGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCCCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGGTTTCCCGCGAAGAGGCCCTAGAACCTGGCCAACCCCTCACgccccgcccc                                                                >  4:364814/1‑148 (MQ=38)
                                                                                    tGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCAccc                                                               <  6:171573/148‑1 (MQ=255)
                                                                                     gagaTAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCAcccaa                                                             >  3:351054/1‑147 (MQ=255)
                                                                                          aaaTGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGt                                                        >  8:365601/1‑149 (MQ=255)
                                                                                            aTGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGc                                                      <  4:351054/149‑1 (MQ=255)
                                                                                               cAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAATCCCGCGCCTACAGGCGGTGCAAAGc                                                                                                                                <  4:341528/72‑1 (MQ=255)
                                                                                               cAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGc                                                                                                                                >  3:341528/1‑72 (MQ=255)
                                                                                               cAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGa                                                    >  1:1143/1‑148 (MQ=255)
                                                                                                agagGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAAACCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGAGGGTTTACCCGCGCAGAGGCCCGATACCCTGGCACACCACTCCCGCCACGCCCCCTATTGGTGCGaca                                                  >  3:41122/1‑149 (MQ=37)
                                                                                                  agGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGacac                                                 >  6:275918/1‑148 (MQ=255)
                                                                                                             tGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAAc                                      <  3:183561/148‑1 (MQ=255)
                                                                                                               cAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCa                                   >  7:89101/1‑149 (MQ=255)
                                                                                                                aGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACCCCCCTCACGCCCCGCACCCTATTGGTGCGACACGATCAAAAAACGCAt                                  >  1:144650/1‑149 (MQ=255)
                                                                                                                aGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCGATTGGTGCGACACGATCAAAAAACGCAt                                  >  2:94260/1‑149 (MQ=255)
                                                                                                                     cTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACACGCCGTGCAACCCCCGTAGGAGCCGGTTAACCCGCGAAGAGGCCCGACCACCTCCCACCCCACGCCCGCCCCGCACCCAAATGGTGCGCCCCGCTCAAAAAACGCCTCGTaa                             >  8:167669/1‑149 (MQ=16)
                                                                                                                        ttCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAGACGCATCGTAAgg                           <  4:119440/148‑1 (MQ=255)
                                                                                                                           cGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTg                        <  7:365601/148‑1 (MQ=255)
                                                                                                                               ccTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACtt                   <  2:23914/149‑1 (MQ=255)
                                                                                                                                      cgGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCg            >  3:341426/1‑149 (MQ=255)
                                                                                                                                              cccGCTCCTACGGGCTTTGCACCGCCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAAccc     >  2:124658/1‑148 (MQ=37)
                                                                                                                                               ccGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTAATGGTGCGACACGATCAAAAAACCCATCGTAAGGGACCACTTACACCCGCTAACCCa    >  7:145349/1‑148 (MQ=255)
                                                                                                                                                cgctCCTACGGGCTTTGCACCGCCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTAACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCAAc  >  2:389059/1‑149 (MQ=37)
                                                                                                                                                   |                                                                                                                                                 
GCGCGCGCAGGGCCTGCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCAAC  >  NC_002947/5751572‑5751864

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