Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,741,935 C→G 7.6% intergenic (‑23/+61) PP_1538 ← / ← PP_1539 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,741,9350CG 7.6% 174.7 / 4.4 62intergenic (‑23/+61)PP_1538/PP_1539hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (26/30);  new base G (2/3);  total (29/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 = 3.17e-01

CACCTCGGTGACTGAAGAGGGAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTACCTGGTGCGGTCCGTCGAA  >  NC_002947/1741787‑1742070
                                                                                                                                                    |                                                                                                                                       
cacCTCGGTGACTGAAGAGGGAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATc                                                                                                                                         >  6:75426/1‑149 (MQ=255)
 acCTCGGTGACTGAAGAGGGAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATcg                                                                                                                                        >  7:200420/1‑149 (MQ=255)
        tGACTGAAGAGGGAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCGATATGGCGGCaaa                                                                                                                                 <  1:189962/149‑3 (MQ=255)
        tGACTGAAGAGGGAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGt                                                                                                                                 <  4:147636/149‑1 (MQ=255)
         gACTGAAGAGGGAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGt                                                                                                                                 >  4:155208/1‑148 (MQ=255)
         gACTGAAGAGGGAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGt                                                                                                                                 <  6:131370/148‑1 (MQ=255)
         gACTGAAGAGGGAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGt                                                                                                                                 >  5:131370/1‑148 (MQ=255)
          aCTGAAGAGGGAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCGGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCGATATGGCGGCaaag                                                                                                                                <  5:414617/148‑4 (MQ=255)
           gtGAAGAGGGAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTCCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATTGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGaa                                                                                                                              <  7:350650/148‑1 (MQ=255)
            tGAAGAGGGAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATcgc                                                                                                                                       <  3:49255/139‑1 (MQ=255)
            tGAAGAGGGAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATcgc                                                                                                                                       >  4:49255/1‑139 (MQ=255)
            tGAAGAGGGAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGaa                                                                                                                              >  3:339272/1‑148 (MQ=255)
            tGAAGAGGGAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGaa                                                                                                                              >  3:112707/1‑148 (MQ=255)
                  gggAAAGCCTGGCTTAACCGGGCTTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGt                                                                                                                        <  4:339272/148‑1 (MQ=255)
                  gggAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTg                                                                                                                       <  3:357080/149‑1 (MQ=255)
                  gggAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTg                                                                                                                       <  5:302922/149‑1 (MQ=255)
                    gAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGga                                                                                                                     <  5:75426/149‑1 (MQ=255)
                      aaGCCTGGCTTAACCGGGCGTCGGGTGACTGGCTTCCGTTCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTATGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGgcga                                                                                                                   <  6:325398/149‑1 (MQ=255)
                       aGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGagag                                                                                                                  >  1:235011/1‑149 (MQ=255)
                       aGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGagag                                                                                                                  <  4:112707/149‑1 (MQ=255)
                          cTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAggg                                                                                                                <  1:238575/148‑1 (MQ=255)
                          cTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAggg                                                                                                                >  2:238575/1‑148 (MQ=255)
                           tGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGt                                                                                                              >  1:132699/1‑149 (MQ=255)
                                    cGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTTCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTTAATAGAGTGGAGAGGGGTTACAGCggg                                                                                                     <  7:239533/149‑1 (MQ=255)
                                         gTCGGGTGACTGGCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCggggtggg                                                                                                <  7:113239/149‑1 (MQ=255)
                                          tCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCggggtgggg                                                                                               <  4:111850/149‑1 (MQ=255)
                                             ggTGACTGTCTTCCGTCCTTCCAGGAAGGCTCGGACCAACGGCCCGCTGAACAGGATCCGGCCTTCCTTTGCTAGAGACATAGCTTCGCCTTTGCCGCCATATAGCGGCAGTCAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTa                                                                                             >  2:36135/1‑148 (MQ=255)
                                              gTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTa                                                                                             >  6:175898/1‑147 (MQ=255)
                                                     tCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAAGAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGc                                                                                    >  2:22812/1‑149 (MQ=255)
                                                       ttCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGctc                                                                                  >  3:20083/1‑149 (MQ=255)
                                                       ttCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGctc                                                                                  >  3:14603/1‑149 (MQ=255)
                                                          cGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCg                                                                                                       <  1:333216/125‑1 (MQ=255)
                                                          cGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCg                                                                                                       >  2:333216/1‑125 (MQ=255)
                                                           gTCCTTCGAGGTTGGCTCGGACCATCGGCCCGCTGAACAGGTTCGGGGGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGGAGTGAATAGAGTGTTGAGGGGTTTCAGCGGGGTGGGGTACAGATGTGCTCTggg                                                                              <  1:36135/149‑1 (MQ=255)
                                                            tCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTgggg                                                                             <  6:312284/149‑1 (MQ=255)
                                                                  cAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTaa                                                                         >  2:271152/1‑147 (MQ=255)
                                                                    ggATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGc                                                                     >  6:25784/1‑149 (MQ=255)
                                                                      aTGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTg                                                                   <  5:25784/149‑1 (MQ=255)
                                                                       tGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTgg                                                                  >  7:386618/1‑149 (MQ=255)
                                                                       tGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTgg                                                                  >  5:272174/1‑149 (MQ=255)
                                                                              gACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCggg                                                           >  2:296100/1‑149 (MQ=255)
                                                                                     ggCCCGCTGAACAGGATCGGTCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTgcg                                                    <  6:5381/149‑1 (MQ=255)
                                                                                          gCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGc                                               <  5:175898/149‑1 (MQ=255)
                                                                                              aaCAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTg                                           >  5:4317/1‑149 (MQ=255)
                                                                                              aaCAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTg                                           >  3:72223/1‑149 (MQ=255)
                                                                                                 aGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGGGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCa                                        <  3:155208/149‑1 (MQ=255)
                                                                                                   gATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGATTAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGc                                               >  2:379830/1‑140 (MQ=255)
                                                                                                   gATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGGGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGc                                               <  1:379830/140‑1 (MQ=255)
                                                                                                        ggCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATgg                                 >  6:85235/1‑149 (MQ=255)
                                                                                                              ccTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGTTGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCgg                            <  8:386618/148‑1 (MQ=255)
                                                                                                               cgttGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCgg                            >  8:47837/3‑147 (MQ=255)
                                                                                                               cgttGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGGGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCgg                            <  7:47837/145‑1 (MQ=255)
                                                                                                               cTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCa                          >  2:353113/1‑149 (MQ=255)
                                                                                                                       gagaCATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTAcc                   <  1:296100/148‑1 (MQ=255)
                                                                                                                        agaCATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGTGGGGGTACTGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTACCTg                 <  2:132699/149‑1 (MQ=255)
                                                                                                                           cATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTACCTGGTg              <  2:235011/149‑1 (MQ=255)
                                                                                                                            aTAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTACCTGGTGc             >  7:120426/1‑149 (MQ=255)
                                                                                                                            aTAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTACCTGGTGc             >  2:70099/1‑149 (MQ=255)
                                                                                                                                 ttCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCggg                                                                                                     >  6:65458/1‑56 (MQ=255)
                                                                                                                                   cGCCTTTGCCGCCATATCGCGGCAGTGAATCGAGGGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTACCTGGTGCGGTCCg       <  8:200420/148‑1 (MQ=255)
                                                                                                                                    gCCTTTGCCGCCATATCGCGGCAGTTAATAGCGTGGTGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTACCTGGTGCGGTCCGTc     <  4:14603/149‑1 (MQ=255)
                                                                                                                                    gCCTTTGCCGCCATATCGCGGCAGTGAACAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCGGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTACCTGGTGCGGTCCGTc     <  1:22812/149‑1 (MQ=255)
                                                                                                                                    gCCTTGGGCGCCTTATGGCGGCAGTTAATAGAGTTGAGAGGGGTTACAGCGGGGTTGGGTACAGATGTGCTGTGGGGTTATGAGCTGGCACCGGTGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTACCTGGTGCGGTCCGTc     <  4:20083/149‑1 (MQ=255)
                                                                                                                                       actGCCGCGATATGGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCTCCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTACCTGGTGCGGTCCGTCGaa  >  6:1418/3‑149 (MQ=255)
                                                                                                                                       actGCCGCGATATGGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTACCTGGTGCGGTCCGTCGaa  >  3:286062/3‑149 (MQ=255)
                                                                                                                                         tGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTACCTGGTGCGGTCCGTCGaa  <  6:272174/147‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                       
CACCTCGGTGACTGAAGAGGGAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTACCTGGTGCGGTCCGTCGAA  >  NC_002947/1741787‑1742070

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑