Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,465,588 A→G V198A (GTT→GCT)  treB ← trehalose‑specific PTS enzyme IIBC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,465,5880AG97.9% 326.8 / ‑3.0 95V198A (GTT→GCT) treBtrehalose‑specific PTS enzyme IIBC component
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (45/48);  total (46/49)
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.98e-01

AGCGCCAGTCCGGCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCA  >  NC_000913/4465456‑4465723
                                                                                                                                    |                                                                                                                                       
aGCGCCAGTCCGGCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAgtgt                                                                                                                                  <  1:1196787/140‑1 (MQ=255)
    ccAGTCCGGCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGTCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACg                                                                                                                              <  1:311805/140‑1 (MQ=255)
    ccAGTCCGGCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACg                                                                                                                              >  2:107084/1‑140 (MQ=255)
    ccAGTCCGGCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACg                                                                                                                              <  2:158906/140‑1 (MQ=255)
    ccAGTCCGGCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACg                                                                                                                              <  1:482014/140‑1 (MQ=255)
      aGTCCGGCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGcc                                                                                                                            <  1:418295/140‑1 (MQ=255)
           ggCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCa                                                                                                                       <  1:488421/140‑1 (MQ=255)
            gCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCa                                                                                                                       <  2:55085/139‑1 (MQ=255)
            gCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCa                                                                                                                       >  1:55085/1‑139 (MQ=255)
            gCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCAc                                                                                                                      >  2:1284439/1‑140 (MQ=255)
            gCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCAc                                                                                                                      >  1:982135/1‑140 (MQ=255)
                  aGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATAcc                                                                                                                <  2:1116272/140‑1 (MQ=255)
                  aGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATAcc                                                                                                                <  2:1484353/140‑1 (MQ=255)
                      ccGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCa                                                                                                                       <  2:1135618/129‑1 (MQ=255)
                      ccGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCa                                                                                                                       >  1:1135618/1‑129 (MQ=255)
                      ccGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCa                                                                                                                       >  2:77479/1‑129 (MQ=255)
                      ccGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCa                                                                                                                       <  1:77479/129‑1 (MQ=255)
                         ggATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATc                                                                                                         <  2:441602/140‑1 (MQ=255)
                           aTCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCgg                                                                                                       >  1:833549/1‑140 (MQ=255)
                              aCCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGt                                                                                                    >  2:576424/1‑140 (MQ=255)
                              aCCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGt                                                                                                    >  2:831660/1‑140 (MQ=255)
                              aCCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGt                                                                                                    >  2:273747/1‑140 (MQ=255)
                              aCCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGt                                                                                                    >  2:52082/1‑140 (MQ=255)
                                   cgcCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTgccgc                                                                                               >  1:51176/1‑140 (MQ=255)
                                     ccTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGccc                                                                                             >  2:911586/1‑140 (MQ=255)
                                     ccTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGccc                                                                                             >  1:1621867/1‑140 (MQ=255)
                                        gATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCAtt                                                                                          >  1:1571499/1‑140 (MQ=255)
                                        gATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCAtt                                                                                          >  1:513899/1‑140 (MQ=255)
                                        gATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCAtt                                                                                          >  1:555513/1‑140 (MQ=255)
                                         aTAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCAttt                                                                                         >  1:980258/1‑140 (MQ=255)
                                          tAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCAtttt                                                                                        >  2:428126/1‑140 (MQ=255)
                                             ccTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCAttttttt                                                                                     >  1:472321/1‑140 (MQ=255)
                                                aCTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCAc                                                                                  >  2:672883/1‑140 (MQ=255)
                                                  tttGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCg                                                                                <  1:428126/140‑1 (MQ=255)
                                                  tttGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCg                                                                                <  2:1275043/140‑1 (MQ=255)
                                                  tttGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCg                                                                                >  1:1557724/1‑140 (MQ=255)
                                                   ttGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGc                                                                               >  2:1294124/1‑140 (MQ=255)
                                                      gCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGa                                                                            <  2:633410/140‑1 (MQ=255)
                                                      gCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGa                                                                            <  2:269571/140‑1 (MQ=255)
                                                      gCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGa                                                                            <  2:980258/140‑1 (MQ=255)
                                                      gCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGa                                                                            <  1:1284439/140‑1 (MQ=255)
                                                      gCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGAGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGa                                                                            <  2:1378580/140‑1 (MQ=255)
                                                       cGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGa                                                                            <  1:693830/139‑1 (MQ=255)
                                                       cGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGa                                                                            >  2:693830/1‑139 (MQ=255)
                                                           gCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGc                                                                       >  2:418079/1‑140 (MQ=255)
                                                           gCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGc                                                                       >  2:199326/1‑140 (MQ=255)
                                                           gCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGc                                                                       >  2:88575/1‑140 (MQ=255)
                                                               aaCATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAAt                                                                   >  2:1420599/1‑140 (MQ=255)
                                                                  aTGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATacc                                                                <  2:1571499/140‑1 (MQ=255)
                                                                  aTGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATacc                                                                <  2:555513/140‑1 (MQ=255)
                                                                    gCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGa                                                                            <  2:1314728/126‑1 (MQ=255)
                                                                    gCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGa                                                                            >  1:1314728/1‑126 (MQ=255)
                                                                    gCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATaccga                                                              >  2:1279932/1‑140 (MQ=255)
                                                                      cAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATaccgacc                                                            >  1:24593/1‑140 (MQ=255)
                                                                        aaGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCgg                                                          <  1:52082/140‑1 (MQ=255)
                                                                        aaGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCgg                                                          <  1:88575/140‑1 (MQ=255)
                                                                        aaGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCgg                                                          <  1:1420599/140‑1 (MQ=255)
                                                                         aGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGc                                                         <  2:1621867/140‑1 (MQ=255)
                                                                           tCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAg                                                       <  2:513899/140‑1 (MQ=255)
                                                                            cccACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAgg                                                      <  2:1350093/140‑1 (MQ=255)
                                                                            cccACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAgg                                                      <  1:672883/140‑1 (MQ=255)
                                                                            cccACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAgg                                                      <  2:165362/140‑1 (MQ=255)
                                                                            cccACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAgg                                                      <  2:51176/140‑1 (MQ=255)
                                                                              cacaCTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTa                                                    >  1:224548/1‑140 (MQ=255)
                                                                              cacaCTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTa                                                    >  2:296383/1‑140 (MQ=255)
                                                                                    tCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTagaagaa                                              <  2:224548/140‑1 (MQ=255)
                                                                                    tCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTagaagaa                                              <  2:332616/140‑1 (MQ=255)
                                                                                    tCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTagaagaa                                              <  1:964778/140‑1 (MQ=255)
                                                                                    tCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTagaagaa                                              <  1:1294124/140‑1 (MQ=255)
                                                                                    tCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGGagaagaa                                              >  2:787914/1‑140 (MQ=255)
                                                                                     ccGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCGCGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCgg                                                          >  1:23290/1‑127 (MQ=255)
                                                                                     ccGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCGCGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCgg                                                          <  2:23290/127‑1 (MQ=255)
                                                                                     ccGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAgg                                                      <  2:970463/131‑1 (MQ=255)
                                                                                     ccGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAgg                                                      >  1:970463/1‑131 (MQ=255)
                                                                                      cGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTagaagaaga                                            <  1:418079/140‑1 (MQ=255)
                                                                                            ctgctgCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCtt                                      <  1:199326/140‑1 (MQ=255)
                                                                                             tgctgcCCGAGCAGATAAGCGTTCATCAGCTTTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTc                                     >  1:612411/1‑140 (MQ=255)
                                                                                               ctgcCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCAc                                   >  1:172897/1‑140 (MQ=255)
                                                                                                      aGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGCCCGGCAGGGAGAAGAAGAACGCCTTACCCGTCaa                            >  1:393099/1‑140 (MQ=255)
                                                                                                      aGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCaa                            <  1:1279932/140‑1 (MQ=255)
                                                                                                      aGCAGATAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCaa                            <  1:831660/140‑1 (MQ=255)
                                                                                                           aTAAGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACcaca                       <  2:982135/140‑1 (MQ=255)
                                                                                                             aaGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGa                     <  2:24593/140‑1 (MQ=255)
                                                                                                             aaGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGa                     <  2:73826/140‑1 (MQ=255)
                                                                                                             aaGCGTTCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGa                     <  1:911586/140‑1 (MQ=255)
                                                                                                                  ttCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAATAAGATCGCTTCACCGATCAACCACAGAAAATc                >  1:508438/1‑140 (MQ=255)
                                                                                                                  ttCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGa                     >  1:899654/1‑135 (MQ=255)
                                                                                                                  ttCATCAGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGa                     <  2:899654/135‑1 (MQ=255)
                                                                                                                     atcaGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACTGCTGAACAGCAAATACCGACCGGCAGGCAGAAAAAGATCGCTTCACCGGACAACCACAGAAAATCGTa             >  1:1436044/1‑140 (MQ=255)
                                                                                                                     atcaGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTa             >  1:755746/1‑140 (MQ=255)
                                                                                                                     atcaGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTa             >  1:429004/1‑140 (MQ=255)
                                                                                                                        aGCTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGAt          <  1:273747/140‑1 (MQ=255)
                                                                                                                          cTGTGGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCg        >  2:1352641/1‑140 (MQ=255)
                                                                                                                             tgGAGAAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGttt     >  1:133189/1‑140 (MQ=255)
                                                                                                                                agaAGCCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCa  <  2:489199/140‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                       
AGCGCCAGTCCGGCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCA  >  NC_000913/4465456‑4465723

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