Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,188,5720CT49.5% 5.6 / 133.4 96R1075C (CGT→TGT) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base C (23/25);  new base T (28/19);  total (51/45)
Fisher's exact test for biased strand distribution p-value = 3.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTG  >  NC_000913/4188428‑4188712
                                                                                                                                                |                                                                                                                                            
gTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTcc                                                                                                                                          <  3:78139/149‑1 (MQ=255)
 ttATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTTACCGGTCTGTCTTCGCTTGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCg                                                                                                                                         <  1:165788/149‑1 (MQ=255)
    tCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGATTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCa                                                                                                                                      <  6:308538/149‑1 (MQ=255)
      aCCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACt                                                                                                                                    >  5:268737/1‑149 (MQ=255)
      aCCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACt                                                                                                                                    >  8:223862/1‑149 (MQ=255)
      aCCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACt                                                                                                                                    >  5:172829/1‑149 (MQ=255)
       ccGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGCACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACt                                                                                                                                    >  5:102516/1‑148 (MQ=255)
       ccGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACt                                                                                                                                    >  8:99373/1‑148 (MQ=255)
       ccGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACt                                                                                                                                    >  3:73516/1‑148 (MQ=255)
        cGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGa                                                                                                                                  <  6:268737/149‑1 (MQ=255)
          aaGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGaa                                                                                                                                 >  4:188066/1‑148 (MQ=255)
          aaGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGaa                                                                                                                                 >  1:9792/1‑148 (MQ=255)
                gCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGTCGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGGCCGGCACTGAAAATCGtt                                                                                                                          <  6:118264/149‑1 (MQ=255)
                gCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGtt                                                                                                                          <  8:158685/149‑1 (MQ=255)
                  ggTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGa                                                                                                                        <  7:223862/149‑1 (MQ=255)
                   gTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGAt                                                                                                                       <  7:308035/149‑1 (MQ=255)
                        gTACGCTTTAATGACATGATCGACGGCCAGACCATTAAGCGTCAGACCGACGAACTGACCGGTCTTGCTTCGCTGGTGGTTCTGGATGCCGCAGAACGTACCGCAGGTGGTAAAGCTCTGTGTCCGGCACTGCAAATCGTTGATGCt                                                                                                                    >  2:298888/1‑147 (MQ=255)
                         tACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGtt                                                                                                                          >  8:165395/1‑140 (MQ=255)
                         tACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGtt                                                                                                                          <  7:165395/140‑1 (MQ=255)
                         tACGCTTTACAGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAg                                                                                                                 >  6:145970/1‑149 (MQ=255)
                           cGCTTTACTGACATGATCGACGGCCAGACCATAACGCGGCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGGGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAggg                                                                                                               <  4:66500/149‑1 (MQ=255)
                                 aCTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGCTCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAacga                                                                                                         >  2:281948/1‑149 (MQ=255)
                                 aCTGACATGAACTACGGCCAGACCATTACGCGTCAGAACGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAacga                                                                                                         >  7:331364/1‑149 (MQ=255)
                                   tGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAacgacg                                                                                                       >  8:90511/1‑149 (MQ=255)
                                     aCATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCa                                                                                                                  >  5:104485/1‑136 (MQ=255)
                                     aCATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCa                                                                                                                  <  6:104485/136‑1 (MQ=255)
                                      cATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTc                                                                                                    <  6:358090/149‑1 (MQ=255)
                                       aTGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCt                                                                                                   <  3:188066/149‑1 (MQ=255)
                                         gATCGACGGCCAGTCCATTACGCGTGAGACCGACTAACTGAGCGGTCTGGCTTCGCTGGTGGTTCGGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGGGTCCGGCACTGAAAATCGTTGATGCTCAGGGTACCGACGTTCTGa                                                                                                 <  1:298888/149‑1 (MQ=255)
                                         gATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGa                                                                                                 >  4:271206/1‑149 (MQ=255)
                                         gATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGtt                                                                                                                          >  3:140171/1‑124 (MQ=255)
                                         gATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGtt                                                                                                                          <  4:140171/124‑1 (MQ=255)
                                         gATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGa                                                                                                 >  1:160112/1‑149 (MQ=255)
                                          aTCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGAt                                                                                                >  1:168751/1‑149 (MQ=255)
                                          aTCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGAt                                                                                                <  2:168751/149‑1 (MQ=255)
                                             gACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATcc                                                                                              >  1:153053/1‑148 (MQ=255)
                                              aCGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCa                                                                                            <  7:90511/149‑1 (MQ=255)
                                                  ccAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTa                                                                                        >  7:152131/1‑149 (MQ=255)
                                                  ccAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTa                                                                                        <  6:172829/149‑1 (MQ=255)
                                                              gcgTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAg                                                                            >  8:300777/1‑149 (MQ=255)
                                                                   aGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCa                                                                        <  3:271206/148‑1 (MQ=255)
                                                                       ctactcaCTGCCAGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCACAACGTAACGCAGGTGGGAAAGATCTGCTTCCGGCACAGAAAATCGGTGATGATCAGGGTAACGACGTTCTGATCCCATGTACAGATATGCCAGCGCAGTACt                                                                   >  4:208404/7‑149 (MQ=255)
                                                                             aCTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGc                                                             >  7:156184/1‑149 (MQ=255)
                                                                             aCTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTg                                                                                                                         <  4:340521/89‑1 (MQ=255)
                                                                             aCTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTg                                                                                                                         >  3:340521/1‑89 (MQ=255)
                                                                                 gtcGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCgg                                                          <  7:265661/146‑1 (MQ=255)
                                                                                 aCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCgg                                                          <  8:337755/148‑1 (MQ=255)
                                                                                 aCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCGGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCggg                                                         >  1:310673/1‑149 (MQ=255)
                                                                                   cGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACAGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGt                                                        >  3:141938/1‑148 (MQ=255)
                                                                                      tctgtctTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTcc                                                                >  7:101063/1‑137 (MQ=255)
                                                                                      tctgtctTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTcc                                                                <  8:101063/137‑1 (MQ=255)
                                                                                      tctgtctTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAg                                                    <  6:102516/149‑1 (MQ=255)
                                                                                       ctgtctTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGc                                                   <  2:9792/149‑1 (MQ=255)
                                                                                          tctTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGAt                                                >  5:71608/1‑149 (MQ=255)
                                                                                          tctTCGCTGGTGGTTATGGATTCCTCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGc                                                   <  8:364026/146‑1 (MQ=255)
                                                                                            ttCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCGGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTg                                              >  1:43871/1‑149 (MQ=255)
                                                                                              cGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGAtt                                               <  1:218673/146‑1 (MQ=255)
                                                                                              cGCTGGGGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGt                                             <  4:73516/148‑1 (MQ=255)
                                                                                               gCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGtt                                            >  5:59785/1‑148 (MQ=255)
                                                                                               gCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAACGCGATTGtt                                            >  6:275046/1‑148 (MQ=255)
                                                                                               gCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGtt                                            >  7:84289/1‑148 (MQ=255)
                                                                                                 tggtggTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATc                                                                                                                             <  1:123436/65‑1 (MQ=255)
                                                                                                 tggtggTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATc                                                                                                                             >  2:123436/1‑65 (MQ=255)
                                                                                                  ggtggtTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGc                                        >  3:172712/1‑149 (MQ=255)
                                                                                                       ttCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGaa                                   <  6:59785/149‑1 (MQ=255)
                                                                                                            gATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAac                                                                                                           <  6:111453/72‑1 (MQ=255)
                                                                                                            gATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAac                                                                                                           >  5:111453/1‑72 (MQ=255)
                                                                                                            gATTCCGCAGAACGTACCGCAGGTGGTAAAGATCCGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGAt                                <  2:310673/147‑1 (MQ=255)
                                                                                                            gATTCCACAGATCGTACAGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCTGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATgg                              >  3:101305/1‑149 (MQ=255)
                                                                                                              ttCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGaa                                   <  8:322209/142‑1 (MQ=255)
                                                                                                              ttCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGaa                                   >  7:322209/1‑142 (MQ=255)
                                                                                                              ttCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCg                            <  2:153053/149‑1 (MQ=255)
                                                                                                                 cGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTAc                         <  2:160112/149‑1 (MQ=255)
                                                                                                                  gCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTa                          >  5:59142/1‑147 (MQ=255)
                                                                                                                  gCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACa                        <  8:152131/149‑1 (MQ=255)
                                                                                                                  gCAGAACGTACCGCAGGTGGTAAAGATATGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACa                        >  5:290953/1‑149 (MQ=255)
                                                                                                                    aGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATTCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGa                      <  4:15380/149‑1 (MQ=255)
                                                                                                                          caCCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGct                >  7:157151/2‑149 (MQ=255)
                                                                                                                            ccGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGctct              >  5:124580/1‑149 (MQ=255)
                                                                                                                            ccGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGctc               <  6:71608/148‑1 (MQ=255)
                                                                                                                             cGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGa                                 <  5:246746/129‑1 (MQ=255)
                                                                                                                             cGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTg             <  7:300777/149‑1 (MQ=255)
                                                                                                                             cGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCATGCCGGGTAAAGCGATTGTTCAGCTGGAAGa                                 >  6:246746/1‑129 (MQ=255)
                                                                                                                               catgtggtAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCTGGTACCGATATGACAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGt           <  7:337624/146‑1 (MQ=255)
                                                                                                                               cAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATTCCAGCGCAGTACTTCCTGCAGGGTAAAGCGATTGATCAGCTGGAAGATGGCGTACAGATCAGATCTgct           >  5:119334/1‑147 (MQ=255)
                                                                                                                               cAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTAAAGATCAGCTCTGGt           <  7:99373/149‑1 (MQ=255)
                                                                                                                               cAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGt           <  6:124580/149‑1 (MQ=255)
                                                                                                                               cAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGt           <  4:172712/149‑1 (MQ=255)
                                                                                                                                    ggtAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGacac      >  1:150802/1‑149 (MQ=255)
                                                                                                                                     gtAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACAcc     >  3:267110/1‑149 (MQ=255)
                                                                                                                                     gtAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACAcc     >  3:230105/1‑149 (MQ=255)
                                                                                                                                       aaaGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCt   >  6:299563/1‑149 (MQ=255)
                                                                                                                                       aaaGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTGCCTGCCGGGAAAAGCGATTGTTCCGCTCGAAGATGGCGTACAGATCAGCTCTGGTGACACCCt   >  3:309136/1‑149 (MQ=255)
                                                                                                                                        aaGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTg  >  8:184080/1‑149 (MQ=255)
                                                                                                                                           aTCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCg                                                           >  4:349254/1‑89 (MQ=255)
                                                                                                                                           aTCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCg                                                           <  3:349254/89‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                            
GTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTG  >  NC_000913/4188428‑4188712

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