Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 127,121 C→G 6.0% intergenic (‑31/+39) thrB ← / ← PP_0122 homoserine kinase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947127,1210CG 6.0% 200.8 / 2.1 62intergenic (‑31/+39)thrB/PP_0122homoserine kinase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/28);  new base G (2/2);  total (32/30)
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.35e-01

AGAAATTGCTGTTCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCTTGCCGTTTT  >  NC_002947/126973‑127269
                                                                                                                                                    |                                                                                                                                                    
agaaATTGCTGTTCTCGGTAACGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCAGGGTCACGGGGGTGAAGACTTACATGATGAAAAATTGCCCATAAGGGCAATTccgc                                                                                                                                                      >  5:181247/1‑149 (MQ=255)
    aTTGCTGTTCTCGGTACAGTCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGTGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCggg                                                                                                                                                  <  6:181247/149‑1 (MQ=255)
             ctcGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCgg                                                                                                                                                   >  6:93887/1‑139 (MQ=255)
             ctcGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCgg                                                                                                                                                   <  5:93887/139‑1 (MQ=255)
                ggTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCggg                                                                                                                                        <  7:41943/147‑1 (MQ=255)
                ggTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGtt                                                                                                                                      <  6:376490/149‑1 (MQ=255)
                ggTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGtt                                                                                                                                      <  6:370201/149‑1 (MQ=255)
                 gTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGt                                                                                                                                       >  4:261634/1‑147 (MQ=255)
                 gTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGt                                                                                                                                       <  3:261634/147‑1 (MQ=255)
                  tACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGa                                                                                                                                    >  2:598/1‑149 (MQ=255)
                    ccGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGAt                                                                                                                                   <  3:171797/148‑1 (MQ=255)
                       gcggcAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGGAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGa                                                                                                                               <  2:393845/149‑1 (MQ=255)
                       gcggcAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGAtgtg                                                                                                                                >  2:230374/1‑148 (MQ=255)
                            aaTGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGGCGGGAAGTGCCGGGTTGATGTGAACGGt                                                                                                                          <  1:230374/149‑1 (MQ=255)
                              tGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTgc                                                                                                                        <  7:19144/149‑1 (MQ=255)
                                  cTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGt                                                                                                                    >  1:54782/1‑149 (MQ=255)
                                     gAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACgg                                                                                                                           <  2:364023/139‑1 (MQ=255)
                                     gAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACgg                                                                                                                           >  1:364023/1‑139 (MQ=255)
                                         tCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTa                                                                                                             <  4:130481/149‑1 (MQ=255)
                                           gAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTAc                                                                                                            >  5:272543/1‑148 (MQ=255)
                                              cAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACt                                                                                                        >  3:10018/1‑149 (MQ=255)
                                                    gCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCg                                                                                                      >  1:367560/1‑145 (MQ=255)
                                                    gCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGa                                                                                                  <  2:54782/149‑1 (MQ=255)
                                                       cAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCCGGCGGAAGTGc                                                                                                                                            <  5:31955/104‑1 (MQ=255)
                                                       cAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCCGGCGGAAGTGc                                                                                                                                            >  6:31955/1‑104 (MQ=255)
                                                          cTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTcc                                                                                            <  1:27178/149‑1 (MQ=255)
                                                          cTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTc                                                                                             >  5:54063/1‑148 (MQ=255)
                                                            cGTACGGCGCCAGAAAGGTTTACAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAg                                                                                          >  7:180185/1‑149 (MQ=255)
                                                               aCGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGAc                                                                                       >  3:402288/1‑149 (MQ=255)
                                                               aCGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGAc                                                                                       >  1:42687/1‑149 (MQ=255)
                                                                   cgcCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGaa                                                                                   >  6:207058/1‑149 (MQ=255)
                                                                      cAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACAACTCGAAGATCTTCCAGGACGGAAt                                                                                  >  4:133292/1‑147 (MQ=255)
                                                                       aGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCAGCCGGGAAGTGCAGGGTTGATGTGAACGGTGAGCGTCAGCTTACCCCTCGAAGATCTTCCAGTACGGAATCAg                                                                               >  7:142083/1‑149 (MQ=255)
                                                                             gTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGc                                                                              <  8:2379/144‑1 (MQ=255)
                                                                             gTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGc                                                                              >  7:2379/1‑144 (MQ=255)
                                                                              tttCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTcc                                                                        >  7:123570/1‑149 (MQ=255)
                                                                                tCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACTACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCg                                                                       <  3:133292/148‑1 (MQ=255)
                                                                                tCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGCCCgg                                                                      >  4:30511/1‑149 (MQ=255)
                                                                                   aGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTg                                                                   <  2:90499/149‑1 (MQ=255)
                                                                                          gCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGAt                                                            <  5:207058/149‑1 (MQ=255)
                                                                                          gCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGAt                                                            >  5:116410/1‑149 (MQ=255)
                                                                                          gCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATAGCCCAAACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGAt                                                            <  4:290424/149‑1 (MQ=255)
                                                                                           ccGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATc                                                           <  4:10018/149‑1 (MQ=255)
                                                                                                tCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACg                                                                                      >  2:173071/1‑117 (MQ=255)
                                                                                                tCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACg                                                                                      <  1:173071/117‑1 (MQ=255)
                                                                                                    gggggTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAAt                                                  >  5:179726/1‑149 (MQ=255)
                                                                                                       ggTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTgg                                               <  4:402288/149‑1 (MQ=255)
                                                                                                            aGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACg                                                                                      >  4:205969/1‑105 (MQ=255)
                                                                                                            aGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACg                                                                                      <  3:205969/105‑1 (MQ=255)
                                                                                                            aGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCa                                          <  6:272543/149‑1 (MQ=255)
                                                                                                               cTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTg                                       <  1:598/149‑1 (MQ=255)
                                                                                                                     tgatgaAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCa                                 <  8:123570/149‑1 (MQ=255)
                                                                                                                         gaAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGATTCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAgcgc                             >  7:76972/1‑149 (MQ=255)
                                                                                                                         gaAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAgcgc                             >  6:425569/1‑149 (MQ=255)
                                                                                                                          aaaaaTTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGcgcg                            >  6:143415/1‑149 (MQ=255)
                                                                                                                          aaaaaTTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGcgcg                            >  2:310931/1‑149 (MQ=255)
                                                                                                                          aaaaaTTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGcgcg                            >  6:257685/1‑149 (MQ=255)
                                                                                                                                   ccATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGaa                    >  8:351693/1‑148 (MQ=255)
                                                                                                                                     aTACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGaaaa                  >  7:359909/1‑148 (MQ=255)
                                                                                                                                      tACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGCCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATa                <  3:30511/149‑1 (MQ=255)
                                                                                                                                       cccggCACTTCCCGGCGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAAt                 >  6:355911/4‑147 (MQ=255)
                                                                                                                                         gggCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAgg             <  6:54063/149‑1 (MQ=255)
                                                                                                                                             aCTTCCGCCGCGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTAGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCAATCAGCGCGTACCAGAAAATAAGGCTTg         <  8:142083/149‑1 (MQ=255)
                                                                                                                                                   gccgGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCtt          >  2:20177/1‑142 (MQ=255)
                                                                                                                                                   gccgGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCtt          <  1:20177/142‑1 (MQ=255)
                                                                                                                                                    ccgGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCTTGCCGtttt  >  7:151170/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGAAATTGCTGTTCTCGGTACCGGCGGCAATGCCCTGGAAGTCGAGCAGACGGCCCAGCTCGTACGGCGCCAGAAAGGTTTCCAGCTCAGGCCGGGTCACGGGGGTGAAGACTGACATGATGAAAAATTGCCCATACGGGCACTTCCGCCGGGAAGTGCCGGGTTGATGTGAACGGTGCGCGTCAGCTTACCACTCGAAGATCTTCCAGGACGGAATCAGCATGTCCGGCTGGTCGGATCGAATGAAATTGGCATCAGTGCCATCAGCGCGTACCAGAAAATAAGGCTTGCCGTTTT  >  NC_002947/126973‑127269

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