Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,181,664 T→C G140G (GGT→GGC rpoB → RNA polymerase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,181,6640TC98.9% 317.0 / ‑6.0 87G140G (GGT→GGCrpoBRNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (46/40);  total (47/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTGCGCGTTAAACTGCGTCTGGTGATCTATGAGCGCGAAGCGCCGGAAGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGA  >  NC_000913/4181528‑4181798
                                                                                                                                        |                                                                                                                                      
cGCTGCGCGTTAAACTGCGTCTGGTGATCTATGAGCGCGAAGCGCCGGAAGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACt                                                                                                                                     >  1:365808/1‑140 (MQ=255)
  cTGCGCGTTAAACTGCGTCTGGTGATCTATGAGCGCGAAGCGCCGGAAGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGa                                                                                                                                   >  1:910061/1‑140 (MQ=255)
  cTGCGCGTTAAACTGCGTCTGGTGATCTATGAGCGCGAAGCGCCGGAAGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGa                                                                                                                                   <  2:910061/140‑1 (MQ=255)
             aCTGCGTCTGGTGATCTATGAGCGCGAAGCGCCGGAAGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCg                                                                                                                        <  2:365808/140‑1 (MQ=255)
               tGCGTCTGGTGATCTATGAGCGCGAAGCGCCGGAAGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGtt                                                                                                                      >  1:237423/1‑140 (MQ=255)
                       gTGATCTATGAGCGCGAAGCGCCGGAAGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCt                                                                                                              >  2:1309148/1‑140 (MQ=255)
                        tGATCTATGAGCGCGAAGCGCCGGAAGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTg                                                                                                             <  1:1010892/140‑1 (MQ=255)
                         gATCTATGAGCGCGAAGCGCCGGAAGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGc                                                                                                            <  1:1286948/140‑1 (MQ=255)
                             tATGAGCGCGAAGCGCCGGAAGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCg                                                                                                        >  1:205035/1‑140 (MQ=255)
                                gAGCGCGAAGCGCCGGAAGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAg                                                                                                     >  1:1239655/1‑140 (MQ=255)
                                 aGCGCGAAGCGCCGGAAGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGt                                                                                                    >  2:237452/1‑140 (MQ=255)
                                 aGCGCGAAGCGCCGGAAGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGt                                                                                                    >  2:539623/1‑140 (MQ=255)
                                 aGCGCGAAGCGCCGGAAGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGt                                                                                                    >  1:1346322/1‑140 (MQ=255)
                                 aGCGCGAAGCGCCGGAAGGCACCGAAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGt                                                                                                    >  2:621105/1‑140 (MQ=255)
                                  gcgcGAAGCGCCGGAAGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTc                                                                                                   >  2:558734/1‑140 (MQ=255)
                                   cgcgAAGCGCCGGAAGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTcc                                                                                                  >  2:1347663/1‑140 (MQ=255)
                                    gcgAAGCGCCGGAAGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCg                                                                                                 >  1:962196/1‑140 (MQ=255)
                                     cgAAGCGCTGGAAGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCgg                                                                                                >  2:149757/1‑140 (MQ=255)
                                     cgAAGCGCCGGAAGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCgg                                                                                                >  2:1437486/1‑140 (MQ=255)
                                     cgAAGCGCCGGAAGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCgg                                                                                                >  1:717437/1‑140 (MQ=255)
                                     cgAAGCGCCGGAAGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCgg                                                                                                >  2:382766/1‑140 (MQ=255)
                                     cgAAGCGCCGGAAGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCgg                                                                                                >  1:1115436/1‑140 (MQ=255)
                                     cgAAGCGCCGGAAGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCgg                                                                                                >  2:325664/1‑140 (MQ=255)
                                      gAAGCGCCGGAAGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCggg                                                                                               >  2:354387/1‑140 (MQ=255)
                                       aaGCGCCGGAAGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGc                                                                                              >  1:878951/1‑140 (MQ=255)
                                        aGCGCCGGAAGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCg                                                                                             <  2:250091/140‑1 (MQ=255)
                                         gcgcCGGAAGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGt                                                                                            <  2:1354947/140‑1 (MQ=255)
                                          cgcCGGAAGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGtc                                                                                           <  1:409842/140‑1 (MQ=255)
                                          cgcCGGAAGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGtc                                                                                           <  2:962196/140‑1 (MQ=255)
                                           gcCGGAAGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGtct                                                                                          <  2:542407/140‑1 (MQ=255)
                                            ccGGAAGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGtctt                                                                                         >  2:1378771/1‑140 (MQ=255)
                                               gAAGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTcttctt                                                                                      <  2:878951/140‑1 (MQ=255)
                                                aaGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCttt                                                                                     <  2:205035/140‑1 (MQ=255)
                                                    cACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACt                                                                                 >  2:263574/1‑140 (MQ=255)
                                                    cACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACt                                                                                 >  1:498210/1‑140 (MQ=255)
                                                     aCCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTc                                                                                >  1:170309/1‑140 (MQ=255)
                                                       cGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCg                                                                              >  1:646135/1‑140 (MQ=255)
                                                         tAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGAc                                                                            >  2:896182/1‑140 (MQ=255)
                                                          aaaaGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACa                                                                           <  2:81062/140‑1 (MQ=255)
                                                          aaaaGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACa                                                                           <  1:1099411/140‑1 (MQ=255)
                                                          aaaaGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACa                                                                           <  1:237452/140‑1 (MQ=255)
                                                           aaaGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACaa                                                                          <  1:960719/140‑1 (MQ=255)
                                                                 aTTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTaa                                                                    <  2:85181/140‑1 (MQ=255)
                                                                 aTTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTaa                                                                    >  2:730605/1‑140 (MQ=255)
                                                                 aTTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTaa                                                                    >  1:85181/1‑140 (MQ=255)
                                                                 aTTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTaa                                                                    <  1:584327/140‑1 (MQ=255)
                                                                          cAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCActc                                                           >  2:191535/1‑140 (MQ=255)
                                                                             gaagTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTc                                                        >  2:1148886/1‑140 (MQ=255)
                                                                             gaagTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTc                                                        <  2:170309/140‑1 (MQ=255)
                                                                               agTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTc                                                        <  1:356237/138‑1 (MQ=255)
                                                                               agTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTc                                                        >  2:356237/1‑138 (MQ=255)
                                                                                  cTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACTGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTa                                                   <  1:492526/140‑1 (MQ=255)
                                                                                  cTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTa                                                   <  1:382766/140‑1 (MQ=255)
                                                                                  cTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTa                                                   <  1:1437486/140‑1 (MQ=255)
                                                                                  cTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTa                                                   <  1:392675/140‑1 (MQ=255)
                                                                                  cTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTa                                                   <  1:580651/140‑1 (MQ=255)
                                                                                  cTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTa                                                   >  2:178971/1‑140 (MQ=255)
                                                                                  cTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTCTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTa                                                   <  1:191535/140‑1 (MQ=255)
                                                                                    aCATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTaaa                                                 <  1:263574/140‑1 (MQ=255)
                                                                                    aCATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTaaa                                                 <  2:717437/140‑1 (MQ=255)
                                                                                        gggCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGc                                             <  1:180926/140‑1 (MQ=255)
                                                                                        gggCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGc                                             >  2:1024656/1‑140 (MQ=255)
                                                                                           cGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGt                                               >  2:127633/1‑135 (MQ=255)
                                                                                           cGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGt                                               <  1:127633/135‑1 (MQ=255)
                                                                                           cGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGt                                          >  1:867668/1‑140 (MQ=255)
                                                                                           cGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGt                                          >  1:813983/1‑140 (MQ=255)
                                                                                           cGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGt                                          >  2:681463/1‑140 (MQ=255)
                                                                                           cGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCACTCTCCGGGTAAAGTGCTGta                                         >  2:1263223/1‑140 (MQ=255)
                                                                                             aaaTTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGtat                                        <  2:498210/140‑1 (MQ=255)
                                                                                              aaTTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGtata                                       <  2:867668/140‑1 (MQ=255)
                                                                                                     cTCATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAAcgcgcg                                <  2:283966/140‑1 (MQ=255)
                                                                                                       cATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTa                              <  2:1115436/140‑1 (MQ=255)
                                                                                                       cATGACAGACAACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTa                              <  1:730605/140‑1 (MQ=255)
                                                                                                            cagacaACGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTatcatc                         <  2:646135/140‑1 (MQ=255)
                                                                                                                  aCGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTAc                   >  1:617731/1‑140 (MQ=255)
                                                                                                                  aCGGTACCTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTAc                   >  1:849412/1‑140 (MQ=255)
                                                                                                                        ccTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGt             <  1:178971/140‑1 (MQ=255)
                                                                                                                         cTTTGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGtt            <  1:1148886/140‑1 (MQ=255)
                                                                                                                           ttGTTATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTcc          >  1:652636/1‑140 (MQ=255)
                                                                                                                              ttATCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCctgg       >  2:820785/1‑140 (MQ=255)
                                                                                                                                aTCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCctggct     >  2:401991/1‑140 (MQ=255)
                                                                                                                                 tCAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCctggctg    <  1:325664/140‑1 (MQ=255)
                                                                                                                                  cAACGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCctggctgg   <  1:820785/140‑1 (MQ=255)
                                                                                                                                   aaCGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGtt            <  1:903169/130‑1 (MQ=255)
                                                                                                                                   aaCGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGtt            >  2:903169/1‑130 (MQ=255)
                                                                                                                                   aaCGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGa  >  1:141556/1‑140 (MQ=255)
                                                                                                                                   aaCGGCACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGa  >  2:922120/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                      
CGCTGCGCGTTAAACTGCGTCTGGTGATCTATGAGCGCGAAGCGCCGGAAGGCACCGTAAAAGACATTAAAGAACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGA  >  NC_000913/4181528‑4181798

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: