Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,081,665 C→T 6.8% intergenic (‑30/‑248) tktA ← / → loiP transketolase 1, thiamine triphosphate‑binding/Phe‑Phe periplasmic metalloprotease, OM lipoprotein; low salt‑inducible; Era‑binding heat shock protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,081,6650CT 6.8% 291.4 / 9.1 88intergenic (‑30/‑248)tktA/loiPtransketolase 1, thiamine triphosphate‑binding/Phe‑Phe periplasmic metalloprotease, OM lipoprotein; low salt‑inducible; Era‑binding heat shock protein
Reads supporting (aligned to +/- strand):  ref base C (41/41);  new base T (3/3);  total (44/44)
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.85e-01

TCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAAT  >  NC_000913/3081523‑3081806
                                                                                                                                              |                                                                                                                                             
tCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCt                                                                                                                                          <  2:1101762/148‑1 (MQ=255)
      aTGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGa                                                                                                                                    <  1:203151/148‑1 (MQ=255)
      aTGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGa                                                                                                                                    >  2:203151/1‑148 (MQ=255)
         tCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACtt                                                                                                                                 >  1:140758/1‑148 (MQ=255)
          cAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTCTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGACGATGAAGGGCACGCCCTTAACGACTTg                                                                                                                                <  4:331401/148‑1 (MQ=255)
               aTACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGccc                                                                                                                                           >  1:1111005/1‑132 (MQ=255)
               aTACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGccc                                                                                                                                           <  2:1111005/132‑1 (MQ=255)
                      tAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTa                                                                                                                                        >  2:246139/1‑128 (MQ=255)
                      tAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTa                                                                                                                                        <  1:246139/128‑1 (MQ=255)
                                  ggTGACCGGATTTGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGc                                                                                                                                             >  2:542757/1‑111 (MQ=255)
                                  ggTGACCGGATTTGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGc                                                                                                                                             <  1:542757/111‑1 (MQ=255)
                                        cGGATTTGGCTTTCTGTACTGCGTCCATGCCCAGCGCACGAATAGCATTGGCGAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGaaa                                                                                                  >  4:44570/1‑148 (MQ=255)
                                             ttGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCCTTGGAAAGCTCTTTACGTGCGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTg                                                                                             >  2:751203/1‑148 (MQ=255)
                                             ttGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTg                                                                                             <  2:140758/148‑1 (MQ=255)
                                                gCTTTCTGTACTGCGTCCATGCCCAGCGCACGAATAGCATTGGCGAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCca                                                                                          <  3:44570/148‑1 (MQ=255)
                                                   ttCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCGTGCCCttc                                                                                                                                        >  2:531770/1‑98 (MQ=255)
                                                   ttCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCGTGCCCttc                                                                                                                                        <  1:531770/99‑2 (MQ=255)
                                                   ttCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGa                                                                                                                                    >  2:354104/1‑103 (MQ=255)
                                                   ttCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGa                                                                                                                                    <  1:354104/103‑1 (MQ=255)
                                                       gTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTg                                                                                                                >  1:500434/1‑119 (MQ=255)
                                                       gTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTg                                                                                                                <  2:500434/119‑1 (MQ=255)
                                                        tACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTa                                                                                  <  1:751203/148‑1 (MQ=255)
                                                         aCTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGAAATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTat                                                                                 >  1:997090/1‑147 (MQ=255)
                                                            gCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTg                                                                                                                                >  2:917848/1‑98 (MQ=255)
                                                            gCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTg                                                                                                                                <  1:917848/98‑1 (MQ=255)
                                                            gCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTg                                                                                             >  3:203466/1‑133 (MQ=255)
                                                            gCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTg                                                                                             <  4:203466/133‑1 (MQ=255)
                                                                   tGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTg                                                                                                                                >  3:887272/1‑91 (MQ=255)
                                                                   tGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTg                                                                                                                                <  4:887272/91‑1 (MQ=255)
                                                                     cTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTgacgac                                                                                                                           >  1:963237/1‑94 (MQ=255)
                                                                     cTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACACTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTgacgac                                                                                                                           <  2:963237/94‑1 (MQ=255)
                                                                      tCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCa                                                                             >  4:314750/1‑139 (MQ=255)
                                                                      tCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCa                                                                             <  3:314750/139‑1 (MQ=255)
                                                                       cAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCGTGCCCttca                                                                                                                                       <  2:216575/80‑3 (MQ=255)
                                                                       cAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCGTGCCCttca                                                                                                                                       >  1:216575/1‑78 (MQ=255)
                                                                       cAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTgac                                                                                                                              <  2:719529/89‑1 (MQ=255)
                                                                       cAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTgac                                                                                                                              >  1:719529/1‑89 (MQ=255)
                                                                       cAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGt                                                                   >  4:197123/1‑148 (MQ=255)
                                                                               gAATAGCATTGGCAAGCTCTTTACGTGAGGAAATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACa                                                           <  2:997090/148‑1 (MQ=255)
                                                                                aaTAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTa                                                                                                                                        <  4:500246/70‑1 (MQ=255)
                                                                                aaTAGCATTGGCAAGCTCTTTACGTGAGGACATTCTGACTCCAGATCGGATGATGAAGGGCACGCCCTTa                                                                                                                                        >  3:500246/1‑70 (MQ=255)
                                                                                   aGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAg                                                                                                                         <  1:68532/82‑1 (MQ=255)
                                                                                   aGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAg                                                                                                                         >  2:68532/1‑82 (MQ=255)
                                                                                   aGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGa                                                       <  3:197123/148‑1 (MQ=255)
                                                                                           cAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCcaaca                                                                                       >  1:654277/1‑108 (MQ=255)
                                                                                           cAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCcaaca                                                                                       <  2:654277/108‑1 (MQ=255)
                                                                                           cAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCAt                                                                >  1:830823/1‑131 (MQ=255)
                                                                                           cAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCAt                                                                <  2:830823/131‑1 (MQ=255)
                                                                                           cAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCtt                                               >  1:1014892/1‑148 (MQ=255)
                                                                                             aGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGt                                                                                                                    <  2:646683/77‑1 (MQ=255)
                                                                                             aGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGt                                                                                                                    >  1:646683/1‑77 (MQ=255)
                                                                                             aGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCcg                                                                        >  1:367298/1‑121 (MQ=255)
                                                                                             aGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCcg                                                                        <  2:367298/121‑1 (MQ=255)
                                                                                             aGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGa                                                       <  3:169269/138‑1 (MQ=255)
                                                                                             aGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGa                                                       >  4:169269/1‑138 (MQ=255)
                                                                                             aGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCtttt                                             >  4:948283/1‑148 (MQ=255)
                                                                                                   ttACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTg                                                                                             >  2:902400/1‑94 (MQ=255)
                                                                                                   ttACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTg                                                                                             <  1:902400/94‑1 (MQ=255)
                                                                                                   ttACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATg                                                         <  1:377723/130‑1 (MQ=255)
                                                                                                   ttACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATg                                                         >  2:377723/1‑130 (MQ=255)
                                                                                                    tACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCtt                                      <  3:948283/148‑1 (MQ=255)
                                                                                                             aCATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTa                                                                                                                                        <  2:855429/41‑1 (MQ=255)
                                                                                                             aCATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTa                                                                                                                                        >  1:855429/1‑41 (MQ=255)
                                                                                                             aCATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCcgcg                                                                      <  2:213270/107‑1 (MQ=255)
                                                                                                             aCATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCcgcg                                                                      >  1:213270/1‑107 (MQ=255)
                                                                                                             aCATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACAt                                                          <  2:458888/119‑1 (MQ=255)
                                                                                                             aCATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACAt                                                          >  1:458888/1‑119 (MQ=255)
                                                                                                             aCATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAAtct                             >  2:826108/1‑148 (MQ=255)
                                                                                                                ttttGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGt                                                              <  3:581256/112‑1 (MQ=255)
                                                                                                                ttttGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGt                                                              >  4:581256/1‑112 (MQ=255)
                                                                                                                ttttGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCt                                       >  4:1024044/1‑135 (MQ=255)
                                                                                                                ttttGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCt                                       <  3:1024044/135‑1 (MQ=255)
                                                                                                                 tttGACTCCAGATCGGATGATGAAGGGCACGCCCTTAGCGACTTGACGACAGCGCGTTTTGGGCTACGc                                                                                                        <  4:975982/69‑1 (MQ=255)
                                                                                                                 tttGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGc                                                                                                        >  3:975982/1‑69 (MQ=255)
                                                                                                                   tGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGAt                       >  1:956563/1‑148 (MQ=255)
                                                                                                                        ccAGATCGGATGATGACGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTg                                                                                                                <  4:694589/54‑1 (MQ=255)
                                                                                                                        ccAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTg                                                                                                                >  3:694589/1‑54 (MQ=255)
                                                                                                                        ccAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCa                     <  2:463606/145‑1 (MQ=255)
                                                                                                                        ccAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCa                     >  1:463606/1‑145 (MQ=255)
                                                                                                                          aGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCa                                                                 >  2:559220/1‑99 (MQ=255)
                                                                                                                          aGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCa                                                                 <  1:559220/99‑1 (MQ=255)
                                                                                                                             tCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGc                                        >  3:357816/1‑121 (MQ=255)
                                                                                                                             tCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGc                                        <  4:357816/121‑1 (MQ=255)
                                                                                                                                 atgatgAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGtt         >  3:204543/1‑148 (MQ=255)
                                                                                                                                  tgatgaAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTg        <  2:1014892/148‑1 (MQ=255)
                                                                                                                                      taaGGGCGTGCCCTTAACGACTTGACGACAGCGCGTTTTggg                                                                                                              <  3:627014/41‑1 (MQ=255)
                                                                                                                                      taaGGGCGTGCCCTTAACGACTTGACGACAGCGCGTTTTggg                                                                                                              >  4:627014/2‑42 (MQ=255)
                                                                                                                                        aGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAAt  <  4:204543/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                             
TCGGCAATGTCAGCCATACCCATAGGGGCACCCGGGTGACCGGATTTGGCTTTCTGTACTGCGTCCATGCTCAGCGCACGAATAGCATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAAT  >  NC_000913/3081523‑3081806

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