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,6050AC98.8% 303.5 / ‑3.5 84H419P (CAC→CCC) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (45/38);  total (45/39)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

GAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTG  >  NC_000913/4186459‑4186751
                                                                                                                                                  |                                                                                                                                                  
gaaAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCg                                                                                                                                                <  2:630579/151‑1 (MQ=255)
gaaAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCg                                                                                                                                                >  1:579542/1‑151 (MQ=255)
 aaaaTGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCgg                                                                                                                                               >  2:150516/1‑151 (MQ=255)
        cACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTActgct                                                                                                                                        <  2:656314/151‑1 (MQ=255)
        cACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTActgct                                                                                                                                        >  1:940834/1‑151 (MQ=255)
         aCTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTActgctg                                                                                                                                       >  2:534227/1‑151 (MQ=255)
         aCTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATACCCTGGACGAAGTTATCCGCGAACCCCCGGTActgctg                                                                                                                                       >  1:793190/1‑151 (MQ=255)
           tGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGa                                                                                                                                      >  1:554214/1‑150 (MQ=255)
             gAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAAcc                                                                                                                                   <  1:1011445/151‑1 (MQ=255)
              aGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCg                                                                                                                                  <  1:857884/151‑1 (MQ=255)
                  gTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTg                                                                                                                                <  1:571279/149‑1 (MQ=255)
                  gTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGc                                                                                                                               >  2:229830/1‑150 (MQ=255)
                  gTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGc                                                                                                                               <  1:705562/150‑1 (MQ=255)
                      aaaCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCg                                                                                                                           >  1:1019244/1‑150 (MQ=255)
                        aCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGAct                                                                                                                        >  1:396466/1‑151 (MQ=255)
                        aCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGAct                                                                                                                        <  1:310471/151‑1 (MQ=255)
                          cGTTCATCTACGGCAAGCTGGAACTTCGTGGTCTTACAACCACCATTAAAGCTGCGAAGAAAATGGCTGCGCGCGAAGAAGCTGTCGGATTGGATATCCTGGACGCAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGAc                                                                                                                         >  2:600645/1‑148 (MQ=255)
                           gTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTg                                                                                                                     <  2:898423/151‑1 (MQ=255)
                           gTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTg                                                                                                                     >  1:332467/1‑151 (MQ=255)
                                   acagCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTc                                                                                                              <  2:793405/147‑1 (MQ=255)
                                   aCGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCt                                                                                                             <  2:208533/151‑1 (MQ=255)
                                    cGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTg                                                                                                            >  2:906739/1‑151 (MQ=255)
                                       cAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGt                                                                                                         >  1:493991/1‑151 (MQ=255)
                                          gCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATc                                                                                                      >  1:275252/1‑151 (MQ=255)
                                            tGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATcc                                                                                                     >  1:89702/1‑150 (MQ=255)
                                               aaCTGCGTGGTCTTGCTACCACCAATAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGc                                                                                                 >  2:508167/1‑151 (MQ=255)
                                                 cTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCAt                                                                                               >  1:204782/1‑151 (MQ=255)
                                                  tGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCAtt                                                                                              >  2:660316/1‑151 (MQ=255)
                                                   gCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCAttt                                                                                             <  1:626298/151‑1 (MQ=255)
                                                   gCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCAttt                                                                                             <  1:126486/151‑1 (MQ=255)
                                                   gCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCAttt                                                                                             >  1:886800/1‑151 (MQ=255)
                                                   gCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGGATCCAGGCAttt                                                                                             >  2:760250/1‑151 (MQ=255)
                                                    cGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTg                                                                                            <  1:722689/151‑1 (MQ=255)
                                                     gTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTg                                                                                            >  1:751480/1‑150 (MQ=255)
                                                         tCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCg                                                                                       >  1:821000/1‑151 (MQ=255)
                                                            tGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCgg                                                                                      <  2:809239/149‑1 (MQ=255)
                                                               tACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTg                                                                                 >  1:824293/1‑151 (MQ=255)
                                                               tACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTg                                                                                 >  1:200457/1‑151 (MQ=255)
                                                                accgccATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACt                                                                                  <  2:357343/149‑1 (MQ=255)
                                                                accaccATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTg                                                                                 >  2:942166/1‑150 (MQ=255)
                                                                accaccATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTg                                                                                 <  1:438582/150‑1 (MQ=255)
                                                                  caccaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATc                                                                              >  1:938599/1‑151 (MQ=255)
                                                                  caccaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATc                                                                              <  2:540959/151‑1 (MQ=255)
                                                                     caTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGa                                                                            >  2:1049954/1‑150 (MQ=255)
                                                                      aTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAg                                                                          <  1:472403/151‑1 (MQ=255)
                                                                        tAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAgg                                                                         <  2:964515/150‑1 (MQ=255)
                                                                         aaaGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGt                                                                        >  1:840086/1‑150 (MQ=255)
                                                                             cTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGc                                                                   >  1:760846/1‑151 (MQ=255)
                                                                             cTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGc                                                                   >  2:355648/1‑151 (MQ=255)
                                                                                 gaagaaAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATc                                                               >  2:81091/1‑151 (MQ=255)
                                                                                     aaaaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGc                                                           >  2:1054451/1‑151 (MQ=255)
                                                                                     aaaaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGc                                                           <  2:820441/151‑1 (MQ=255)
                                                                                        aTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGc                                                        <  1:983948/151‑1 (MQ=255)
                                                                                           gTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGCACCCACGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCa                                                       >  2:1062708/1‑149 (MQ=255)
                                                                                             tGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCAccc                                                    <  1:531953/150‑1 (MQ=255)
                                                                                             tGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCAccc                                                    <  1:556136/150‑1 (MQ=255)
                                                                                              gAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGACCCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCg                                                   <  2:965231/150‑1 (MQ=255)
                                                                                                 cgcgAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTgg                                               <  2:383004/151‑1 (MQ=255)
                                                                                                  gcgAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGt                                              <  1:544667/151‑1 (MQ=255)
                                                                                                   cgAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGtt                                             <  2:396257/151‑1 (MQ=255)
                                                                                                    gaagaaGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGttt                                            >  1:498753/1‑151 (MQ=255)
                                                                                                    gaagaaGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGttt                                            <  2:56068/151‑1 (MQ=255)
                                                                                                            tGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGc                                     >  2:119233/1‑150 (MQ=255)
                                                                                                            tGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGc                                     >  2:739601/1‑150 (MQ=255)
                                                                                                            tGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGc                                     >  1:922044/1‑150 (MQ=255)
                                                                                                                gTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCatat                                 >  1:76043/1‑150 (MQ=255)
                                                                                                                    gggATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACg                             <  1:504843/150‑1 (MQ=255)
                                                                                                                    gggATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACg                             <  2:222060/150‑1 (MQ=255)
                                                                                                                    gggATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACg                             <  1:150682/150‑1 (MQ=255)
                                                                                                                    gggATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGc                            >  1:488748/1‑151 (MQ=255)
                                                                                                                      gATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCg                          <  2:442773/151‑1 (MQ=255)
                                                                                                                           ccTGGACGACGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTTAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTc                     >  2:884876/1‑151 (MQ=255)
                                                                                                                           ccTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTc                     <  1:760544/151‑1 (MQ=255)
                                                                                                                             tGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCg                    >  1:91061/1‑150 (MQ=255)
                                                                                                                               gACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCGCCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATg                 <  2:333967/151‑1 (MQ=255)
                                                                                                                               gACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATg                 <  1:421078/151‑1 (MQ=255)
                                                                                                                               gACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATg                 <  2:940495/151‑1 (MQ=255)
                                                                                                                                  gAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTg              >  2:398245/1‑151 (MQ=255)
                                                                                                                                    aGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGAc            >  1:995530/1‑151 (MQ=255)
                                                                                                                                       tATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGAcc           >  2:960386/1‑149 (MQ=255)
                                                                                                                                          ccGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGa        <  2:275074/149‑1 (MQ=255)
                                                                                                                                          ccGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATg      <  2:852987/151‑1 (MQ=255)
                                                                                                                                          ccGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATg      <  2:430111/151‑1 (MQ=255)
                                                                                                                                              gAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTg  >  1:155175/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTG  >  NC_000913/4186459‑4186751

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