Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,651,210 T→C 6.2% intergenic (‑18/+195) PP_4959 ← / ← fba diguanylate cyclase/fructose‑bisphosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,651,2100TC 6.2% 187.5 / 2.3 64intergenic (‑18/+195)PP_4959/fbadiguanylate cyclase/fructose‑bisphosphate aldolase
Reads supporting (aligned to +/- strand):  ref base T (29/31);  new base C (2/2);  total (31/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 = 9.67e-01

GCCGTCAATACTTGCCAATCCTGCCGGGCCAGCAAGGCCCGCATGCTGATGAGGTTCTCGGGGTAGTCGTCTACCACCAGAAGAACTGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCC  >  NC_002947/5651062‑5651357
                                                                                                                                                      |                                                                                                                                                   
gCCGTAAATACTTGCCAATCCTGCATGGCCAGCAAGTCCCGCATGCTGATGAGGTTCTCGGGGTAGTCTTCTACCACCAGAACAACTGAGCTCCCATCCAGTGGCTGATG‑‑TTGAGGGTAAGAGCCATCCATGCTGATTCTCTGTTAtgt                                                                                                                                                     >  4:150199/1‑149 (MQ=255)
    tCAATACTTGCCAATCCTGCCGGGCCAGCAAGGCCCGCATGCTGATGAGGTTCTCGGGGTAGTCGTCTACCACCAGAATAACTGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGAcc                                                                                                                                                 >  5:258676/1‑149 (MQ=255)
    tCAATACTTGCCAATCCTGCCGGGCCAGCAAGGCCCGCATGCTGATGAGGTTCTCGGGGTAGTCGTCTACCACCAGAAGAACTGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGccc                                                                                                                                                 >  8:351426/1‑149 (MQ=255)
       aTACTTGCAAATCCTGCAGGGACAGCAAGGCCCGCATGCTGATGAGGTTCTAGGGGTAGTCGTCTACCACCAGAAGAACTGAGCTGCCATACAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCAAGCTGCTTATCTGTTATCTGACCGAc                                                                                                                                              >  8:258927/1‑149 (MQ=255)
                 aTCCTTCCGGGCCAGCAAGGCCCGCATGCTGATGAGGTTCTCGGGGTAGTCGTCTACCACCAGAAGAACTGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGCGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATAc                                                                                                                                    <  7:258927/149‑1 (MQ=255)
                     tGCCGGGCCAGCAAGGCCCGCATGCTGATGAGGTTCTCGGGGTAGTCGTCTACCACCAGAAGAACTGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGtt                                                                                                                                >  5:89202/1‑149 (MQ=255)
                       ccGGGCCAGCAAGGCCCGCATGCTGATGAGGTTCTCGGGGTAGTCGTCTACCACCAGAAGAACTGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGtt                                                                                                                                <  4:203485/147‑1 (MQ=255)
                       ccGGGCCAGCAAGGCCCGCATGCTGATGAGGTTCTCGGGGTAGTCGTCTACCACCAGAAGAACTGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTgg                                                                                                                              <  8:4279/149‑1 (MQ=255)
                       ccGGGCCAGCAAGGCCCGCATGCTGATGAGGTTCTCGGGGTAGTCGTCTACCACCAGAAGAACAGAGCTGCCATCCAGTGGCTGCTGTTTTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGtt                                                                                                                                <  1:106649/149‑1 (MQ=255)
                        cGGGCCAGCAAGGCCCGCATGCTGATGAGGTTCTCGGGGTAGTCGTCTACCACCAGAAGAACTGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGa                                                                                                                             >  4:300566/1‑149 (MQ=255)
                           gCCAGCAAGGCCCGCATGCTGATGAGGTTCTCGGGGTAGTCGTCTACCACCAGAAGAACTGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATc                                                                                                                           >  6:62916/1‑148 (MQ=255)
                            ccAGCAAGGCCCGCATGCTGATGAGGTTCTCGGGGTAGTCGTCTACCACCAGAAGAACTGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCaaa                                                                                                                                        >  8:104856/1‑134 (MQ=255)
                             cAGCAAGGCCCGCATGCTGATGAGGTTCTCGGGGTAGTCGTCTACCACCAGAAGAACTGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCAt                                                                                                                        <  2:335876/149‑1 (MQ=255)
                               gCAAGGCCCGCATGCTGATGAGGTTCTCGGGGTAGTCGTCTACCACCAGAAGAACTGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTa                                                                                                                      <  6:35153/149‑1 (MQ=255)
                               gCAAGGCCCCCATGCTGATGAGGTTCTCGGGGTACTAGTCTACCACCAGAAGAACAGAGCTGCCATCCAGTGCCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGATGCATCTCTGTTATCTGACCGACTTCAAAATAAAGTTGGATCCATTa                                                                                                                      >  8:60032/1‑149 (MQ=255)
                                 aaGGCCCGCATGCTGATGAGGTTCTCGGGGTAGTCGTCTACCACCAGAAGAACTGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTg                                                                                                                               <  5:1805/138‑1 (MQ=255)
                                 aaGGCCCGCATGCTGATGAGGTTCTCGGGGTAGTCGTCTACCACCAGAAGAACTGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTg                                                                                                                               >  6:1805/1‑138 (MQ=255)
                                 aaGGCCCGCATGCTGATGAGGTTCTCGGGGTAGTCGTCTACCACCAGAAGAACTGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAg                                                                                                                     <  5:47902/148‑1 (MQ=255)
                                  aGGCCCGCATGCTGATGAGGTTCTCGGGGTAGTCGTCTACCACCAGAAGAACTGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTa                                                                                                                   <  7:153329/149‑1 (MQ=255)
                                     cccGCATGCTGATGAGGTTCTCGGGGTAGTCGTCTACCACCAGAAGAACTGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTActc                                                                                                                <  6:89202/149‑1 (MQ=255)
                                      ccGCATGCTGATGAGGTTCTCGGGGTAGTCGTCTACCACCAGAAGAACTGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTActct                                                                                                               <  5:361796/149‑1 (MQ=255)
                                      ccGCATGCTGATGAGGTTCTCGGGGTAGTCGTCTACCACCAGAAGAACTGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTActc                                                                                                                >  1:71465/1‑148 (MQ=255)
                                            gCTGATGAGGTTCTCGGGGTAGTCGTCTACCACCAGAAGAACTGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGAc                                                                                                                                              >  4:199337/1‑112 (MQ=255)
                                            gCTGATGAGGTTCTCGGGGTAGTCGTCTACCACCAGAAGAACTGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGAc                                                                                                                                              <  3:199337/112‑1 (MQ=255)
                                              tgatgaGGTTCTCGGGGTAGTCGTCTACCACCAGAAGAACTGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTActct                                                                                                               <  8:150515/141‑1 (MQ=255)
                                              tgatgaGGTTCTCGGGGTAGTCGTCTACCACCAGAAGAACTGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTActct                                                                                                               >  7:150515/1‑141 (MQ=255)
                                                atgaGGTTCTCGGGGTAGTCGTCTACCACCAGAAGAACTGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGa                                                                                                       <  7:210419/147‑1 (MQ=255)
                                                 tgaGGTTCTCGGGGTAGTCGTCTACCACCAGAAGAACTGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAgg                                                                                                     <  3:245956/148‑1 (MQ=255)
                                                     gTTCTCGGGGTAGTCGTCTACCACCAGAAGAACTGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGc                                                                                                 <  3:300566/148‑1 (MQ=255)
                                                              gTAGTCGTCTACCACCAGAAGAACTGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACgcaat                                                                                       <  7:351426/149‑1 (MQ=255)
                                                                  tcgtcTACCACCAGAAGAACTGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACgcaatgcaa                                                                                   <  2:71465/149‑1 (MQ=255)
                                                                     tcTACCACCAGAAGAACTGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCa                                                                                <  5:62916/149‑1 (MQ=255)
                                                                           accaGAAGAACTGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGc                                                                          >  4:267562/1‑149 (MQ=255)
                                                                              agaagaACTGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTAt                                                                       >  8:167235/1‑149 (MQ=255)
                                                                                    aCTGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGa                                                                 <  6:258676/149‑1 (MQ=255)
                                                                                      tGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGAt                                                                >  7:159155/1‑148 (MQ=255)
                                                                                           tGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGCGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACg                                                                                           >  6:372313/1‑116 (MQ=255)
                                                                                           tGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGCGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACg                                                                                           <  5:367331/116‑1 (MQ=255)
                                                                                           tGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGCGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACg                                                                                           >  6:367331/1‑116 (MQ=255)
                                                                                           tGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGCGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACg                                                                                           <  5:372313/116‑1 (MQ=255)
                                                                                            gCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTActct                                                                                                               >  6:373273/1‑95 (MQ=255)
                                                                                               aTCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTAc                                                      >  6:188437/1‑149 (MQ=255)
                                                                                                tCCATTGGCTTCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTTGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACg                                                     <  3:150199/149‑1 (MQ=255)
                                                                                                tCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCt                                                                                                         <  1:240977/97‑1 (MQ=255)
                                                                                                tCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCt                                                                                                         >  2:240977/1‑97 (MQ=255)
                                                                                                    gTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCaa                                                                                               <  2:5797/103‑1 (MQ=255)
                                                                                                    gTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCaa                                                                                               >  1:5797/1‑103 (MQ=255)
                                                                                                         tgctg‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGAc                                            >  2:157088/1‑149 (MQ=255)
                                                                                                         tgctg‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATAAGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGAc                                            >  2:135770/1‑149 (MQ=255)
                                                                                                                     ggTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAAcggtcgg                                  >  6:75032/1‑149 (MQ=255)
                                                                                                                     ggTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAAcggtcgg                                  >  2:324480/1‑149 (MQ=255)
                                                                                                                      gTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAAcggtcggt                                 <  7:167235/149‑1 (MQ=255)
                                                                                                                       tAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACgcaatgcaat                                                                                  >  3:289819/1‑99 (MQ=255)
                                                                                                                       tAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACgcaatgcaat                                                                                  <  4:289819/99‑1 (MQ=255)
                                                                                                                         aGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTa                                                                                                                      >  2:122403/1‑61 (MQ=255)
                                                                                                                         aGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTa                                                                                                                      <  1:122403/61‑1 (MQ=255)
                                                                                                                              catccatGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCgg                          >  3:23738/1‑148 (MQ=255)
                                                                                                                                    tgctgcTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATg                   <  3:267562/149‑1 (MQ=255)
                                                                                                                                       tgctTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGc                             <  6:215748/136‑1 (MQ=255)
                                                                                                                                       tgctTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGc                             >  5:215748/1‑136 (MQ=255)
                                                                                                                                        gctTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTa                 >  7:306912/1‑147 (MQ=255)
                                                                                                                                            ctctGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAg              <  7:60032/146‑1 (MQ=255)
                                                                                                                                               tGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACgcaatg                                                                                      >  7:219734/1‑71 (MQ=255)
                                                                                                                                               tGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACgcaatg                                                                                      <  8:219734/71‑1 (MQ=255)
                                                                                                                                                   aTGTGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCAtc     <  1:324480/148‑1 (MQ=255)
                                                                                                                                                      tgACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTcc  <  8:159155/148‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
GCCGTCAATACTTGCCAATCCTGCCGGGCCAGCAAGGCCCGCATGCTGATGAGGTTCTCGGGGTAGTCGTCTACCACCAGAAGAACTGAGCTGCCATCCAGTGGCTGCTG‑‑TTGAGGGTAAGCGCCATCCATGCTGCTTCTCTGTTATGTGACCGACTTCAAAATACGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCC  >  NC_002947/5651062‑5651357

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: