Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 802,772 T→G 8.0% 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.0% 184.3 / 7.4 63intergenic (+30/‑56)obg/proBGTPase Obg/glutamate 5‑kinase
Reads supporting (aligned to +/- strand):  ref base T (28/30);  new base G (3/2);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.44e-01

CCGTACCCTGCGTCGTACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAGGGCCTCGACCGCGGT  >  NC_002947/802625‑802917
                                                                                                                                                   |                                                                                                                                                 
ccGTACCCTGCGTCGTACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTg                                                                                                                                                  >  8:225195/1‑149 (MQ=255)
   tACCCTGCGTCGTACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAg                                                                                                                                               <  1:171973/149‑1 (MQ=255)
       cTGCGTCGTACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAgcggc                                                                                                                                           >  1:258896/1‑149 (MQ=255)
        tGCGTCGTACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAgcggcg                                                                                                                                          >  4:315440/1‑149 (MQ=255)
        tGCGTCGTACCGGCATCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGAGCGACTGATCGGGTGCAGTACACTAAACGCCGCTCTCTGGAgcgggg                                                                                                                                          >  8:297253/1‑147 (MQ=255)
             cgtACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCACGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGttttt                                                                                                                                     <  3:273150/149‑1 (MQ=255)
              gtACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGttttt                                                                                                                                     <  4:200163/148‑1 (MQ=255)
              gtACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGttttt                                                                                                                                     <  5:270845/148‑1 (MQ=255)
              gtACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTg                                                                                                                                    <  7:225195/149‑1 (MQ=255)
               tACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGATTTTg                                                                                                                                    >  5:160288/1‑148 (MQ=255)
                    gcatCAAGAGCGTGCATGCCATCGGCGCCGACGATGGTTGGGACGATGATTGCGAGGACCACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCa                                                                                                                              <  7:297253/146‑1 (MQ=255)
                    gCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCa                                                                                                                              <  3:303204/149‑1 (MQ=255)
                    gCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCa                                                                                                                              <  3:315440/149‑1 (MQ=255)
                       tCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACgg                                                                                                                           >  5:249730/1‑149 (MQ=255)
                        cAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGt                                                                                                                          <  6:310595/149‑1 (MQ=255)
                           gagCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATc                                                                                                                       >  8:30208/1‑149 (MQ=255)
                            agCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCt                                                                                                                      <  7:30208/149‑1 (MQ=255)
                                  cATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTAc                                                                                                                    <  2:241132/145‑1 (MQ=255)
                                  cATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTAc                                                                                                                    >  1:241132/1‑145 (MQ=255)
                                   aTGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCt                                                                                                                      <  5:175817/142‑1 (MQ=255)
                                   aTGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCt                                                                                                                      >  6:175817/1‑142 (MQ=255)
                                   aTGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGAtt                                                                                                               >  8:74004/1‑149 (MQ=255)
                                      aCATCGGGGACGACGATGGTTGGGACGAGGATTTCGTGGACGACGAAGACGGTGCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACTCTAAACGCCGCTGTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTct                                                                                                             <  4:114835/148‑2 (MQ=255)
                                      aCATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCg                                                                                                             <  5:255094/148‑1 (MQ=255)
                                        aTCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACa                                                                                                          >  1:212183/1‑149 (MQ=255)
                                              gacgacgaTGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGtt                                                                                                    <  8:195752/149‑1 (MQ=255)
                                              gacgacgaTGGTTGGGACGATGATTTAGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGAGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCAACAGTATCTACAGATTCGACaaaattt                                                                                                    >  4:246926/1‑143 (MQ=255)
                                               acgacgaTGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTg                                                                                                   <  2:122263/149‑1 (MQ=255)
                                                cgacgaTGGTTGTGACGATGATTTCGATGACGACGAAGACGGTCCGGAAATCATTTACGTGAGCGACTGATCGTTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTg                                                                                                   <  3:246926/148‑1 (MQ=255)
                                                        gTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGAt                                                                                          >  3:28371/1‑149 (MQ=255)
                                                                  tgatTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAg                                                                                <  2:212183/149‑1 (MQ=255)
                                                                   gatTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAgg                                                                               >  2:185567/1‑149 (MQ=255)
                                                                   gatTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAgg                                                                               >  3:301534/1‑149 (MQ=255)
                                                                    atTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGgaaga                                                                                              >  6:164191/1‑133 (MQ=255)
                                                                     tttCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGt                                                                              <  4:28371/148‑1 (MQ=255)
                                                                     tttCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTg                                                                             <  7:74004/149‑1 (MQ=255)
                                                                       tCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGgaaga                                                                                              >  3:220044/1‑130 (MQ=255)
                                                                       tCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGgaaga                                                                                              <  4:220044/130‑1 (MQ=255)
                                                                               gacgaAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTgg                                                                                                  <  8:234482/118‑1 (MQ=255)
                                                                                acgaAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTAAAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCaa                                                                  >  7:208726/1‑149 (MQ=255)
                                                                                   aaGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAgcg                                                               <  1:185567/149‑1 (MQ=255)
                                                                                        ggTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTggg                                                          <  4:301534/149‑1 (MQ=255)
                                                                                          tCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTc                                                        >  5:236349/1‑149 (MQ=255)
                                                                                               aaaTCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTaa                                                   <  6:237394/149‑1 (MQ=255)
                                                                                                            gcgcGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGct                                      >  4:227637/1‑149 (MQ=255)
                                                                                                               cgACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTctg                                   <  2:258896/149‑1 (MQ=255)
                                                                                                                gACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTctgc                                  >  1:156926/1‑149 (MQ=255)
                                                                                                                         gTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAg                                          >  1:10606/1‑132 (MQ=255)
                                                                                                                         gTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAg                                          <  2:10606/132‑1 (MQ=255)
                                                                                                                            gCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATg                      >  6:285873/1‑149 (MQ=255)
                                                                                                                             cAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATg                      <  8:208726/148‑1 (MQ=255)
                                                                                                                               gTACACTAAACGCCGCGCTCTGGAGCGGCGTTTTTGTATCCACGGGATCGACAGAGTCGACAGAGGTTGGAAGAAGATGCGAAGCAAGGTGCCGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCGCTGCTGACCGCCGATGGc                    >  5:160129/1‑148 (MQ=255)
                                                                                                                                   acTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAggg               >  4:209131/1‑149 (MQ=255)
                                                                                                                                   acTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCGGGGGCGGTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAggg               >  8:259046/1‑149 (MQ=255)
                                                                                                                                    aaaaaCGCCGCTCCAGAGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGAt                                                                                          <  6:154300/71‑1 (MQ=25)
                                                                                                                                    aaaaaCGCCGCTCCAGAGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGAt                                                                                          >  5:154300/3‑73 (MQ=25)
                                                                                                                                    aaaaaCGCCGCTCCAGAGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTa                                                    <  3:241692/109‑1 (MQ=39)
                                                                                                                                    aaaaaCGCCGCTCCAGAGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTa                                                    >  4:241692/3‑111 (MQ=39)
                                                                                                                                    aaaaaCGCCGCTCCAGAGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAGGGc              >  1:65214/3‑149 (MQ=255)
                                                                                                                                       aaCGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAGGGCCTc           <  6:236349/149‑1 (MQ=255)
                                                                                                                                          gccgcTCTCTGGAGCGGCGTTTTTGTATGCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAGGGCCTCGAc        >  1:180427/1‑149 (MQ=255)
                                                                                                                                          gccgcTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAGGGCCTCGAc        >  4:65340/1‑149 (MQ=255)
                                                                                                                                             gcTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAGGGCCTCGACcgc     <  6:249730/149‑1 (MQ=255)
                                                                                                                                                 tctGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAGGGCCTCGACCGCGGt  >  5:250929/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCGTACCCTGCGTCGTACCGGCGTCAAGAGCGTGCATGACATCGGCGACGACGATGGTTGGGACGATGATTTCGAGGACGACGAAGACGGTCCGGAAATCATTTACGTGCGCGACTGATCGGTTGCAGTACACTAAACGCCGCTCTCTGGAGCGGCGTTTTTGTATCCACGGTATCTACAGATTCGACATAGGTTGGAAGAAGATGCGAAGCAAGGTGACGGGCGCCAAGCGCTGGGTCGTTAAGATTGGCAGTGCTCTGCTGACCGCCGATGGCAAGGGCCTCGACCGCGGT  >  NC_002947/802625‑802917

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