Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,577,428 Δ2 bp coding (164‑165/528 nt) gntK ← gluconate kinase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,577,4280T.98.7% 344.4 / ‑1.8 78coding (165/528 nt)gntKgluconate kinase 2
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (35/42);  total (36/42)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01
*NC_0009133,577,4290T.98.7% 344.4 / ‑1.8 78coding (164/528 nt)gntKgluconate kinase 2
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (35/42);  total (36/42)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

TTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATG  >  NC_000913/3577288‑3577566
                                                                                                                                            ||                                                                                                                                         
ttACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACg                                                                                                                                  <  3:173806/149‑1 (MQ=255)
       cACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCAtttt                                                                                                                           >  3:358717/1‑149 (MQ=255)
         cGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCAtttttt                                                                                                                         <  2:285377/149‑1 (MQ=255)
         cGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACAC‑TTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTc                                                                                                                        <  8:222042/149‑1 (MQ=255)
         cGCAGCAAGGCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCAtttttt                                                                                                                         <  4:331961/149‑1 (MQ=255)
          gcagcaAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTc                                                                                                                        <  4:249207/149‑1 (MQ=255)
            agcaAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGa                                                                                                                      <  8:288456/149‑1 (MQ=255)
                aGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATAtt                                                                                                                  <  7:356406/149‑1 (MQ=255)
                  tCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCAtttttt                                                                                                                         >  1:119498/1‑140 (MQ=255)
                  tCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCAtttttt                                                                                                                         <  2:119498/140‑1 (MQ=255)
                  tCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTgc                                                                                                                >  2:191352/1‑149 (MQ=255)
                          gtgtTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGtgga                                                                                                        <  3:122251/149‑1 (MQ=255)
                           tgtTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGtggat                                                                                                       >  1:274557/1‑149 (MQ=255)
                             ttttttCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTGTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGtggatg                                                                                                      <  4:289916/148‑1 (MQ=255)
                             ttttttCAATGCAGAACAGACGATCAGCGACACCTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCGAGTCGCCATTTTTTCGATATTGCGCCGtggatg                                                                                                      <  2:114396/148‑1 (MQ=255)
                              tttttCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGtggatg                                                                                                      >  1:356585/1‑147 (MQ=255)
                                    aaTGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACTCGGCGTCGTTCAGCGCATGCAACCACGGTTTGCGATCTTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGaaa                                                                                               >  3:147942/1‑148 (MQ=255)
                                    aaTGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGTGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGaaa                                                                                               <  7:344948/148‑1 (MQ=255)
                                    aaTGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAAt                                                                                              >  3:142196/1‑149 (MQ=255)
                                         aGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCa                                                                                         >  8:215059/1‑149 (MQ=255)
                                           aaCAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCAt                                                                                        >  4:286461/1‑148 (MQ=255)
                                           aaCAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCAt                                                                                        >  2:350974/1‑148 (MQ=255)
                                           aaCAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATc                                                                                       <  4:142196/149‑1 (MQ=255)
                                            aCAGAAGATAAGAGACACTTTATTAGTTCGCTGCATATCAAACGAGGCGTCTTTCAGCGCCTGCAACCACGGTGTGCGATAGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTATATATTGCGCCGTGGATGGAGGAAATCGCCATc                                                                                       <  4:147942/148‑1 (MQ=255)
                                               gACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGa                                                                                   <  5:412243/149‑1 (MQ=255)
                                                aCGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGaa                                                                                  <  7:215059/149‑1 (MQ=255)
                                                  gATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAAc                                                                                >  5:319007/1‑149 (MQ=255)
                                                   aTCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑GGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAAcg                                                                               >  5:228324/1‑149 (MQ=255)
                                                    tCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAAcgc                                                                              <  2:356585/149‑1 (MQ=255)
                                                    tCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAAcgc                                                                              >  8:151919/1‑149 (MQ=255)
                                                      aGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAAcgcg                                                                             <  6:281890/148‑1 (MQ=255)
                                                        cGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGata                                                                                                                    >  5:174567/1‑107 (MQ=255)
                                                        cGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGata                                                                                                                    <  6:174567/107‑1 (MQ=255)
                                                        cGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCa                                                                          >  3:244399/1‑149 (MQ=255)
                                                            acTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCAAGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATg                                                                        >  5:173619/1‑147 (MQ=255)
                                                             cTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAg                                                                     >  5:96554/1‑149 (MQ=255)
                                                                tattaGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAAGCGCCATCAAGAAACGCGGCATGAAGt                                                                    <  1:191352/147‑1 (MQ=255)
                                                                  ttaGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGa                                                                 <  7:88109/148‑1 (MQ=255)
                                                                  ttaGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGAt                                                                <  2:274557/149‑1 (MQ=255)
                                                                      tGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTgcgc                                                                                                              <  4:285977/99‑1 (MQ=255)
                                                                      tGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTgcgc                                                                                                              >  3:285977/1‑99 (MQ=255)
                                                                      tGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCa                                                                                         <  6:37773/120‑1 (MQ=255)
                                                                      tGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCa                                                                                         >  5:37773/1‑120 (MQ=255)
                                                                      tGCGCTGCATAGCAAACGCGGCGGCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCCGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCACCATCAAGAAACGCGGCATGAAGGTGATgcgc                                                            <  6:355771/149‑1 (MQ=255)
                                                                           tGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACtt                                                       >  6:270396/1‑149 (MQ=255)
                                                                           tGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACtt                                                       >  8:408935/1‑149 (MQ=255)
                                                                              aTAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCAc                                                    >  5:91491/1‑149 (MQ=255)
                                                                              aTAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCAc                                                    <  7:408935/149‑1 (MQ=255)
                                                                                aGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGtt                                                                   <  1:23835/132‑1 (MQ=255)
                                                                                aGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGtt                                                                   >  2:23835/1‑132 (MQ=255)
                                                                                      cgcgGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGAcc                                            >  8:289456/1‑149 (MQ=255)
                                                                                            gtcgtTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGAt                                      <  3:72078/149‑1 (MQ=255)
                                                                                             tcgtTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGAtt                                     <  7:289456/149‑1 (MQ=255)
                                                                                                ttCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGc                                  >  2:84963/1‑149 (MQ=255)
                                                                                                ttCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGc                                  >  6:87797/1‑149 (MQ=255)
                                                                                                 tCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGcc                                 >  4:78076/1‑149 (MQ=255)
                                                                                                     cgcCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTg                             <  3:78076/149‑1 (MQ=255)
                                                                                                       ccTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGcc                           >  8:181432/1‑149 (MQ=255)
                                                                                                        cTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAACTCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGccc                          <  7:181432/149‑1 (MQ=255)
                                                                                                        cTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCCCTTCACTGGCGACCGCAGATTTGCCGCTGccc                          >  8:148788/1‑149 (MQ=255)
                                                                                                        cTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGccc                          >  5:363359/1‑149 (MQ=255)
                                                                                                         tGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATg                                                               <  2:17247/111‑1 (MQ=255)
                                                                                                         tGCAACCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATg                                                               >  1:17247/1‑111 (MQ=255)
                                                                                                            aaCCACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATa                      <  6:319007/149‑1 (MQ=255)
                                                                                                              ccACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACg                    >  1:181288/1‑149 (MQ=255)
                                                                                                               cACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGc                   <  7:151919/149‑1 (MQ=255)
                                                                                                               cACGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGc                   <  6:91491/149‑1 (MQ=255)
                                                                                                                aCGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGc                                                                           <  7:257576/92‑1 (MQ=255)
                                                                                                                aCGGTTTGCGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGc                                                                           >  8:257576/1‑92 (MQ=255)
                                                                                                                       gCGATCGTCGTCATTCAGTGG‑‑CGCCAGACTCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCaa            <  5:87797/148‑1 (MQ=255)
                                                                                                                        cGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCg                         <  3:71902/134‑1 (MQ=255)
                                                                                                                        cGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCg                         >  4:71902/1‑134 (MQ=255)
                                                                                                                        cGATCGTCGTCATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGa          <  4:358717/149‑1 (MQ=255)
                                                                                                                          aTCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGa          >  1:10751/1‑149 (MQ=255)
                                                                                                                             gtcgtcATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTaa     <  1:84963/149‑1 (MQ=255)
                                                                                                                                gtcATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGTCGACCTCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAAtg  <  4:244399/149‑1 (MQ=255)
                                                                                                                                 tcATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAg                                                                                    <  6:249360/66‑1 (MQ=255)
                                                                                                                                 tcATTCAGTGG‑‑CGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAg                                                                                    >  5:249360/1‑66 (MQ=255)
                                                                                                                                            ||                                                                                                                                         
TTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATG  >  NC_000913/3577288‑3577566

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