Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,998,172 A→C 6.4% intergenic (‑72/‑144) PP_4405 ← / → PP_4406 sensory box protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,998,1720AC 6.4% 195.3 / 4.4 63intergenic (‑72/‑144)PP_4405/PP_4406sensory box protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (33/26);  new base C (2/2);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

ATCGAACAAGGTGTGCAATTGCTGTTCACTGAAGCCGGCAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGCCAACGAGATTTCCAAGGAGGAAACAT  >  NC_002947/4998035‑4998317
                                                                                                                                          |                                                                                                                                                 
aTCGAACAAGGTGTGCAATTGCTGTTCACTGAAGCCGGCAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAgat                                                                                                                                         >  3:42520/1‑148 (MQ=255)
 tCGAACAAGGTGTGCAATTGCTGTTCACTGAAGCCGGCAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCTCGCGGCaatcaac                                                                                                                                       <  4:184067/149‑6 (MQ=255)
   gAACAAGGTGTGCAATTGCTGTTCACTGAAGCCGGCAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAgatgat                                                                                                                                      <  5:183174/148‑1 (MQ=255)
    aaCAAGGTGTGCAATTGCTGTTCACTGAAGCCGGCAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGc                                                                                                                                    <  4:42520/149‑1 (MQ=255)
             tgCAATTGCTGTTCACTGAAGCCGGCAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACgg                                                                                                                           >  1:240768/1‑149 (MQ=255)
               cAATTGCTGTTCACTGAAGCCGGCAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACggg                                                                                                                          <  7:306773/148‑1 (MQ=255)
                    gCTGTTCACTGAAGCCGGCAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGc                                                                                                                     <  2:240768/148‑1 (MQ=255)
                        ttCACTGAAGCCGGCAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATg                                                                                                                >  4:58311/1‑149 (MQ=255)
                             tGAAGCCGGCAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCAGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATTATGCCTAGCACGGGCAAGCATAAGcccac                                                                                                           >  8:293210/1‑149 (MQ=255)
                                  ccGGCAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACg                                                                                                      >  8:103637/1‑149 (MQ=255)
                                      cAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCccacccac                                                                                                       >  1:299821/1‑144 (MQ=255)
                                      cAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCccacccac                                                                                                       <  2:299821/144‑1 (MQ=255)
                                             cGGCCTGGTGTCTTGACTCATGGACGCTCTCTCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGa                                                                                             >  1:35339/1‑148 (MQ=255)
                                                ccTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAg                                                                                        >  5:63513/1‑149 (MQ=255)
                                                 cTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGTAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGt                                                                                       <  7:28448/149‑1 (MQ=255)
                                                 cTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGc                                                                                                                     >  7:259991/1‑119 (MQ=255)
                                                 cTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGc                                                                                                                     <  8:259991/119‑1 (MQ=255)
                                                 cTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGt                                                                                       <  7:293210/149‑1 (MQ=255)
                                                    gtgtCTTGTCTCTTGGAGGCTGTCCGTGATGACCATCAACTTCGCTCAGATCAAGTGGCTGACGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCTAGGATATGCCCACCCACGTGCACCCGAAATAGTGc                                                                                     <  5:157428/148‑1 (MQ=255)
                                                    gtgtCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGcc                                                                                    <  2:35339/149‑1 (MQ=255)
                                                      gtCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCaa                                                                                  >  7:224127/1‑149 (MQ=255)
                                                         ttGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATg                                                                                                                >  6:174393/1‑116 (MQ=255)
                                                         ttGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATg                                                                                                                <  5:174393/116‑1 (MQ=255)
                                                         ttGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGaaa                                                                                           <  6:316633/137‑1 (MQ=255)
                                                         ttGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGc                                                                               >  6:336283/1‑149 (MQ=255)
                                                         ttGACTCATGGACGCTCTCACTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGaaa                                                                                           >  5:316633/1‑137 (MQ=255)
                                                            aCTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCt                                                                            >  3:219976/1‑149 (MQ=255)
                                                                aTGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGAtt                                                                        >  1:175063/1‑149 (MQ=255)
                                                                   gACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGTGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTc                                                                      <  5:294814/148‑1 (MQ=255)
                                                                             ccTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGcaca                                                            >  4:69139/1‑149 (MQ=255)
                                                                             ccTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGca                                                              <  5:336283/147‑1 (MQ=255)
                                                                             ccTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGca                                                              >  2:267712/1‑147 (MQ=255)
                                                                                gATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAgcg                                                         >  2:118507/1‑149 (MQ=255)
                                                                                gATGCCCAGCAACTGCGCTCCGATCAAGGCGCTGCAGACAGTGGTCTTGATTGCCGCGAGAGCCAGGATGATGCCTATCACGGGCAATCATATGCCCACCCACGTGCAACCGaaa                                                                                           >  3:299140/1‑115 (MQ=255)
                                                                                  tGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCt                                                       >  5:167090/1‑149 (MQ=255)
                                                                                     ccAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTgg                                                     >  8:155204/1‑148 (MQ=255)
                                                                                     ccAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGc                                                    <  4:219976/149‑1 (MQ=255)
                                                                                      cAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCa                                                   >  6:115929/1‑149 (MQ=255)
                                                                                                 cTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATg                                                 >  2:199385/1‑140 (MQ=255)
                                                                                                 cTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATg                                                 <  1:199385/140‑1 (MQ=255)
                                                                                                  tCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTgcgc                                       >  1:638/1‑149 (MQ=255)
                                                                                                      aTCAAGGGGCTGCCGACAGTTGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTgg                                   <  6:63513/149‑1 (MQ=255)
                                                                                                        cAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTa                                 >  1:179118/1‑149 (MQ=255)
                                                                                                          aGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCTCAGCGATGACATGGATTTTGCGCTTGGTAcc                               >  8:126971/1‑149 (MQ=255)
                                                                                                           ggggCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCt                              <  7:155204/149‑1 (MQ=255)
                                                                                                                   cGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCt                                      >  2:298174/1‑133 (MQ=255)
                                                                                                                   cGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCt                                      <  1:298174/133‑1 (MQ=255)
                                                                                                                   cGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGCCAACga                      <  1:118507/149‑1 (MQ=255)
                                                                                                                      cAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGCCAACGAGAt                   >  1:6139/1‑149 (MQ=255)
                                                                                                                       aGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCATGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTg                                    >  6:40581/1‑131 (MQ=255)
                                                                                                                       aGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGCCAACGAGAtt                  >  7:44295/1‑149 (MQ=255)
                                                                                                                         tggtgTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGCCAACGAGATTTc                >  7:254945/1‑149 (MQ=255)
                                                                                                                           gtgtTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGCCAACGAGATTTcc               >  3:3470/1‑148 (MQ=255)
                                                                                                                           gtgtTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGCCAACGAGATTTcc               <  6:246389/148‑1 (MQ=255)
                                                                                                                              ttGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGCCAACGAGATTTCCAag            <  8:224127/148‑1 (MQ=255)
                                                                                                                                tcttGCCGCTCGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGCCAACGAGATTTCCAaggag         >  1:246396/3‑149 (MQ=255)
                                                                                                                                gATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGCCAACGAGATTTCCAagga          <  8:254945/148‑1 (MQ=255)
                                                                                                                                  ttGCCGCTCGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTg                                          >  7:176360/1‑114 (MQ=255)
                                                                                                                                  ttGCCGCTCGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTg                                          <  8:176360/114‑1 (MQ=255)
                                                                                                                                    gCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGCCAACGAGATTTCCAAGGAGGaa      <  7:103637/148‑1 (MQ=255)
                                                                                                                                     ccGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGCCAACGAGATTTCCAAGGAGGaaa     >  6:250334/1‑148 (MQ=255)
                                                                                                                                     ccGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGCCAACGAGATTTCCAAGGAGGaaa     <  3:58311/148‑1 (MQ=255)
                                                                                                                                       gcgAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGCCAACGAGATTTCCAAGGAGGAAACAt  <  4:3470/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                 
ATCGAACAAGGTGTGCAATTGCTGTTCACTGAAGCCGGCAGGTACCGGCCTGGTGTCTTGACTCATGGACGCTCTCCCTGATGCCCAGCAACTGCGCTCCGATCAAGGGGCTGCCGACAGTGGTGTTGATTGCCGCGAGCGGCAAGATGATGCCTATCACGGGCAAGCATATGCCCACCCACGTGCAACCGAAATAGTGCCAACGCCCTGATTTCGACCGGCACAGCGCTGGCATGGATTTTGCGCTTGGTACCTGCCAACGAGATTTCCAAGGAGGAAACAT  >  NC_002947/4998035‑4998317

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: