Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,577,416 G→A 6.9% D59D (GAC→GAT)  gntK ← gluconate kinase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,577,4160GA 6.9% 182.5 / 2.4 58D59D (GAC→GAT) gntKgluconate kinase 2
Reads supporting (aligned to +/- strand):  ref base G (26/28);  new base A (2/2);  total (28/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.99e-01

GAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACG  >  NC_000913/3577272‑3577560
                                                                                                                                                |                                                                                                                                                
gAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCAt                                                                                                                                              >  6:311467/1‑149 (MQ=255)
      gagaTTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTg                                                                                                                                        <  1:159292/149‑1 (MQ=255)
             ggATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCATCGTCATTCAGTGGTTCGcc                                                                                                                                 <  3:277942/149‑1 (MQ=255)
             ggATTACCTTCACGCAGCAAGTCGCGATAGTGATTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGc                                                                                                                                  >  2:175887/1‑148 (MQ=255)
                ttACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGa                                                                                                                              >  4:315361/1‑149 (MQ=255)
                ttACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGa                                                                                                                              <  6:131953/149‑1 (MQ=255)
                   ccTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCAtt                                                                                                                                             <  8:18915/131‑1 (MQ=255)
                   ccTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCAtt                                                                                                                                             >  7:18915/1‑131 (MQ=255)
                    cTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGc                                                                                                                           >  8:51404/1‑148 (MQ=255)
                      tCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGTCAt                                                                                                                        >  2:295372/1‑149 (MQ=255)
                      tCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCAt                                                                                                                        >  6:898/1‑149 (MQ=255)
                              caAGTAGCGATAGTGTTTTTTCAATGCAGAACACACGATCAGCGACACTTTATTAGTGCGCTGCATATCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGa                                                                                                                >  4:204830/1‑149 (MQ=255)
                                      gATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGt                                                                                                                                         <  6:39793/116‑1 (MQ=255)
                                      gATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGt                                                                                                                                         >  5:39793/1‑116 (MQ=255)
                                      gATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTgcgc                                                                                                        >  2:84317/1‑149 (MQ=255)
                                          gtgtTTTTTCAATGCAGCACAGACGATCCGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCAAGACGCCATTTTTTCGATATTGCGCCGtg                                                                                                    <  1:124223/149‑1 (MQ=255)
                                          gtgtTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGt                                                                                                     <  8:289452/148‑1 (MQ=255)
                                             ttttttCAATGCAGAACATACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGtgga                                                                                                  <  7:121859/148‑1 (MQ=255)
                                             ttttttCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGtggat                                                                                                 >  3:232134/1‑149 (MQ=255)
                                             ttttttCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGtgga                                                                                                  <  8:226787/148‑1 (MQ=255)
                                                  tCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATggag                                                                                             <  3:206379/148‑1 (MQ=255)
                                                        cAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATc                                                                                       <  1:39856/148‑1 (MQ=255)
                                                          gAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGcc                                                                                    <  1:91145/149‑1 (MQ=255)
                                                          gAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGc                                                                                     >  4:114364/1‑148 (MQ=255)
                                                           aaCAGACGATCAGCGACACTTTATTAGTGCGCTGCTTAGCATAGGAGGCGTCGTTCAGCGCCTGCAAGCACGGTTTGCGATCGTCGTAATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCa                                                                                   <  8:254693/149‑1 (MQ=255)
                                                             cAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATc                                                                                 <  4:199125/149‑1 (MQ=255)
                                                               gACGATCAGCGACACTTTATTAGTGCGCTGCATGGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCaa                                                                               <  7:335043/149‑1 (MQ=255)
                                                                  gATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGa                                                                             >  5:91800/1‑148 (MQ=255)
                                                                      aGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGTTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGTCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAAcgc                                                                        <  1:295372/149‑1 (MQ=255)
                                                                       gCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCATCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCa                                                                                <  2:160905/140‑1 (MQ=255)
                                                                       gCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCATCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCa                                                                                >  1:160905/1‑140 (MQ=255)
                                                                             cTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGa                                                                 <  3:156256/149‑1 (MQ=255)
                                                                              tttATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATg                                                                  >  1:168524/1‑147 (MQ=255)
                                                                              tttATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATg                                                                  >  1:316685/1‑147 (MQ=255)
                                                                                     gTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGAt                                                          <  5:311467/148‑1 (MQ=255)
                                                                                      tGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATgc                                                        <  1:84317/149‑1 (MQ=255)
                                                                                             cATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCAAtt                                                 >  3:18096/1‑149 (MQ=255)
                                                                                                  cAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTg                                            <  8:107671/149‑1 (MQ=255)
                                                                                                  cAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCATCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTg                                            >  3:98886/1‑149 (MQ=255)
                                                                                                   aaaCGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTgg                                           >  7:31633/1‑149 (MQ=255)
                                                                                                       gcggcgTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTGATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGGCGTGGATGGAGGAAATCGCCATCTAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGAc                                       >  5:280304/1‑149 (MQ=255)
                                                                                                         ggcgTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGCAACGCGGCATGCAGTTGCTGCGCCAGTTCACTGGCGACCg                                     >  6:110005/1‑149 (MQ=255)
                                                                                                            gtcgtTCAGCGCCTGCAACCACGGTTTGCGATCGTCATCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAg                                  <  4:98886/149‑1 (MQ=255)
                                                                                                             tcgtTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTgg                                           >  6:276890/1‑139 (MQ=255)
                                                                                                             tcgtTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTgg                                           <  5:276890/139‑1 (MQ=255)
                                                                                                                ttCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCCCTTCACTGGCGACCGCAGAttt                              >  1:248057/1‑149 (MQ=255)
                                                                                                                   aGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCCTGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGc                            <  7:51404/148‑1 (MQ=255)
                                                                                                                   aGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGc                            <  4:313893/148‑1 (MQ=255)
                                                                                                                      gcCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCt                        <  8:31633/149‑1 (MQ=255)
                                                                                                                         tGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGATGATGCGCCACTTCACTGTCGACCGCAGATTTGCCGCTGcc                     >  2:296568/1‑149 (MQ=255)
                                                                                                                          gCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTAGACAGACGCCATTTTTTCGATATTGCGCCGTGGATGGTGGAAATCGCCATCAAGAAACGCGGCATGATGTTGATGAGAGACTTCa                                               >  6:80061/1‑122 (MQ=255)
                                                                                                                               cACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATAc               >  5:31407/1‑149 (MQ=255)
                                                                                                                                aCGCTTTGCGATCGTCGTCATTCAGTGGTTCCCCACACCCCATTATTTAGATATTGCTCCGTGGATGGAGGAAATCGCCAACAAGAAACGCGGCATGAACTTGATGCGCCACATCAATGGCGACCGCAGATTTGCCGCTGACCGATACg              >  8:83474/1‑149 (MQ=255)
                                                                                                                                 cGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGc             <  5:898/149‑1 (MQ=255)
                                                                                                                                    tttGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCa          >  2:41138/1‑149 (MQ=255)
                                                                                                                                       gCGATCGGCGGCATTCAGTGGATCGCCAGACGCCATTTTTTCGATATAGCTCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATc        <  1:175887/148‑1 (MQ=255)
                                                                                                                                        cGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTGCGATATTGCGCCGTGGATGGAGGAAAGCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCaa      <  5:312499/149‑1 (MQ=255)
                                                                                                                                         gATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAg     <  3:114364/149‑1 (MQ=255)
                                                                                                                                             gtcgtcATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGACACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGACGCTGCCCGATACGCCCATCAAGACg  >  4:293913/1‑148 (MQ=255)
                                                                                                                                               cgtcATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCg                          <  8:211158/122‑1 (MQ=255)
                                                                                                                                               cgtcATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCg                          >  7:211158/1‑122 (MQ=255)
                                                                                                                                                |                                                                                                                                                
GAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACG  >  NC_000913/3577272‑3577560

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