Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 802,772 T→G 8.2% intergenic (+30/‑56) obg → / → proB GTPase Obg/glutamate 5‑kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947802,7720TG 8.2% 178.7 / 7.8 61intergenic (+30/‑56)obg/proBGTPase Obg/glutamate 5‑kinase
Reads supporting (aligned to +/- strand):  ref base T (27/29);  new base G (3/2);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.75e-01

CCGTACCCTGCGTCGTACCGGCGTCAAGAGCGTGCATGACATCGGC‑GACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGT‑GACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAGGGCCTCGACCGCG  >  NC_002947/802625‑802915
                                                                                                                                                    |                                                                                                                                                
ccGTACCCTGCGTCGTACCGGCGTCAAGAGCGTGCATGACATCGGCGGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGctctct                                                                                                                                                  >  1:59323/1‑149 (MQ=255)
 cGTACCCTGCGTCGTACCGGCGTCAAGAGCGTGCATGACATCGGC‑GACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTgg                                                                                                                                                <  3:208552/149‑1 (MQ=255)
   gtacgggCGTCGTACCGGCGTCAAGAGCGTGCATGACATCGGC‑GACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAg                                                                                                                                              <  2:246435/143‑1 (MQ=255)
    aCCCTGCGTCGTACCGGCGTCAAGAGCGTGAATGACATCGGC‑GACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAgc                                                                                                                                             >  3:103839/1‑149 (MQ=255)
        tGCGTCGTACCGGCGTCAAGAGCGTGCATGACATCGGC‑GACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAgcggcg                                                                                                                                         >  8:172408/1‑149 (MQ=255)
          cgtcgtACCGGCGTCAAGAGCGTGCATGACATCGGC‑GACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCCAGAGAGCGGCGt                                                                                                                                        <  1:211052/148‑1 (MQ=255)
           gtcgtACCGGCGTCAAGAGCGTGCATGACATCGGC‑GACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGttt                                                                                                                                      >  2:57131/1‑149 (MQ=255)
             cgtACCGGCGTCAAGAGCGTGCATGACATCGGC‑GACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGttttt                                                                                                                                    <  7:85497/149‑1 (MQ=255)
              gtACCGGCGTCAAGAGCGTGCATGACATCGGC‑GACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGttttt                                                                                                                                    <  2:206690/148‑1 (MQ=255)
                aCCGGCGTCAAGAGCGTGCATGACATCGGC‑GACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTa                                                                                                                                 >  6:86266/1‑149 (MQ=255)
                 ccGGCGTCAAGAGCGTGCATGACATCGGC‑GACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTAt                                                                                                                                >  8:160660/1‑149 (MQ=255)
                  cGGCGTCAAGAGCGTGCATGACATCGGC‑GACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATc                                                                                                                               >  5:85228/1‑149 (MQ=255)
                      gTCAAGAGCGTGCATGACATCGGC‑GACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACg                                                                                                                           >  7:148821/1‑149 (MQ=255)
                      gTCAAGAGCGTGCATGACATCGGC‑GACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCCAGAGAGCGGCGtttagt                                                                                                                                   <  1:73141/141‑4 (MQ=255)
                      gTCAAGAGCGTGCATGACATCGGC‑GACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCCAGAGAGCGGCGtttagt                                                                                                                                   >  2:73141/1‑138 (MQ=255)
                       tCAAGAGCGTGCATGACATCGGC‑GACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACgg                                                                                                                          >  5:195568/1‑149 (MQ=255)
                         aaGAGCGTGCATGACATCGGC‑GACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTa                                                                                                                        >  3:67462/1‑149 (MQ=255)
                          agagCGTGCATGACATCGGC‑GACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTAt                                                                                                                       <  1:57131/149‑1 (MQ=255)
                              cGTGCATGACATCGGC‑GACGACGATGGTTGGGACGATGATTTCGAGGACGATGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTAc                                                                                                                   <  1:38281/149‑1 (MQ=255)
                                 gCATGACATCGGC‑GACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGt                                                                                                                         >  7:69150/1‑140 (MQ=255)
                                 gCATGACATCGGC‑GACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGt                                                                                                                         >  7:65507/1‑140 (MQ=255)
                                 gCATGACATCGGC‑GACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGt                                                                                                                         <  8:65507/140‑1 (MQ=255)
                                 gCATGACATCGGC‑GACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGt                                                                                                                         <  8:69150/140‑1 (MQ=255)
                                  cATGACATCGGC‑GACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGAt                                                                                                               >  8:134874/1‑149 (MQ=255)
                                      aCATCGGC‑GACGACGATGGTAGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTTTATACACCGTATATAAACATTCg                                                                                                            >  2:135421/1‑148 (MQ=255)
                                        aTCGGC‑GACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACa                                                                                                         <  7:160660/149‑1 (MQ=255)
                                        aTCGGC‑GACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTCATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCG‑TTTTGTATCCACGGTATCTACAGATTCGACa                                                                                                         >  8:59161/1‑148 (MQ=255)
                                           ggC‑GACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAg                                                                                                      <  7:134874/149‑1 (MQ=255)
                                                acgacgaTGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTg                                                                                                  <  1:81132/149‑1 (MQ=255)
                                                 ctacgtTGGTTGGGACGATTATTTCTAGGACGACTAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGTTTGCAGTACACTAAACGCCGCTCTGTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTgg                                                                                                 <  1:135421/147‑1 (MQ=255)
                                                              gACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGCAGAAGATGCGaa                                                                                    <  2:234876/149‑1 (MQ=255)
                                                                 gatgatTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGt                                                                                                                         >  6:131983/1‑109 (MQ=255)
                                                                 gatgatTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGt                                                                                                                         <  5:131983/109‑1 (MQ=255)
                                                                 gatgatTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGc                                                                                  >  5:124492/1‑148 (MQ=255)
                                                                  atgatTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCaa                                                                                <  6:85228/149‑1 (MQ=255)
                                                                   tgatTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAg                                                                               <  6:195568/149‑1 (MQ=255)
                                                                         ctaGGACGACGAAGACGTTCCGGAAATCATTTTCGTGCGCGACTGATAGGTTGCAGTACACTATACGCGGGTCTCTGGAGCGGCGTTTTTGTATACACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGT‑GAc                                                                         <  4:103839/146‑1 (MQ=255)
                                                                         cGATGACGACGAAGACGGTCCGGAAATCATTTACGTTCGCTACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGa                                                                                                                                               <  7:139546/79‑1 (MQ=255)
                                                                         cGAGGACGAAGAAGACGGTCCGGAAATCATTTACGTGCTCGACTGATCGGTTGAAGTACACTAAACGCCCCTCTCTGGa                                                                                                                                               >  8:139546/1‑79 (MQ=255)
                                                                          gAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGT‑GACg                                                                        <  2:59323/148‑1 (MQ=255)
                                                                              acgacgAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGT‑GACGGgcgc                                                                   <  4:67462/149‑1 (MQ=255)
                                                                                  cgaAGACGGTCCGGAAAACAATTACCACCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTAGGAAGAAGATGCGAAGCAAGGT‑GACGGGCGCCAAg                                                               >  4:44316/1‑149 (MQ=255)
                                                                                           tCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGT‑GACGGGCGCCAAGCGCTGGGTc                                                      <  8:148821/149‑1 (MQ=255)
                                                                                               gAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCG‑TTTTGTATCCACTGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGT‑GACGGGCGCCAAGCGCTGGGTCGTTaa                                                 <  7:59161/149‑1 (MQ=255)
                                                                                                      tttACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGT‑GACGGGCGCCAAGCGCTGGGTCGTTAAGATTgg                                           <  5:211121/149‑1 (MQ=255)
                                                                                                        tACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTCGACGGGCGCCAAGCGCTGGGTCGTTAAGATTgg                                           <  1:32718/148‑1 (MQ=255)
                                                                                                            tGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACa                                                                                                                  <  3:164144/73‑1 (MQ=255)
                                                                                                            tGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACa                                                                                                                  >  4:164144/1‑73 (MQ=255)
                                                                                                              cgcgACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCTGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGT‑GACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGct                                    <  5:86266/148‑1 (MQ=255)
                                                                                                               gcgACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATGCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGT‑GACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGctct                                  >  1:259838/1‑149 (MQ=255)
                                                                                                               gcgACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGT‑GACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGctct                                  >  1:194557/1‑149 (MQ=255)
                                                                                                                cgACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGT‑GACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTctg                                 >  5:20104/1‑149 (MQ=255)
                                                                                                                      aTCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGT‑GACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGAcc                           <  2:92299/149‑1 (MQ=255)
                                                                                                                      aTCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGT‑GACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGAcc                           >  2:121074/1‑149 (MQ=255)
                                                                                                                         ggTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGT‑GACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGAccgcc                        >  7:179347/1‑149 (MQ=255)
                                                                                                                          gTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGT‑GACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGAccgcc                        >  7:51769/1‑148 (MQ=255)
                                                                                                                          gTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGT‑GACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGAccgcc                        >  8:67985/1‑148 (MQ=255)
                                                                                                                              cAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGT‑GACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATg                    <  3:44316/148‑1 (MQ=255)
                                                                                                                              cAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGT‑GACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATg                    <  6:20104/148‑1 (MQ=255)
                                                                                                                                   acaaaaaCGCCGCTCCAGAGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGT‑GACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAgg              >  2:170967/5‑149 (MQ=255)
                                                                                                                                     cTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGT‑GACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAGGGc            >  2:21233/1‑149 (MQ=255)
                                                                                                                                     aaaaaCGCCGCTCCAGAGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGT‑GACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAGGGc            >  7:15262/3‑149 (MQ=255)
                                                                                                                                               ctctctGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGT‑GACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAGGGCCTCGACcgcg  <  6:124492/149‑1 (MQ=255)
                                                                                                                                               ctctctGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGT‑GACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAGGGCCTCGACcgcg  <  1:121074/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CCGTACCCTGCGTCGTACCGGCGTCAAGAGCGTGCATGACATCGGC‑GACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGT‑GACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAGGGCCTCGACCGCG  >  NC_002947/802625‑802915

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: ‑