Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,751,725 A→G 9.9% intergenic (‑163/‑308) 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,7250AG 9.9% 160.3 / 8.2 63intergenic (‑163/‑308)thiI/glnAtRNA sulfurtransferase/glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base A (30/24);  new base G (2/4);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.87e-01

CTGCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCAACGCCAAA  >  NC_002947/5751585‑5751870
                                                                                                                                            |                                                                                                                                                  
cTGCCGAAAATGAGGGGCGCGAATTAAATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGc                                                                                                                                            <  8:20911/149‑1 (MQ=255)
     gAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCTTTGCaccgc                                                                                                                                       >  4:89319/1‑145 (MQ=38)
     gAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGc                                                                                                                                       >  4:271897/1‑149 (MQ=255)
     gAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGc                                                                                                                                       <  1:233886/149‑1 (MQ=255)
     gAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGc                                                                                                                                       <  1:238369/149‑1 (MQ=255)
     gAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGc                                                                                                                                       >  6:289199/1‑149 (MQ=255)
     gAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCCTTCGCGGGTAAACCCGCTCCTACGGGCTTTGCaccgc                                                                                                                                       <  3:133978/149‑5 (MQ=38)
     gAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCCTTCGCGGGTAAACCCGCTCCTACGGGCTTTGCaccgc                                                                                                                                       <  3:129164/149‑5 (MQ=38)
      aaaaTGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGcc                                                                                                                                      >  1:130694/1‑149 (MQ=255)
       aaaTGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACGAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCGAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCCGTGGAAAGccc                                                                                                                                     >  3:39714/1‑149 (MQ=255)
          tGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTa                                                                                                                                  <  2:273568/149‑1 (MQ=255)
            aGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAg                                                                                                                                <  6:206377/148‑1 (MQ=255)
                  gcgAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCggg                                                                                                                         <  1:197765/149‑1 (MQ=255)
                  gcgAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCggg                                                                                                                         <  7:159524/149‑1 (MQ=255)
                  gcgAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCggg                                                                                                                         >  2:7400/1‑149 (MQ=255)
                        tataTCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTAcc                                                                                                                   >  6:65075/1‑149 (MQ=255)
                               gAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCcgcg                                                                                                               >  6:58960/1‑146 (MQ=255)
                                aaaTTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAaga                                                                                                           <  1:7400/149‑1 (MQ=255)
                                  aTTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCTTTGCACCGCCTGTAGGCGCGGGTTTACCCGCGaaaggg                                                                                                         <  3:89319/149‑5 (MQ=18)
                                     gCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGATAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGcc                                                                                                       <  1:257219/148‑1 (MQ=255)
                                        cAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCgg                                                                                                                          <  6:45351/126‑1 (MQ=255)
                                        cAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCgg                                                                                                                          >  5:45351/1‑126 (MQ=255)
                                        cAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAAACCGCGCCTACAGGCGGTGCAAATCCCGTAGGAGCGGGTTTACCCGCGAAGAGGccc                                                                                                      >  5:133292/1‑146 (MQ=255)
                                         aGGATTTGACCAACATTTGATCAGAAGCTTTTAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGTTTTACCCGCGAAGAGGCCCTAGa                                                                                                  <  4:39714/149‑1 (MQ=255)
                                               tGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTg                                                                                             <  6:272938/148‑1 (MQ=255)
                                                gACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTgg                                                                                            <  1:238293/148‑1 (MQ=255)
                                                gACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTgg                                                                                            <  5:65075/148‑1 (MQ=255)
                                                    aaCTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCccac                                                                                       >  5:243113/1‑149 (MQ=255)
                                                    aaCTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCacac                                                                                       <  5:289199/149‑1 (MQ=255)
                                                    aaCTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCACTTTCGAGGGTAAACCCCAGCCTACAGGCGGTTCAAAGCCCGAAGGAGCGGGTTTAACCGCGAAGAGGCCCTAGAACCAGGCacac                                                                                       >  3:26237/1‑149 (MQ=255)
                                                             gcAGAAGCTTTTAGTAAAATGCAGCGCCCAGGTTTTGCAGGGCCTGTTCCGGCCCTTTCGGGGGTAAACCCGCGCCTACTGGGGGTGAAAAGCCCGTCGGAGGGGTTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTcacgc                                                                              <  2:323710/148‑1 (MQ=255)
                                                                gAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCTTTGCACCGCCTGTAGGCGCgg                                                                                                                          <  3:142595/102‑1 (MQ=18)
                                                                gAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCTTTGCACCGCCTGTAGGCGCgg                                                                                                                          >  4:142595/1‑102 (MQ=18)
                                                                gAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTcacgccac                                                                           <  4:131700/149‑1 (MQ=255)
                                                                 aaGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTcacgccacg                                                                          >  6:155873/1‑149 (MQ=255)
                                                                          gaTAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTAtt                                                                 >  4:90351/1‑149 (MQ=255)
                                                                              aaTGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTg                                                             >  2:233048/1‑149 (MQ=255)
                                                                              aaTGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTg                                                             >  2:237518/1‑149 (MQ=255)
                                                                              aaTGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTg                                                             <  1:237518/149‑1 (MQ=255)
                                                                              aaTGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTg                                                             <  1:233048/149‑1 (MQ=255)
                                                                               aTGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGc                                                            <  2:209805/149‑1 (MQ=255)
                                                                                  cAGATGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGac                                                         <  6:126362/149‑1 (MQ=255)
                                                                                    gagGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGt                                                                                                                        >  8:33075/1‑84 (MQ=255)
                                                                                       gCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCAACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGAt                                                    >  2:80677/1‑149 (MQ=255)
                                                                                        ccAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGCCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATc                                                   <  3:271897/149‑1 (MQ=255)
                                                                                                         tGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAg                                  >  1:74881/1‑149 (MQ=255)
                                                                                                          gTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACCCCCCTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAgg                                 >  7:52883/1‑149 (MQ=255)
                                                                                                           ttCCGGCCCTTTCGCGGGTAAACCAGCGCCTAAAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGAGAAGAGGCCCTAGAACCTGGCACAACACTCACGCCACCCACCCTATTGGTGCGACAAGAACAAAAAACCCATAGTAAcgg                                >  2:211887/1‑149 (MQ=38)
                                                                                                              cGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGAGAACCTGCCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGc                             >  5:68622/1‑149 (MQ=255)
                                                                                                                  cttttCGGGGTAAAACCCGCGCGTACCGGGGGTCCAACCCCGTATGTTGGGGTTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAgg                                 <  5:40478/139‑1 (MQ=255)
                                                                                                                      tCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTc                       <  5:155873/147‑1 (MQ=255)
                                                                                                                              aaaCCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAAc             >  8:282297/1‑149 (MQ=255)
                                                                                                                               aaCCCGCGCCTACTGGCGGTGAAAAGCCGGTAGGTGCGGTTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAAcc            <  4:26237/149‑1 (MQ=255)
                                                                                                                                  ccGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCAcc                    >  1:35671/1‑138 (MQ=255)
                                                                                                                                  ccGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCAcc                    <  2:35671/138‑1 (MQ=255)
                                                                                                                                  ccGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCaa         >  6:18225/1‑149 (MQ=255)
                                                                                                                                  ccGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCaa         >  1:22609/1‑149 (MQ=255)
                                                                                                                                  ccGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCaa         >  4:150470/1‑149 (MQ=255)
                                                                                                                                  ccGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCaa         >  8:89205/1‑149 (MQ=255)
                                                                                                                                    gcgcCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCAACg       >  3:60942/1‑149 (MQ=255)
                                                                                                                                        cTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCa                                         >  1:60421/1‑111 (MQ=255)
                                                                                                                                        cTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCa                                         <  2:60421/111‑1 (MQ=255)
                                                                                                                                         tACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCAACGCCaaa  >  4:1569/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
CTGCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCAACGCCAAA  >  NC_002947/5751585‑5751870

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: