Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,751,719 G→T 9.6% 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 9.6% 179.7 / 7.2 60intergenic (‑157/‑314)thiI/glnAtRNA sulfurtransferase/glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base G (28/26);  new base T (2/4);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.39e-01

GCGCGCAGGGCCTGCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCAACG  >  NC_002947/5751574‑5751865
                                                                                                                                                 |                                                                                                                                                   
gcgcgcAGGGCCTGCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGcc                                                                                                                                                   >  3:131700/1‑148 (MQ=255)
           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                                                                                                                    >  2:7400/1‑149 (MQ=255)
                             gcgAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCggg                                                                                                                    <  7:159524/149‑1 (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                                                                                                                     >  5:45351/1‑126 (MQ=255)
                                                   cAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCgg                                                                                                                     <  6:45351/126‑1 (MQ=255)
                                                   cAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAAACCGCGCCTACAGGCGGTGCAAATCCCGTAGGAGCGGGTTTACCCGCGAAGAGGccc                                                                                                 >  5:133292/1‑146 (MQ=255)
                                                    aGGATTTGACCAACATTTGATCAGAAGCTTTTAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGTTTTACCCGCGAAGAGGCCCTAGa                                                                                             <  4:39714/149‑1 (MQ=255)
                                                          tGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTg                                                                                        <  6:272938/148‑1 (MQ=255)
                                                           gACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTgg                                                                                       <  5:65075/148‑1 (MQ=255)
                                                           gACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTgg                                                                                       <  1:238293/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                                                                                                                     >  4:142595/1‑102 (MQ=18)
                                                                           gAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCTTTGCACCGCCTGTAGGCGCgg                                                                                                                     <  3:142595/102‑1 (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                                                        <  1:237518/149‑1 (MQ=255)
                                                                                         aaTGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTg                                                        <  1:233048/149‑1 (MQ=255)
                                                                                         aaTGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTg                                                        >  2:237518/1‑149 (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    >  8:89205/1‑149 (MQ=255)
                                                                                                                                             ccGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCaa    >  6:18225/1‑149 (MQ=255)
                                                                                                                                             ccGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCaa    >  1:22609/1‑149 (MQ=255)
                                                                                                                                             ccGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCaa    >  4:150470/1‑149 (MQ=255)
                                                                                                                                               gcgcCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCAACg  >  3:60942/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCGCGCAGGGCCTGCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCAACG  >  NC_002947/5751574‑5751865

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: