Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 802,773 G→A 8.0% intergenic (+31/‑55) 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,7730GA 8.0% 171.9 / 4.8 62intergenic (+31/‑55)obg/proBGTPase Obg/glutamate 5‑kinase
Reads supporting (aligned to +/- strand):  ref base G (27/30);  new base A (3/2);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01

CGTACCCTGCGTCGTACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGT‑GACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAGGGCCTCGACCGCG  >  NC_002947/802626‑802915
                                                                                                                                                   |                                                                                                                                               
cGTACCCTGCGTCGTACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTgg                                                                                                                                                <  3:208552/149‑1 (MQ=255)
  gtacgggCGTCGTACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAg                                                                                                                                              <  2:246435/143‑1 (MQ=255)
   aCCCTGCGTCGTACCGGCGTCAAGAGCGTGAATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAgc                                                                                                                                             >  3:103839/1‑149 (MQ=255)
       tGCGTCGTACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAgcggcg                                                                                                                                         >  8:172408/1‑149 (MQ=255)
         cgtcgtACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCCAGAGAGCGGCGt                                                                                                                                        <  1:211052/148‑1 (MQ=255)
          gtcgtACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGttt                                                                                                                                      >  2:57131/1‑149 (MQ=255)
            cgtACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGttttt                                                                                                                                    <  7:85497/149‑1 (MQ=255)
             gtACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGttttt                                                                                                                                    <  2:206690/148‑1 (MQ=255)
               aCCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTa                                                                                                                                 >  6:86266/1‑149 (MQ=255)
                ccGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTAt                                                                                                                                >  8:160660/1‑149 (MQ=255)
                 cGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATc                                                                                                                               >  5:85228/1‑149 (MQ=255)
                     gTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACg                                                                                                                           >  7:148821/1‑149 (MQ=255)
                     gTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCCAGAGAGCGGCGtttagt                                                                                                                                   <  1:73141/141‑4 (MQ=255)
                     gTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCCAGAGAGCGGCGtttagt                                                                                                                                   >  2:73141/1‑138 (MQ=255)
                      tCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACgg                                                                                                                          >  5:195568/1‑149 (MQ=255)
                        aaGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTa                                                                                                                        >  3:67462/1‑149 (MQ=255)
                         agagCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTAt                                                                                                                       <  1:57131/149‑1 (MQ=255)
                             cGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGATGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTAc                                                                                                                   <  1:38281/149‑1 (MQ=255)
                                gCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGt                                                                                                                         >  7:69150/1‑140 (MQ=255)
                                gCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGt                                                                                                                         >  7:65507/1‑140 (MQ=255)
                                gCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGt                                                                                                                         <  8:65507/140‑1 (MQ=255)
                                gCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGt                                                                                                                         <  8:69150/140‑1 (MQ=255)
                                 cATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGAt                                                                                                               >  8:134874/1‑149 (MQ=255)
                                     aCATCGGCGACGACGATGGTAGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTTTATACACCGTATATAAACATTCg                                                                                                            >  2:135421/1‑148 (MQ=255)
                                       aTCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACa                                                                                                         <  7:160660/149‑1 (MQ=255)
                                       aTCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTCATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCG‑TTTTGTATCCACGGTATCTACAGATTCGACa                                                                                                         >  8:59161/1‑148 (MQ=255)
                                          ggCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAg                                                                                                      <  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                                                                                                                         <  5:131983/109‑1 (MQ=255)
                                                               gatgatTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGt                                                                                                                         >  6:131983/1‑109 (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)
                                                                                                                                                   ggAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCa                                                                                 <  2:179060/65‑1 (MQ=255)
                                                                                                                                                   ggAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCa                                                                                 >  1:179060/1‑65 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CGTACCCTGCGTCGTACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGT‑GACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAGGGCCTCGACCGCG  >  NC_002947/802626‑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: ‑