Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,186,578 A→G 65.8% D410G (GAT→GGT)  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,5780AG65.8% 280.3 / 223.5 225D410G (GAT→GGT) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (37/40);  new base G (77/71);  total (114/111)
Fisher's exact test for biased strand distribution p-value = 5.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 227 total aligned reads displayed.

GCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACG  >  NC_000913/4186432‑4186724
                                                                                                                                                  |                                                                                                                                                  
gCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGgtat                                                                                                                                                 >  1:891407/1‑150 (MQ=255)
        aTCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACg                                                                                                                                         >  1:2701614/1‑150 (MQ=255)
         tCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGa                                                                                                                                        <  1:947160/150‑1 (MQ=255)
          cAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGaa                                                                                                                                       <  1:274497/150‑1 (MQ=255)
           aGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAg                                                                                                                                      <  1:1013155/150‑1 (MQ=255)
           aGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAg                                                                                                                                      <  2:2195066/150‑1 (MQ=255)
             tGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGtt                                                                                                                                    >  2:1647426/1‑150 (MQ=255)
             tGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGtt                                                                                                                                    >  2:2964231/1‑150 (MQ=255)
                    tGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAATTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCcgcg                                                                                                                             >  1:43947/1‑150 (MQ=255)
                    tGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCcgcg                                                                                                                             <  1:848461/150‑1 (MQ=255)
                    tGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCcgcg                                                                                                                             <  1:2406763/150‑1 (MQ=255)
                    tGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCcgcg                                                                                                                             >  2:567513/1‑150 (MQ=255)
                      ccGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGaa                                                                                                                           >  1:317650/1‑150 (MQ=255)
                      ccGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGaa                                                                                                                           >  1:68488/1‑150 (MQ=255)
                        gaagaaAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAaca                                                                                                                         <  2:2645650/150‑1 (MQ=255)
                        gaagaaAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAaca                                                                                                                         >  2:1040460/1‑150 (MQ=255)
                        gaagaaAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAaca                                                                                                                         <  1:2957130/150‑1 (MQ=255)
                         aagaaAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAacac                                                                                                                        <  1:2390722/150‑1 (MQ=255)
                          agaaAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACAcc                                                                                                                       <  2:2908415/150‑1 (MQ=255)
                           gaaAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACAccc                                                                                                                      <  2:205544/150‑1 (MQ=255)
                            aaaaTGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCg                                                                                                                     >  1:3054668/1‑150 (MQ=255)
                             aaaTGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCgg                                                                                                                    <  1:2249421/150‑1 (MQ=255)
                             aaaTGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCgg                                                                                                                    <  2:329935/150‑1 (MQ=255)
                             aaaTGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCgg                                                                                                                    <  2:1494587/150‑1 (MQ=255)
                                tGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTAc                                                                                                                 >  2:2089212/1‑150 (MQ=255)
                                    aCTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTActgct                                                                                                             <  2:2775888/150‑1 (MQ=255)
                                     cTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTActgctg                                                                                                            >  2:2416749/1‑150 (MQ=255)
                                     cTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTActgctg                                                                                                            >  2:1291775/1‑150 (MQ=255)
                                          gCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCg                                                                                                       >  1:480670/1‑150 (MQ=255)
                                          gCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCg                                                                                                       <  1:1500224/150‑1 (MQ=255)
                                           cTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGt                                                                                                      <  1:2512279/150‑1 (MQ=255)
                                           cTGTTCAAACCGGTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGt                                                                                                      <  1:874635/150‑1 (MQ=255)
                                            tGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTg                                                                                                     <  1:1322648/150‑1 (MQ=255)
                                               tCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCAc                                                                                                  <  2:230336/150‑1 (MQ=255)
                                                 aaaCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCg                                                                                                >  2:2968938/1‑150 (MQ=255)
                                                 aaaCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCg                                                                                                >  1:3019110/1‑150 (MQ=255)
                                                 aaaCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCg                                                                                                >  1:2183017/1‑150 (MQ=255)
                                                 aaaCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCg                                                                                                >  2:1228563/1‑150 (MQ=255)
                                                  aaCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGa                                                                                               <  2:807157/150‑1 (MQ=255)
                                                  aaCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGa                                                                                               >  1:2162002/1‑150 (MQ=255)
                                                           tCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCAcc                                                                                      <  2:865550/150‑1 (MQ=255)
                                                             tACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGt                                                                                    >  1:2137781/1‑150 (MQ=255)
                                                             tACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGt                                                                                    >  1:1592878/1‑150 (MQ=255)
                                                               cGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCt                                                                                  >  2:3097391/1‑150 (MQ=255)
                                                                ggCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTg                                                                                 <  1:1039157/150‑1 (MQ=255)
                                                                ggCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTg                                                                                 >  2:930921/1‑150 (MQ=255)
                                                                 gCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTgg                                                                                <  1:2483877/150‑1 (MQ=255)
                                                                  cAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGTAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTggg                                                                               >  2:2428258/1‑150 (MQ=255)
                                                                   aaGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGt                                                                              >  2:2543097/1‑150 (MQ=255)
                                                                     gCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTAt                                                                            >  1:2839812/1‑150 (MQ=255)
                                                                         gAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAg                                                                        <  2:2318689/150‑1 (MQ=255)
                                                                           aCTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGc                                                                      <  2:2143549/150‑1 (MQ=255)
                                                                           aCTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGc                                                                      >  2:1834467/1‑150 (MQ=255)
                                                                             tGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCAt                                                                    >  1:400172/1‑150 (MQ=255)
                                                                             tGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCAt                                                                    >  1:1590272/1‑150 (MQ=255)
                                                                               cGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCAttt                                                                  <  1:51248/150‑1 (MQ=255)
                                                                                   gTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAAc                                                              <  1:610600/150‑1 (MQ=255)
                                                                                   gTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAAc                                                              <  2:60387/150‑1 (MQ=255)
                                                                                   gTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGCTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAAc                                                              <  1:610695/150‑1 (MQ=255)
                                                                                     ctcgcgaccaccATTAAAGCTGCGAAGACAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCGGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCg                                                            <  1:499196/144‑1 (MQ=255)
                                                                                      ttGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCgg                                                           <  1:2769359/150‑1 (MQ=255)
                                                                                      ttGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCgg                                                           >  1:748085/1‑150 (MQ=255)
                                                                                            ccaccaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGa                                                     >  2:1717997/1‑150 (MQ=255)
                                                                                            ccaccaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGa                                                     >  2:1004241/1‑150 (MQ=255)
                                                                                             caccaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGAt                                                    >  1:1960191/1‑150 (MQ=255)
                                                                                              accaTTAAAGCTGCGAAGAAAATGGTTGTGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATACGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATc                                                   >  2:1282305/1‑150 (MQ=255)
                                                                                              accaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATc                                                   <  2:1373394/150‑1 (MQ=255)
                                                                                                 aTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGaa                                                <  2:2059560/150‑1 (MQ=255)
                                                                                                  tttaaGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTTGACGAAGTTATCCGCGAACCCCCGGTACTTCTGAACCGTGCACCGACTCTGCACCGTCTGGTTATCCAGGCATTTGAACCGGTACTGATCTAAg                                               <  1:2894908/147‑1 (MQ=255)
                                                                                                  ttAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAg                                               <  1:796887/150‑1 (MQ=255)
                                                                                                   tAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAgg                                              >  2:2761957/1‑150 (MQ=255)
                                                                                                     aaGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTa                                            >  1:2299614/1‑150 (MQ=255)
                                                                                                     aaGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCAGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTa                                            >  1:2323610/1‑150 (MQ=255)
                                                                                                      aGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACACTGCACCGTCTGGGTATACAGGCAGTTGAACCGGTACTGATCGAAGGTaa                                           >  1:1476413/1‑150 (MQ=255)
                                                                                                      aGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGGTATCCGCGAACACCCGGTACAGCTGAACCCGGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGACCCGGTACGGATCGAAGGGaa                                           >  1:1702175/1‑150 (MQ=255)
                                                                                                       gCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTaaa                                          <  2:317650/150‑1 (MQ=255)
                                                                                                       gCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTaaa                                          <  2:2212409/150‑1 (MQ=255)
                                                                                                       gCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTaaa                                          <  2:2045009/150‑1 (MQ=255)
                                                                                                        cTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCCGTACAGCTGAACCGTGCACCCACGCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAg                                         >  1:1108238/1‑150 (MQ=255)
                                                                                                          gCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATAAACTGGACGAAGTTATACGCGAACACCCCGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCt                                       >  1:1400060/1‑150 (MQ=255)
                                                                                                           cGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTa                                      >  2:2597671/1‑150 (MQ=255)
                                                                                                               gaaAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATACGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATACAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTaaaca                                  >  1:1721367/1‑146 (MQ=255)
                                                                                                                aaaaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAg                                 <  2:68488/150‑1 (MQ=255)
                                                                                                                aaaaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAg                                 <  1:2616319/150‑1 (MQ=255)
                                                                                                                 aaaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGc                                <  1:257444/150‑1 (MQ=255)
                                                                                                                  aaTGGTTGAGCGCGAAGAAGCTGTCGTTTTTGGTATCCTTGAAGAAGTTATAACCGAAAAAACCGTAATGCTGAACCCTGCACCGACACGGCACCCTCTTGGTATACACGCAATTGAAACCGTACTGATAGAAGGTAAATCTaaaaagca                               >  1:3061229/1‑143 (MQ=255)
                                                                                                                  aaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCt                               <  2:396140/150‑1 (MQ=255)
                                                                                                                     ggTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAgcagca                            >  2:1547089/1‑150 (MQ=255)
                                                                                                                     ggTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCa                            >  2:271068/1‑150 (MQ=255)
                                                                                                                      gTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCCGTACAGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATACAGGCAATTGAACCGGTACTGATAGAAGGTAAAGCTATCCAGCTGCAc                           >  1:175899/1‑150 (MQ=255)
                                                                                                                      gTTGAGCGCGAAGAAGCTGTCGTTTGGGGTATCCAGGACGAAGTTATACGCGAACACCCCGTACTGCATAAACCTGCACCGACACATCAACCTCAGGGTATACAGGCATTTGAAACGGTACAGATAGAAGGTAAAGCTATACAGcatcac                           >  1:3076768/1‑145 (MQ=255)
                                                                                                                            cgcgAAGAAGCTGTCGTTTGGGGGATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTg                     <  2:1424267/150‑1 (MQ=255)
                                                                                                                                    agCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGc             <  1:1870546/150‑1 (MQ=255)
                                                                                                                                    agCTGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGc             <  1:1805778/150‑1 (MQ=255)
                                                                                                                                    agCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGc             <  2:1761383/150‑1 (MQ=255)
                                                                                                                                      cTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCgg           <  2:3063725/150‑1 (MQ=255)
                                                                                                                                       tGTCGTTTGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGc          <  2:2394075/150‑1 (MQ=255)
                                                                                                                                            tttGGGGTATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCAtata     >  1:2687700/1‑150 (MQ=255)
                                                                                                                                             ttGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATaa    >  1:368811/1‑150 (MQ=255)
                                                                                                                                               gggATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACg  >  2:2128037/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACG  >  NC_000913/4186432‑4186724

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