Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,186,607 C→T 67.5% P420S (CCG→TCG)  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,6070CT67.5% 90.0 / 71.5 80P420S (CCG→TCG) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base C (15/11);  new base T (27/27);  total (42/38)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

TGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTG  >  NC_000913/4186464‑4186751
                                                                                                                                               |                                                                                                                                                
tGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCATGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGt                                                                                                                                              <  8:130378/148‑1 (MQ=255)
  gCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTAc                                                                                                                                            >  1:60563/1‑148 (MQ=255)
   cACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCTTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTAct                                                                                                                                           >  6:139957/1‑148 (MQ=255)
     cTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTActgc                                                                                                                                         >  5:236686/1‑148 (MQ=255)
      tGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTActgct                                                                                                                                        >  3:116349/1‑148 (MQ=255)
        gAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGa                                                                                                                                      <  8:270752/148‑1 (MQ=255)
         aGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGaa                                                                                                                                     >  3:259835/1‑148 (MQ=255)
         aGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGaa                                                                                                                                     >  4:164900/1‑148 (MQ=255)
              ttaaaaCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGTTTGAGCGCGATGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGAGAACACCCGGTCCTGCTGAACCGTg                                                                                                                                >  1:285879/4‑148 (MQ=255)
               tCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTg                                                                                                                                >  7:317747/1‑147 (MQ=255)
                  aaCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCAcc                                                                                                                            >  4:43489/1‑148 (MQ=255)
                   aCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCACCg                                                                                                                           <  3:43489/148‑1 (MQ=255)
                   aCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGt                                                                                                                                 >  8:212049/1‑142 (MQ=255)
                   aCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGt                                                                                                                                 <  7:212049/142‑1 (MQ=255)
                    ccGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCACCGa                                                                                                                          >  5:301643/1‑148 (MQ=255)
                    ccGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGa                                                                                                                          >  1:96265/1‑148 (MQ=255)
                            cTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCAc                                                                                                                  >  7:349563/1‑148 (MQ=255)
                             tACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTCGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTTCACCGACTCTGcccc                                                                                                                 >  1:296074/1‑148 (MQ=255)
                              aCGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCg                                                                                                                >  7:146065/1‑148 (MQ=255)
                               cGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCa                                                                                                                   <  4:199454/144‑1 (MQ=255)
                               cGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCa                                                                                                                   <  2:60563/144‑1 (MQ=255)
                               cGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCAcc                                                                                                                 <  7:55217/146‑1 (MQ=255)
                                      cTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTa                                                                                                        >  2:372856/1‑148 (MQ=255)
                                         gAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATc                                                                                                      >  2:195342/1‑147 (MQ=255)
                                          aaCTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATc                                                                                                      >  5:45787/1‑146 (MQ=255)
                                          aaCTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATACGCGAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATc                                                                                                      >  1:102497/1‑146 (MQ=255)
                                          aaCTGCGTGGTCATGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATCTCCTGGCCGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCa                                                                                                    <  5:304382/148‑1 (MQ=255)
                                                 tGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCAttt                                                                                             <  6:315374/148‑1 (MQ=255)
                                                    tCTTGCTACCACCATTAAAGCTGCTACGAAAATGGTGGAGAGCGAAGAAGCGGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCAACGACTCTGGAGCGTCTGGGTATCCAGGCATTTgta                                                                                          <  2:296074/148‑3 (MQ=255)
                                                    tCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGaa                                                                                          >  7:46584/1‑148 (MQ=255)
                                                      ttGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGACGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGACTGGGTATCCAGGCATATGaaac                                                                                        >  2:309117/1‑146 (MQ=255)
                                                        gCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAAc                                                                                         <  4:116349/145‑1 (MQ=255)
                                                         cTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATAGCCTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGt                                                                                     <  7:353371/148‑1 (MQ=255)
                                                           accaccATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAgg                                                                                                  <  3:121077/133‑1 (MQ=255)
                                                           accaccATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAgg                                                                                                  >  4:121077/1‑133 (MQ=255)
                                                             caccaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACt                                                                                  >  5:4893/1‑147 (MQ=255)
                                                              accaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGa                                                                                >  2:370187/1‑148 (MQ=255)
                                                              accaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGa                                                                                <  6:4893/148‑1 (MQ=255)
                                                              accaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGCCTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGa                                                                                >  1:76365/1‑148 (MQ=255)
                                                              aacaTTAAAGCTGCGAAGAAAATGGTTGAGCTCGAAGAAGCTGTCGTTTGGGATATCCTGGACTAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGa                                                                                <  4:98190/146‑1 (MQ=255)
                                                                caTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATc                                                                              >  1:46902/1‑148 (MQ=255)
                                                                   tAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCACCGAct                                                                                                                        <  7:157394/103‑1 (MQ=255)
                                                                   tAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCACCGAct                                                                                                                        >  8:157394/1‑103 (MQ=255)
                                                                    aaaGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGaa                                                                           <  2:62891/147‑1 (MQ=255)
                                                                    aaaGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAg                                                                          >  1:62891/1‑148 (MQ=255)
                                                                     aaGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAgg                                                                         <  6:45787/148‑1 (MQ=255)
                                                                           cGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTaa                                                                      <  6:301643/145‑1 (MQ=255)
                                                                           cGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTaaa                                                                     <  8:349563/146‑1 (MQ=255)
                                                                               gaaAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTa                                                                                    >  2:103307/1‑127 (MQ=255)
                                                                               gaaAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTa                                                                                    <  1:103307/127‑1 (MQ=255)
                                                                                 aaaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATc                                                               >  1:386141/1‑146 (MQ=255)
                                                                                  aaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAg                                                            >  3:269614/1‑148 (MQ=255)
                                                                                      gTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGc                                                        <  2:386141/148‑1 (MQ=255)
                                                                                              cgAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTTGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTg                                                <  8:351968/148‑1 (MQ=255)
                                                                                                 agacgCTCTCGGTTGGGATATCCTGGACGAAGCAATCCGCGAACACTCGGTACTGCCGCACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGtt                                             <  4:259835/148‑1 (MQ=255)
                                                                                                       tGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTgg                                                                                                           <  5:106974/80‑1 (MQ=255)
                                                                                                       tGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTgg                                                                                                           >  6:106974/1‑80 (MQ=255)
                                                                                                          cGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCa                                    >  1:98399/1‑148 (MQ=255)
                                                                                                            tttGGGATATCCTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCata                                  <  8:317747/148‑1 (MQ=255)
                                                                                                                ggATATCCTGGACGAAGTTAGCCGCGAACACCCGGTACTGCTGAACCGTCCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATaa                               <  6:262843/147‑1 (MQ=255)
                                                                                                                      ccTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCAtata                                >  6:372559/1‑140 (MQ=255)
                                                                                                                      ccTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCAtata                                <  5:372559/140‑1 (MQ=255)
                                                                                                                       cTGGACGAAGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACt                       >  7:240734/1‑148 (MQ=255)
                                                                                                                         ggACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGaa                                                                                          >  1:14567/1‑79 (MQ=255)
                                                                                                                         ggACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGaa                                                                                          <  2:14567/79‑1 (MQ=255)
                                                                                                                          gACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCg                    <  1:372856/148‑1 (MQ=255)
                                                                                                                            cGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGAt                  <  8:146065/148‑1 (MQ=255)
                                                                                                                              aaGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTaaa                                                                     <  1:158915/95‑1 (MQ=255)
                                                                                                                              aaGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTaaa                                                                     >  2:158915/1‑95 (MQ=255)
                                                                                                                              aaGTTATCCGCGAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATgg                <  5:139957/148‑1 (MQ=255)
                                                                                                                              aaGTTATCCGCGAACACTCGGTACTGCTGAACCGTGAACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATgg                >  8:285074/1‑148 (MQ=255)
                                                                                                                                  tATCCGCGAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGAc            >  4:211405/1‑148 (MQ=255)
                                                                                                                                      cgctAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCACCTTATGGGTATCCAGTCATTTGATCCGGTTCTGATCGAAGGTAAAGCTATCCTGCTGCACCCGCt                                                 <  6:344013/107‑1 (MQ=255)
                                                                                                                                      cgcgCACCCTCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTACAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGGGACCAg         <  6:236686/147‑1 (MQ=255)
                                                                                                                                      cgcgAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAg         <  1:370187/147‑1 (MQ=255)
                                                                                                                                       gcgAACACTCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGAt       >  7:222905/1‑148 (MQ=255)
                                                                                                                                       gcgAACACTCGGTACTGCTGAACAGTGCACCGACACTGCACCGACTGGGTATCCAGGCATTTGAACCGGTACAGATCGAAGGAAAAGCTATCCAGCTGCACCCGCTg                                                >  5:344013/1‑107 (MQ=255)
                                                                                                                                        cgAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATg      >  8:274915/1‑148 (MQ=255)
                                                                                                                                          aaCACTCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGc    <  3:284971/148‑1 (MQ=255)
                                                                                                                                            cacCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTg  <  8:46584/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
TGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTG  >  NC_000913/4186464‑4186751

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