Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,186,605 A→C H419P (CAC→CCC)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,186,6050AC97.3% 267.8 / ‑3.0 74H419P (CAC→CCC) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (50/21);  total (51/23)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

GGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGT  >  NC_000913/4186471‑4186752
                                                                                                                                      |                                                                                                                                                   
ggAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGa                                                                                                                                       <  1:514762/149‑1 (MQ=255)
  aGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCg                                                                                                                                   >  1:434682/1‑151 (MQ=255)
      gTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTg                                                                                                                                 <  2:1116888/149‑1 (MQ=255)
          aaaCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCg                                                                                                                            >  1:92851/1‑150 (MQ=255)
              cGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGActct                                                                                                                       >  2:54528/1‑151 (MQ=255)
                  cATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCa                                                                                                                    <  1:492771/150‑1 (MQ=255)
                      tACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGGCGGTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGt                                                                                                                >  2:606966/1‑150 (MQ=255)
                       aCGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCt                                                                                                              <  1:601626/151‑1 (MQ=255)
                           cAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGt                                                                                                          >  1:315043/1‑151 (MQ=255)
                            aaGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTa                                                                                                         <  2:265743/151‑1 (MQ=255)
                               cTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATcc                                                                                                      >  1:469874/1‑151 (MQ=255)
                                tGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATcc                                                                                                      >  1:921191/1‑150 (MQ=255)
                                 ggAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAg                                                                                                    <  2:448801/151‑1 (MQ=255)
                                    aCTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGc                                                                                                  >  1:950267/1‑150 (MQ=255)
                                    aCTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGc                                                                                                  >  1:1059003/1‑150 (MQ=255)
                                       gCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCAttt                                                                                              <  2:814551/151‑1 (MQ=255)
                                        cGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTg                                                                                             >  2:956525/1‑151 (MQ=255)
                                          tGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGa                                                                                            >  2:937993/1‑150 (MQ=255)
                                                  cTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTAc                                                                                    >  1:587862/1‑150 (MQ=255)
                                                   tACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTg                                                                                  >  2:392566/1‑151 (MQ=255)
                                                    accaccATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACt                                                                                   <  1:645496/149‑1 (MQ=255)
                                                       gccaTTAAAGCTGTGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACGCTCGGTACTGCTGAACCGTGCACCGACTCTGTACCGTCTGGGTATCCAGGCATTGCAACCGGTaat                                                                                   <  2:725045/145‑3 (MQ=255)
                                                       accaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCg                                                                              <  2:345641/151‑1 (MQ=255)
                                                        ccaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCg                                                                              >  2:344896/1‑150 (MQ=255)
                                                        ccaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCg                                                                              >  1:889975/1‑150 (MQ=255)
                                                         caTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGa                                                                             >  1:499084/1‑150 (MQ=255)
                                                         caTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGa                                                                             >  2:521840/1‑150 (MQ=255)
                                                          aTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGa                                                                             >  1:682242/1‑149 (MQ=255)
                                                           ttAAAGGTGCGAAGAAAATGGTTGACCGCGAAGAAGCTGTCGATTGGGATATCCTGGACGAGGTTATCTGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCATCGTCTGGGAATCCGGGCATTTGAACCGGTACTGATCGAAg                                                                           >  2:25352/1‑150 (MQ=255)
                                                            tAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGt                                                                         >  2:777192/1‑151 (MQ=255)
                                                            tAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTCGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAgg                                                                          >  1:266111/1‑150 (MQ=255)
                                                                gCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAg                                                                     <  2:112037/151‑1 (MQ=255)
                                                                 cTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGc                                                                    >  2:594306/1‑151 (MQ=255)
                                                                 cTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGc                                                                    >  2:560011/1‑151 (MQ=255)
                                                                   gCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCt                                                                   <  2:375646/150‑1 (MQ=255)
                                                                        gaaAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAg                                                             <  2:1023014/151‑1 (MQ=255)
                                                                         aaaaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAg                                                             <  2:260683/150‑1 (MQ=255)
                                                                         aaaaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGc                                                            <  2:950242/151‑1 (MQ=255)
                                                                         aaaaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGc                                                            >  2:1052784/1‑151 (MQ=255)
                                                                          aaaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCTTTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCt                                                           >  1:340482/1‑151 (MQ=255)
                                                                           aaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCt                                                           >  2:913042/1‑150 (MQ=255)
                                                                           aaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTg                                                          >  2:185966/1‑151 (MQ=255)
                                                                           aaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTg                                                          >  1:566547/1‑151 (MQ=255)
                                                                              ggTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCAc                                                       >  2:673205/1‑151 (MQ=255)
                                                                                 tGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCg                                                    >  2:928368/1‑151 (MQ=255)
                                                                                  gAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCg                                                    <  1:775876/150‑1 (MQ=255)
                                                                                  gAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCg                                                    >  1:523997/1‑150 (MQ=255)
                                                                                   aGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCt                                                  >  1:130278/1‑151 (MQ=255)
                                                                                        gaagaaGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCTTTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGttt                                             >  1:423377/1‑151 (MQ=255)
                                                                                               cTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGc                                      >  2:699439/1‑151 (MQ=255)
                                                                                                 gTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCat                                    <  1:247187/151‑1 (MQ=255)
                                                                                                 gTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCat                                    >  2:558346/1‑151 (MQ=255)
                                                                                                 gTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCat                                    >  1:76313/1‑151 (MQ=255)
                                                                                                   cGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCatat                                  >  1:255231/1‑151 (MQ=255)
                                                                                                       tGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAAc                               <  2:1061720/150‑1 (MQ=255)
                                                                                                       tGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACg                              >  2:875673/1‑151 (MQ=255)
                                                                                                        gggATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACg                              <  2:526800/150‑1 (MQ=255)
                                                                                                          gATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCg                           >  2:1057509/1‑151 (MQ=255)
                                                                                                             atCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGAc                         <  2:434670/150‑1 (MQ=255)
                                                                                                             atCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACt                        >  1:726475/1‑151 (MQ=255)
                                                                                                             atCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACt                        >  1:70347/1‑151 (MQ=255)
                                                                                                              tCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACtt                       >  1:279916/1‑151 (MQ=255)
                                                                                                                cTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCg                     >  2:405784/1‑151 (MQ=255)
                                                                                                                  ggACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGa                    <  2:809045/150‑1 (MQ=255)
                                                                                                                  ggACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGAt                   >  1:1012415/1‑151 (MQ=255)
                                                                                                                    aCGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATgg                 >  1:388962/1‑151 (MQ=255)
                                                                                                                     cGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGt                >  2:450109/1‑151 (MQ=255)
                                                                                                                      gAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTg               >  2:533419/1‑151 (MQ=255)
                                                                                                                        aGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGAc             <  1:152766/151‑1 (MQ=255)
                                                                                                                          ttATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGAcc            >  1:709010/1‑150 (MQ=255)
                                                                                                                           tATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGAcc            >  1:330228/1‑149 (MQ=255)
                                                                                                                                  gAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTg   >  1:502722/1‑151 (MQ=255)
                                                                                                                                   aaCCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTg   <  2:635470/150‑1 (MQ=255)
                                                                                                                                   aaCCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGt  >  1:54726/1‑151 (MQ=255)
                                                                                                                                      |                                                                                                                                                   
GGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGT  >  NC_000913/4186471‑4186752

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