Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,079,856 C→G 5.7% A594P (GCT→CCT)  tktA ← transketolase 1, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,079,8560CG 5.7% 315.1 / 9.2 87A594P (GCT→CCT) tktAtransketolase 1, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base C (43/39);  new base G (3/2);  total (46/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.42e-01

AGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACC  >  NC_000913/3079708‑3079994
                                                                                                                                                    |                                                                                                                                          
agcTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAagcagc                                                                                                                                            >  4:626807/1‑149 (MQ=255)
agcTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAagcagc                                                                                                                                            >  1:457751/1‑149 (MQ=255)
 gcTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAgcagca                                                                                                                                           <  3:397201/149‑1 (MQ=255)
      gCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTg                                                                                                                                      <  4:381215/149‑1 (MQ=255)
          gAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTg                                                                                                                                  <  1:599700/149‑1 (MQ=255)
           aGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGt                                                                                                                                 <  4:537001/149‑1 (MQ=255)
            gCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGGATGCTGCttacc                                                                                                                                <  3:407304/149‑4 (MQ=255)
             cAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGc                                                                                                                                     <  3:358472/143‑1 (MQ=255)
             cAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGc                                                                                                                                     >  4:358472/1‑143 (MQ=255)
                aTTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAAt                                                                                                                            >  2:232497/1‑149 (MQ=255)
                aTTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAAt                                                                                                                            <  3:626807/149‑1 (MQ=255)
                aTTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAAt                                                                                                                            <  1:690945/149‑1 (MQ=255)
                  tCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGc                                                                                                                          >  3:235273/1‑149 (MQ=255)
                  tCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGc                                                                                                                          >  3:512770/1‑149 (MQ=255)
                       gAAGGTGGTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGc                                                                                                                                     >  4:489699/1‑133 (MQ=255)
                       gAAGGTGGTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGc                                                                                                                                     <  3:489699/133‑1 (MQ=255)
                                cATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCa                                                                                                            <  2:409049/149‑1 (MQ=255)
                                 aTACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGGATGCTGCttaccgt                                                                                                                              <  3:223840/130‑6 (MQ=255)
                                 aTACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGGATGCTGCttaccgt                                                                                                                              >  4:223840/1‑125 (MQ=255)
                                 aTACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCAt                                                                                                           >  4:285060/1‑149 (MQ=255)
                                  tACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATg                                                                                                          >  2:624166/1‑149 (MQ=255)
                                    ccGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGa                                                                                                        <  4:512770/149‑1 (MQ=255)
                                     cgacgaTAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGac                                                                                                       >  1:644513/1‑149 (MQ=255)
                                               aCCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCAcgcgcg                                                                                             <  2:194936/149‑1 (MQ=255)
                                                    tcaagCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACGGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTc                                                                                        <  4:377951/145‑1 (MQ=255)
                                                    tCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTc                                                                                        <  4:104865/149‑1 (MQ=255)
                                                     cAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCa                                                                                       <  2:680000/149‑1 (MQ=255)
                                                        gCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGc                                                                                    >  2:70787/1‑149 (MQ=255)
                                                           aaCATACGTGTACCAGTAGTGGAAAGTACCCGCTTCTACGGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCtt                                                                                 >  4:456603/1‑149 (MQ=255)
                                                              aTACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGt                                                                                                                        >  3:28530/1‑107 (MQ=255)
                                                              aTACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGt                                                                                                                        <  4:28530/107‑1 (MQ=255)
                                                                 cTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGcag                                                                           >  1:671107/1‑149 (MQ=255)
                                                                 cTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGcag                                                                           <  3:67230/149‑1 (MQ=255)
                                                                  ttGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCa                                                                                       <  3:599325/136‑1 (MQ=255)
                                                                  ttGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCa                                                                                       >  4:33601/1‑136 (MQ=255)
                                                                  ttGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCa                                                                                       >  4:599325/1‑136 (MQ=255)
                                                                  ttGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCa                                                                                       <  3:33601/136‑1 (MQ=255)
                                                                  ttGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGcagt                                                                          >  3:565279/1‑149 (MQ=255)
                                                                   tGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGcagtc                                                                         <  2:671107/149‑1 (MQ=255)
                                                                    gTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGcagtca                                                                        >  4:35517/1‑149 (MQ=255)
                                                                      aCCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGcagtcagt                                                                      <  1:387145/149‑1 (MQ=255)
                                                                       ccAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGtt                                                                     <  1:454461/149‑1 (MQ=255)
                                                                        cAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGttt                                                                    >  1:626868/1‑149 (MQ=255)
                                                                        cAGTAGTCAGCAATACCCACTTGTAGAGCAACGCGTGCAGTAACCGCTTTCAGCAGTACCGATTCACGGTATGCAGCATCCTGCTTGTCAAATGCGTCGGTGGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGttt                                                                    >  4:581508/1‑149 (MQ=255)
                                                                         agtagtCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGtttt                                                                   >  1:228115/1‑149 (MQ=255)
                                                                           tagtCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTc                                                                 <  1:70787/149‑1 (MQ=255)
                                                                             gtCAGCAATACCCGCTTCAACAGCAACGCGTGCAGTAACCACTTTCGGCAGTACGGATTCACGGTAAGCAGGATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGt                                                               >  1:200848/1‑149 (MQ=255)
                                                                              tCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTATCCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTa                                                              <  4:467523/149‑1 (MQ=255)
                                                                              tCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGca                                                                            <  4:526576/135‑1 (MQ=255)
                                                                              tCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGca                                                                            >  3:526576/1‑135 (MQ=255)
                                                                              tCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTa                                                              >  1:469142/1‑149 (MQ=255)
                                                                              tCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTa                                                              >  2:488113/1‑149 (MQ=255)
                                                                                aGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAgg                                                            >  2:510705/1‑149 (MQ=255)
                                                                                aGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAgg                                                            <  2:333868/149‑1 (MQ=255)
                                                                                 gCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGgc                                                           <  1:232497/149‑1 (MQ=255)
                                                                                 gCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGgc                                                           >  3:365403/1‑149 (MQ=255)
                                                                                   aaTACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGgcag                                                         >  1:35320/1‑149 (MQ=255)
                                                                                   aaTACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGgcag                                                         >  4:477941/1‑149 (MQ=255)
                                                                                        ccGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACa                                                    >  3:265930/1‑149 (MQ=255)
                                                                                        ccGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACa                                                    <  1:436192/149‑1 (MQ=255)
                                                                                        ccGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACa                                                    >  4:381776/1‑149 (MQ=255)
                                                                                          gCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGc                                                  >  2:464876/1‑149 (MQ=255)
                                                                                          gCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGc                                                  <  4:365403/149‑1 (MQ=255)
                                                                                            ttCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCa                                                <  1:344593/149‑1 (MQ=255)
                                                                                              cTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGt                                              >  1:687753/1‑149 (MQ=255)
                                                                                               tACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGtt                                             >  3:14599/1‑149 (MQ=255)
                                                                                                  aGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAg                                                             <  3:590831/130‑1 (MQ=255)
                                                                                                  aGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAg                                                             >  4:590831/1‑130 (MQ=255)
                                                                                                  aGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCaa                                          <  2:228115/149‑1 (MQ=255)
                                                                                                   gCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCttt                                                                                         <  2:48517/101‑1 (MQ=255)
                                                                                                   gCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCttt                                                                                         >  1:48517/1‑101 (MQ=255)
                                                                                                   gCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAAc                                         <  3:285060/149‑1 (MQ=255)
                                                                                                      aCGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTc                                      >  4:627643/1‑149 (MQ=255)
                                                                                                          gTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGaa                                  >  2:532006/1‑149 (MQ=255)
                                                                                                               gTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGt                             >  4:673691/1‑149 (MQ=255)
                                                                                                                  aCCGCTTTCGTCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCGTGGACAGCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTATTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCCGTAGc                          >  4:448421/1‑149 (MQ=255)
                                                                                                                  aCCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGc                          >  1:46918/1‑149 (MQ=255)
                                                                                                                    cGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGca                                                                            >  2:630708/1‑97 (MQ=255)
                                                                                                                    cGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGca                                                                            <  1:630708/97‑1 (MQ=255)
                                                                                                                    cGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGa                        <  4:222883/149‑1 (MQ=255)
                                                                                                                    cGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGa                        >  3:117233/1‑149 (MQ=255)
                                                                                                                           ggCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAAt                 <  2:687753/149‑1 (MQ=255)
                                                                                                                            gCAGTACGGATTCACGGTAAGCAGCATCCTGCGTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGTGATGAAAATc                >  4:591247/1‑149 (MQ=255)
                                                                                                                             cAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCa               <  3:477941/149‑1 (MQ=255)
                                                                                                                               gTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGt             >  3:415362/1‑149 (MQ=255)
                                                                                                                                tACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGtt            >  4:386874/1‑149 (MQ=255)
                                                                                                                                  cGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTcc          <  2:46918/149‑1 (MQ=255)
                                                                                                                                   ggATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCg         <  1:624166/149‑1 (MQ=255)
                                                                                                                                        cACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGa    <  4:235273/149‑1 (MQ=255)
                                                                                                                                          tgacaaGCAGGATGCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGAcc  >  2:126677/5‑149 (MQ=255)
                                                                                                                                                 cagcaTCCTTCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGgcagc                                                        >  3:368195/1‑88 (MQ=255)
                                                                                                                                                 cagcaTCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGgcagc                                                        <  4:368195/88‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                          
AGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGTAGGCAGCAACAGCCAGTTCAACTTCTGAACCGGTAGCGATGAAAATCAGTTCCGGCTGACC  >  NC_000913/3079708‑3079994

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: