Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,951,635 C→T 7.9% V636M (GTG→ATG)  mnmC ← methyltransferase/FAD‑dependent oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,951,6350CT 7.9% 180.1 / 5.2 63V636M (GTG→ATG) mnmCmethyltransferase/FAD‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (31/27);  new base T (2/3);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

TTAGGGCACCTTAAGGAACTGACCTATGTGCTTTTCGAGAGTTTTAGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGC  >  NC_002947/1951488‑1951783
                                                                                                                                                   |                                                                                                                                                      
ttAGGGCACCTTAAGGAACTGACCTATGTGCTTTTCGAGAGTTTTAGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACt                                                                                                                                                       <  3:135183/149‑1 (MQ=255)
  aGGGCACCTTAAGGAACTGACCTATGTGCTTTTCGAGAGTTTTAGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGc                                                                                                                                                     >  7:230463/1‑149 (MQ=255)
   gggCACCTTAAGGAACTGACCTATGTGCTTTTCGAGAGTTTTAGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGcc                                                                                                                                                    <  2:27815/149‑1 (MQ=255)
    ggCACCTTAAGGAACTGACCTATGTGCTTTTCGAGAGTTTTAGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCATTGccg                                                                                                                                                   >  4:87190/1‑148 (MQ=255)
    ggCACCTTAAGGAACTGACCTATGTGCTTTTCGAGAGTTTTAGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGccc                                                                                                                                                   >  4:224444/1‑149 (MQ=255)
         cTTAAGGAACTGACCTATGTGCTTTTCGAGAGTTTTAGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCgg                                                                                                                                               >  5:261224/1‑148 (MQ=255)
                aaCTGACCTATGTGCTTTTCCAGAGTTTTAGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCACACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCa                                                                                                                                       <  3:230749/149‑1 (MQ=255)
                  cTGACCTATGTGCTTTTCGAGAGTTTTAGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCATTGCCGCGGGCAGTGGcgca                                                                                                                                     <  3:87190/149‑4 (MQ=255)
                     aCCTATGTGCTTTTCGAGAGTTTTAGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCATTGCCGCGGGCAGTGGcgcagg                                                                                                                                   <  4:215558/148‑6 (MQ=255)
                       cTATGTGCTTTTCGAGAGTTTTAGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTc                                                                                                                              >  7:281575/1‑149 (MQ=255)
                       cTATGTGCTTTTCGAGAGTTTTAGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTc                                                                                                                              <  8:230463/149‑1 (MQ=255)
                          tgtgCTTTTCGAGAGTTTTAGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTcccc                                                                                                                           <  8:182594/149‑1 (MQ=255)
                            tgCTTTTCGAGAGTTTTAGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCca                                                                                                                          >  4:321539/1‑148 (MQ=255)
                               ttttCGAGAGTTTTAGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCCCATa                                                                                                                      >  4:169190/1‑149 (MQ=255)
                               ttttCGAGAGTTTTAGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAAACCAAACCGGTTTGGATGGCAAGCCTGCGCCATTGCCGCGGGCAGTGGcgca                                                                                                                                     <  4:354439/136‑4 (MQ=255)
                               ttttCGAGAGTTTTAGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAAACCAAACCGGTTTGGATGGCAAGCCTGCGCCATTGCCGCGGGCAGTGGcgca                                                                                                                                     >  3:354439/1‑133 (MQ=255)
                                    gagaGTTTTCGAATGTTGGGTTAAATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCac                                                                                                                 <  8:105577/149‑2 (MQ=255)
                                    gagaGTTTTAGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATGAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAg                                                                                                                 >  7:175151/1‑149 (MQ=255)
                                     agagTTTTAGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAgg                                                                                                                >  1:362132/1‑149 (MQ=255)
                                      gagTTTTAGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGgc                                                                                                               <  6:384149/149‑1 (MQ=255)
                                         ttttAGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGgccgc                                                                                                            >  3:171313/1‑149 (MQ=255)
                                         ttttAGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGgccgc                                                                                                            >  5:285128/1‑149 (MQ=255)
                                          tttAGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGgccgcc                                                                                                           >  1:327415/1‑149 (MQ=255)
                                             aGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGccac                                                                                                         <  6:261224/148‑1 (MQ=255)
                                             aGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGccac                                                                                                         >  6:227796/1‑148 (MQ=255)
                                             aGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGccac                                                                                                         <  3:224444/148‑1 (MQ=255)
                                               aaTGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGccac                                                                                                         >  1:68181/1‑146 (MQ=255)
                                                       gTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGAAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCAGCCACCAGCTCGCCa                                                                                               >  6:99772/1‑148 (MQ=255)
                                                         tataTGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACt                                                                                             >  1:203246/1‑148 (MQ=255)
                                                                   caggcaggCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCg                                                                                  >  8:204082/1‑149 (MQ=255)
                                                                   caggcaggCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCg                                                                                  >  5:52890/1‑149 (MQ=255)
                                                                    aggcaggcAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGt                                                                                 <  8:140077/149‑1 (MQ=255)
                                                                      gcaggcaACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCa                                                                                               >  6:294789/1‑133 (MQ=255)
                                                                      gcaggcaACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCa                                                                                               <  5:294789/133‑1 (MQ=255)
                                                                      gcaggcaACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGa                                                                               <  3:321539/149‑1 (MQ=255)
                                                                       caggcaACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGAt                                                                              <  2:68181/149‑1 (MQ=255)
                                                                         ggcaACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGAAATGGCAGCGCGAGTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGAt                                                                              >  1:319849/1‑149 (MQ=255)
                                                                            aaCCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAAcc                                                                         <  2:327415/149‑1 (MQ=255)
                                                                                 gATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGt                                                                    <  8:175151/149‑1 (MQ=255)
                                                                                 gATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGt                                                                    <  3:69279/149‑1 (MQ=255)
                                                                                      aaTCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTTGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCGACCAGCTCGCGACTCAACGGTGCCGTGATCAACCCGCGTGAGcc                                                               <  5:99772/149‑1 (MQ=255)
                                                                                      aaTCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGcc                                                               <  2:362132/149‑1 (MQ=255)
                                                                                        tCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCAt                                                             <  3:217089/149‑1 (MQ=255)
                                                                                        tCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCAt                                                             <  8:281575/149‑1 (MQ=255)
                                                                                           tttCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGcc                                                          >  1:230184/1‑149 (MQ=255)
                                                                                             tCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGccc                                                         <  6:285128/148‑1 (MQ=255)
                                                                                                    gATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTca                                                 >  4:279842/1‑149 (MQ=255)
                                                                                                          ccGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATAc                                            >  6:27319/1‑148 (MQ=255)
                                                                                                            cttgccaTCCAAACCGGTTTGGATTGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAAc                                         >  2:284500/7‑149 (MQ=255)
                                                                                                                ccaTCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCGCGCTGTTCACATACAACCCGc                                     >  2:142246/3‑149 (MQ=255)
                                                                                                                   tCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGcc                                    >  1:146689/1‑147 (MQ=255)
                                                                                                                      aaaCCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGc                                >  7:166529/1‑148 (MQ=255)
                                                                                                                      aaaCCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGc                                >  8:115468/1‑148 (MQ=255)
                                                                                                                          cGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAagg                           >  8:139560/1‑149 (MQ=255)
                                                                                                                                ggATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAagggcagg                      >  6:21312/1‑148 (MQ=255)
                                                                                                                                ggATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAagggcagg                      <  5:21312/148‑1 (MQ=255)
                                                                                                                                 gATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAagggcaggg                     >  8:203299/1‑148 (MQ=255)
                                                                                                                                      cAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATc               <  2:230184/149‑1 (MQ=255)
                                                                                                                                       aaGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATcc              >  7:218055/1‑149 (MQ=255)
                                                                                                                                         gCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCgg            <  5:227796/149‑1 (MQ=255)
                                                                                                                                          ccTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGt           <  8:218055/149‑1 (MQ=255)
                                                                                                                                           cTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTa          <  7:139560/149‑1 (MQ=255)
                                                                                                                                              cgcCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCGCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACtt       <  1:142246/149‑1 (MQ=255)
                                                                                                                                                   cTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCcgc  <  4:156188/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
TTAGGGCACCTTAAGGAACTGACCTATGTGCTTTTCGAGAGTTTTAGAATGTTTGGTTATATGTGGCCAGGCAGGCAACCCGATGCAATCATTTCCCCCGGATCAACCGGCGCAATCCAAACCGGTTTGGATGGCAAGCCTGCGCCACTGCCCGCGGCAATGGCAGCG‑G‑TTCCCCACATACCCAGGCCGCCACCAGCTCGCCACTCAACGGTGCCGTGATCAACCCGCGTGAGCCATGCCCGCTGTTCACATACAACCCGCCCAGCCAAGGGCAGGGCACATCCGGTACTTGCCGC  >  NC_002947/1951488‑1951783

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