Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,859,571 T→A 8.4% intergenic (‑121/+35) PP_1662 ← / ← PP_1663 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,859,5710TA 8.4% 178.0 / 7.2 61intergenic (‑121/+35)PP_1662/PP_1663hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/26);  new base A (2/3);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

ATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAGTTGCACCAGCAGGTAGTCCCA  >  NC_002947/1859428‑1859713
                                                                                                                                               |                                                                                                                                              
aTGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCgg                                                                                                                                          <  2:982255/150‑1 (MQ=255)
    cTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCGGCCTCAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCTAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCCCAGTGCGGGGCt                                                                                                                                      <  2:1229330/150‑1 (MQ=255)
      gCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTc                                                                                                                                    <  2:577526/150‑1 (MQ=255)
          cAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCCCAGTGCGGGGCttttctt                                                                                                                                <  1:71805/150‑6 (MQ=255)
                 tCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTgcag                                                                                                                         >  1:534097/1‑150 (MQ=255)
                   gCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTgcaggc                                                                                                                       >  2:989431/1‑150 (MQ=255)
                     aCGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTgcaggcag                                                                                                                     >  2:1299216/1‑150 (MQ=255)
                        tcttctTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTgc                                                                                                                  >  2:1010255/1‑150 (MQ=255)
                             tttATCGATTGTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGAt                                                                                                             <  1:733673/150‑1 (MQ=255)
                             tttATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGAt                                                                                                             <  2:241796/150‑1 (MQ=255)
                               tataGAGTTTGTTAGCAGACTAACTAAATAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGGAAAGACGAGCTGCTCACTAGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCa                                                                                                           <  1:602986/150‑1 (MQ=255)
                                   gagTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTg                                                                                                       <  1:1060920/150‑1 (MQ=255)
                                    agTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGc                                                                                                      <  2:383835/150‑1 (MQ=255)
                                      tttGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGTGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCa                                                                                                           >  1:385232/1‑143 (MQ=255)
                                      tttGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGTGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCa                                                                                                           <  2:385232/143‑1 (MQ=255)
                                      tttGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTgcagg                                                                                                                        >  1:904995/1‑130 (MQ=255)
                                      tttGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTgcagg                                                                                                                        <  2:904995/130‑1 (MQ=255)
                                               aGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGc                                                                                           <  2:774729/150‑1 (MQ=255)
                                                   taagtaagTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCCCAGTGCGGGGCtttt                                                                                                                                   <  1:1074658/106‑3 (MQ=255)
                                                   taagtaagTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCCCAGTGCGGGGCtttt                                                                                                                                   >  2:1074658/1‑104 (MQ=255)
                                                       taagtaaAACATGTATGAAACAATTTTCATGTAGCTCGAAGCGATGCTGCCTGGGTGGGCCTTTTCGCAGCAAAGAAAAGCCCCGTACTCGGCGGTGCTTCAGGTATTGCAGGCAGTGCGCCATCTGTTGGCCTGCAATGGCTTGCCCTc                                                                                   <  2:152444/150‑1 (MQ=255)
                                                         agtaaAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTGGCGATAGa                                                                                 >  1:60484/1‑150 (MQ=255)
                                                               aCAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGc                                                                           >  1:478531/1‑150 (MQ=255)
                                                               aCAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGc                                                                           <  1:989431/150‑1 (MQ=255)
                                                                    cATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGGTGCCCTTCACCGGCTTGCCATCGACCGTGCCATcc                                                                      >  2:760815/1‑150 (MQ=255)
                                                                     aTGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCt                                                                     >  2:363338/1‑150 (MQ=255)
                                                                           cAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGtt                                                                                                        >  2:1113454/1‑109 (MQ=255)
                                                                           cAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGtt                                                                                                        <  1:1113454/109‑1 (MQ=255)
                                                                           cAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTgcagca                                                               >  1:911857/1‑150 (MQ=255)
                                                                            aaGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCAc                                                              >  2:1033705/1‑150 (MQ=255)
                                                                                   aTGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGACACGAACAGCCGCAGGCATGAAGCGAGGTATTGCAGGCAGGGCGAGATCAGTTGCCCTTCACCGCAGCGCCAGCGACCGGGCCATCCTGCAGCACGATcaca                                                       <  1:1033705/150‑1 (MQ=255)
                                                                                   aTGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATcaca                                                       >  2:926860/1‑150 (MQ=255)
                                                                                   aTGAAGCTCGAAGGCACGATTGCGGTGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATcaca                                                       >  2:59705/1‑150 (MQ=255)
                                                                                     gAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATa                                                     >  2:133731/1‑150 (MQ=255)
                                                                                        gCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCa                                                                                                <  1:168774/104‑1 (MQ=255)
                                                                                        gCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCa                                                                                                >  2:168774/1‑104 (MQ=255)
                                                                                          tCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCTTTCGTGCCATCCTGCAGCACGATCACATACTCCt                                                <  1:1045469/150‑1 (MQ=255)
                                                                                            gAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTg                                                                            >  1:1060468/1‑120 (MQ=255)
                                                                                            gAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTg                                                                            <  2:1060468/120‑1 (MQ=255)
                                                                                              aGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATAc                                                    >  2:283422/1‑142 (MQ=255)
                                                                                              aGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATAc                                                    <  1:283422/142‑1 (MQ=255)
                                                                                                  aTGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGcctgtg                                        <  1:760815/150‑5 (MQ=255)
                                                                                                                     tttCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAGtt                     <  1:641317/150‑1 (MQ=255)
                                                                                                                      ttCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGAcc                                                                               <  1:464280/91‑1 (MQ=255)
                                                                                                                      ttCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGAcc                                                                               >  2:464280/1‑91 (MQ=255)
                                                                                                                          cagcaAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAGTTGCAcc                >  1:96047/1‑150 (MQ=255)
                                                                                                                          cagcaAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAGTTGCAcc                >  1:96109/1‑150 (MQ=255)
                                                                                                                             caAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAGTTGCACcagc             >  2:639739/1‑150 (MQ=255)
                                                                                                                             caAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTAGCACAGTTACACcagc             >  2:397005/1‑150 (MQ=255)
                                                                                                                               aaGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAGTTGCACcagcag           <  2:534097/150‑1 (MQ=255)
                                                                                                                               aaGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAGTTGCACcagcag           >  1:661149/1‑150 (MQ=255)
                                                                                                                                  aaaaGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTATACCTgg                            >  2:744460/1‑130 (MQ=255)
                                                                                                                                  aaaaGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTATACCTgg                            <  1:744460/130‑1 (MQ=255)
                                                                                                                                     aGCCCCGCCCAGTGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAGTTGCACCAGCAGGTAGTc     >  2:715978/1‑150 (MQ=255)
                                                                                                                                      gCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAGTTGCACCAGCAGGTAGTcc    <  2:478531/150‑1 (MQ=255)
                                                                                                                                       ccccGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAGTTGCACCAGCAGGTAGTccc   <  1:1010255/150‑1 (MQ=255)
                                                                                                                                        cccGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTgca                                                                  <  2:948919/86‑1 (MQ=255)
                                                                                                                                        cccGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTgca                                                                  >  1:948919/1‑86 (MQ=255)
                                                                                                                                        cccGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAGTTGCACCAGCAGGTAGTCCCa  >  1:504632/1‑150 (MQ=255)
                                                                                                                                             aCTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATc                                                                       >  2:477755/1‑76 (MQ=255)
                                                                                                                                             aCTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATc                                                                       <  1:477755/76‑1 (MQ=255)
                                                                                                                                               tGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATa                                                     <  1:1275991/92‑1 (MQ=255)
                                                                                                                                               tGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATa                                                     >  2:1275991/1‑92 (MQ=255)
                                                                                                                                               |                                                                                                                                              
ATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAGTTGCACCAGCAGGTAGTCCCA  >  NC_002947/1859428‑1859713

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