Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,751,725 A→G 7.1% 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 7.1% 157.4 / 2.8 59intergenic (‑163/‑308)thiI/glnAtRNA sulfurtransferase/glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base A (23/28);  new base G (2/2);  total (26/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-01

GGGCCTGCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAG‑‑GCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCAACGC  >  NC_002947/5751581‑5751866
                                                                                                                                                |                                                                                                                                               
gggCCTGCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCg                                                                                                                                             <  7:291208/149‑1 (MQ=255)
gggCCTGCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCg                                                                                                                                             <  6:177023/149‑1 (MQ=255)
  gCCTGCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGt                                                                                                                                           <  2:10805/149‑1 (MQ=255)
    cTGCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGc                                                                                                                                         <  5:9329/149‑1 (MQ=255)
    cTGCCGAAAAACAGCGGCGCGAATAATATCGCACATTGCGCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGc                                                                                                                                         <  7:274608/149‑1 (MQ=38)
           aaaTGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGccc                                                                                                                                  <  8:129164/149‑1 (MQ=255)
           aaaTGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGccc                                                                                                                                  <  1:202182/149‑1 (MQ=255)
           aaaTGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGccc                                                                                                                                  <  2:86555/149‑1 (MQ=255)
            aaTGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCg                                                                                                                                 <  6:162407/149‑1 (MQ=255)
             aTGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCg                                                                                                                                 >  8:71638/1‑148 (MQ=255)
             aTGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGt                                                                                                                                <  5:147203/149‑1 (MQ=255)
                     cgcgAATTATAGCGGAAATTGCTCAGGATATGACCAACTTTTGATCAGCAGCTTTGAGATCAATGCAGAGGCCAGGTTGTGCCGCGCCGGTTCCGGCCTCTTCGCGGGTAAACCCGCGCCTACGGGCTTTGCCCCGCCTGTAGGCGCgg                                                                                                                        <  5:258615/149‑1 (MQ=21)
                      gcgAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCggg                                                                                                                       <  7:71638/149‑1 (MQ=255)
                      gcgAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCggg                                                                                                                       >  4:305890/1‑149 (MQ=255)
                            tataTCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCTTTGCACCGCCTGTAGgcgc                                                                                                                          >  4:337576/1‑140 (MQ=37)
                            tataTCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCTTTGCACCGCCTGTAGgcgc                                                                                                                          <  3:337576/140‑1 (MQ=37)
                              tatCGGAAATTCCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGTAGAGGCCAGGTTTTGCAGCGCCTGATCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCcg                                                                                                               <  8:86710/149‑1 (MQ=255)
                                tCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCcgc                                                                                                              <  8:244218/148‑1 (MQ=255)
                                  ggAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGACGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACCGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACTCGCGaa                                                                                                           <  7:314917/149‑1 (MQ=255)
                                   gAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCcgcg                                                                                                             >  3:50126/1‑146 (MQ=255)
                                    aaaTTGCTCAGGATTTTACAAACTTTTGATCACAAGCTTGTAGATAAATGCAGAGGCCCGGTTTGGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGGGGGTTAACCCGCGAaga                                                                                                         <  6:54118/149‑1 (MQ=38)
                                    aaaTTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAaga                                                                                                         <  6:154748/149‑1 (MQ=255)
                                       ttGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAG‑‑Gc                                                                                                    >  8:274832/1‑149 (MQ=255)
                                           tCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAG‑‑Gccc                                                                                                  >  2:22108/1‑147 (MQ=255)
                                                aTTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAG‑‑GCCCTAGAAcc                                                                                           >  6:80229/1‑149 (MQ=255)
                                                 tttGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAG‑‑GCCCTAGAACCt                                                                                          <  2:216310/149‑1 (MQ=255)
                                                     aCCAACTTATCATCAGAAGCTTTGACAAAAATGCAGAGGCCAGGTTTTGCAGCCCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAG‑‑GCCCTAGAACCTGGca                                                                                      <  4:305680/149‑1 (MQ=38)
                                                       cAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAG‑‑GCCCt                                                                                                 >  5:117633/1‑136 (MQ=255)
                                                       cAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAG‑‑GCCCt                                                                                                 <  6:117633/136‑1 (MQ=255)
                                                       cAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAG‑‑GCCCTAGAACCTGGcaca                                                                                    >  5:376672/1‑149 (MQ=255)
                                                       cAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGACTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGGAAAGCCCGTAGGAGCGGGTTTACACGCGAAGAG‑‑GCCCTAGACCCTGGcac                                                                                     >  8:245577/1‑148 (MQ=255)
                                                          cTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAG‑‑GCCCTAGAACCTGGCAcacca                                                                                 <  7:114323/149‑1 (MQ=255)
                                                             ttGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAG‑‑GCCCTAGAACCTGGCACACCACt                                                                               >  3:346749/1‑148 (MQ=255)
                                                                         tttGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAG‑‑GCCCTAGAACCTGGCACACCACTCACGCCACGCAcc                                                                  >  1:100645/1‑149 (MQ=255)
                                                                         tttGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAG‑‑GCCCTAGAACCTGGCACACCACTCACGCCACGCAcc                                                                  >  3:270968/1‑149 (MQ=255)
                                                                           tGAGATAATGGCAGAGGCCCCGTTTTGCCGCGCCTGTCCCGGCCCTTTCGCGGGTAAACCCGCGCATACAGGCGGTGCAAAGCCCTTAGGAGCGGGTTTACCCGCGAAGAG‑‑GCCCTAGAACCTGGCACACCACTCACGCCACGCAcc                                                                  <  8:328371/147‑1 (MQ=37)
                                                                                 aaaTGCAGCGGCCAGTTTTTGCAGCGCGTGTCCCGGCCTTTGCGGGGGAAACCCCGCGCCTACCGGCGGTGAAACGCCGGTAGGCGCGGTTTAACCCGCGAAGAG‑‑GCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGt                                                          <  1:183387/149‑1 (MQ=255)
                                                                                           gCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAG‑‑GCCCTAGACCCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGAt                                                >  7:79533/1‑149 (MQ=255)
                                                                                           gCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAG‑‑GCCCTAGAACCTGGCAcacc                                                                                  <  8:29216/115‑1 (MQ=255)
                                                                                           gCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAG‑‑GCCCTAGAACCTGGCAcacc                                                                                  >  7:29216/1‑115 (MQ=255)
                                                                                           gCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAG‑‑GCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGAt                                                >  3:268966/1‑149 (MQ=255)
                                                                                           gCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTAACCGCGAAGAGTAAACGTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGacac                                                   >  1:123571/1‑148 (MQ=38)
                                                                                            ccAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAG‑‑GCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATc                                               <  4:346749/149‑1 (MQ=255)
                                                                                            ccAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGGTTACCCGCGAAGAG‑‑GCCCGAGAACCTGGCCCCCCGCTGCCGCCACGCCCCCTATTGGTGCGACACGag                                                >  4:60173/1‑147 (MQ=38)
                                                                                                        gcgcCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAG‑‑GCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATc                                   <  1:22108/149‑1 (MQ=255)
                                                                                                           ccTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGGAGGAGCGGGTTTACCCGCGAAGAG‑‑GCCCTAGAACCTGGCACACCACACCCGCCCCGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGt                                 >  2:230366/1‑148 (MQ=255)
                                                                                                            cTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACCGGCGGTGCAAAGCCCCTAGGAGCGGGTTTACCCGCGAAGAG‑‑GCCCGAGAACCTGGCACACCACTCACGCCACGCACCCAATTGCTGCGACACGATCAAAAACCGCATCGt                                 >  6:77768/1‑147 (MQ=38)
                                                                                                            cTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAG‑‑GCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTaa                               <  3:305890/149‑1 (MQ=255)
                                                                                                              gTTCCGGCCCTTTCGCGGGTAAACCCGCGCCGACAGGCGGTGCCAAGCCCGTAGGAGCGGGTTTACCCGGGAAGAG‑‑GCCCTTGAACCTGGCACACCCCGACCGCCACGCCCCCTATTGGTGGGACACGCTCAAAAACCGCATCGGAAgg                             >  8:222004/1‑149 (MQ=37)
                                                                                                                          tCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAG‑‑GCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTc                   <  4:50126/147‑1 (MQ=255)
                                                                                                                               ggTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAG‑‑GCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGc             >  4:368029/1‑148 (MQ=255)
                                                                                                                                   aaCCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAG‑‑GCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAAc         <  4:270968/148‑1 (MQ=255)
                                                                                                                                   aaCCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAG‑‑GCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAAc         <  5:80229/148‑1 (MQ=255)
                                                                                                                                      ccGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAG‑‑GCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTACGGGTGCACTTTCACCCGCTACCCCaa     >  7:33162/1‑149 (MQ=255)
                                                                                                                                      ccGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAG‑‑GCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCaa     >  3:67059/1‑149 (MQ=255)
                                                                                                                                       ggccggaaCAGGCGCTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAG‑‑GCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCAAc    <  3:172410/142‑1 (MQ=255)
                                                                                                                                       cgcgCCTACCGGGGGTGCAAAGCCCGTAGTAGCGGGTTTACCCGCGAAGAG‑‑GCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCTACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACTCGCTAACCCAAc    <  5:62694/149‑1 (MQ=255)
                                                                                                                                        gctccTACGGGCTTTGCACCGCCTGTAGGCGCGGGTTTACCCGCGAAAGG‑‑GCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCCTAAGGGTGCACTTTCACCCGCTAACCCAACg   >  7:255986/4‑149 (MQ=25)
                                                                                                                                          gcCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTATACCCGCGAAGAG‑‑GCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCAACGc  >  3:257025/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                               
GGGCCTGCCGAAAATGAGGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAG‑‑GCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCAACGC  >  NC_002947/5751581‑5751866

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: