Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,185,997 A→G T1024A (ACC→GCC)  rpoC → RNA polymerase beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,185,9970AG89.3% 755.5 / 71.8 299T1024A (ACC→GCC) rpoCRNA polymerase beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (12/20);  new base G (134/132);  total (147/152)
Fisher's exact test for biased strand distribution p-value = 1.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.68e-01

Only 100 of 302 total aligned reads displayed.

AGCTACAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGAC  >  NC_000913/4185907‑4186088
                                                                                          |                                                                                           
  cTACAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCAt                                                                                         <  1:5496914/93‑1 (MQ=255)
   tACAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATg                                                                                        <  2:3894303/93‑1 (MQ=255)
    aCAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGc                                                                                       <  2:4545140/93‑1 (MQ=255)
     cAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGcc                                                                                      <  2:3796144/93‑1 (MQ=255)
     cAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGcc                                                                                      >  2:1220/1‑93 (MQ=255)
     cAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGcc                                                                                      >  1:962027/1‑93 (MQ=255)
      aaaGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCg                                                                                     >  2:5424080/1‑93 (MQ=255)
      aaaGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCg                                                                                     >  1:5022664/1‑93 (MQ=255)
      aaaGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCg                                                                                     >  1:5346282/1‑93 (MQ=255)
       aaGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCgg                                                                                    >  2:6700148/1‑93 (MQ=255)
       aaGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCgg                                                                                    >  1:4092896/1‑93 (MQ=255)
        aGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGt                                                                                   <  2:3604734/93‑1 (MQ=255)
        aGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGt                                                                                   <  2:250101/93‑1 (MQ=255)
         gTACCTTACGGTGCGTTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGtt                                                                                  >  2:6546838/1‑93 (MQ=255)
           aCCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTAt                                                                                >  1:1826847/1‑93 (MQ=255)
           aCCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGTTAt                                                                                >  1:4563241/1‑93 (MQ=255)
             cTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCa                                                                              >  2:5729467/1‑93 (MQ=255)
                 cGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACCGGGACCCGCACGCCATGCCGGTTATCACCGa                                                                          <  2:4669930/93‑1 (MQ=255)
                   gTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAg                                                                        <  2:4296265/93‑1 (MQ=255)
                   gTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAg                                                                        >  1:4121741/1‑93 (MQ=255)
                     gCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTa                                                                      >  2:6408114/1‑93 (MQ=255)
                      cGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTaa                                                                     <  1:4194103/93‑1 (MQ=255)
                      cGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTaa                                                                     <  1:5786727/93‑1 (MQ=255)
                       ggTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAg                                                                    <  1:829409/93‑1 (MQ=255)
                        gTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGc                                                                   <  2:3346565/93‑1 (MQ=255)
                         tACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCg                                                                  >  2:304366/1‑93 (MQ=255)
                         tACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCg                                                                  >  1:6992598/1‑93 (MQ=255)
                         tACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCACCGAAGTAAGCg                                                                  >  1:715296/1‑93 (MQ=255)
                          aCTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCgg                                                                 <  1:352153/93‑1 (MQ=255)
                           cTGGCGAAAGGCGATGGCGGACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCggg                                                                >  1:7324674/1‑92 (MQ=255)
                           cTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGt                                                                >  1:4875601/1‑93 (MQ=255)
                            tGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGtt                                                               >  2:1688820/1‑93 (MQ=255)
                            tGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGtt                                                               <  1:3328722/93‑1 (MQ=255)
                             ggCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGttt                                                              >  2:3100263/1‑93 (MQ=255)
                             ggCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGttt                                                              <  1:3625852/93‑1 (MQ=255)
                             ggCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCACCGAAGTAAGCGGttt                                                              <  2:4584391/93‑1 (MQ=255)
                              gCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGtttt                                                             <  1:3144524/93‑1 (MQ=255)
                               cGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTg                                                            <  1:7482831/93‑1 (MQ=255)
                                gAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGt                                                           <  1:1436951/93‑1 (MQ=255)
                                 aaaGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTa                                                          <  2:4850671/93‑1 (MQ=255)
                                 aaaGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTa                                                          <  1:7378035/93‑1 (MQ=255)
                                 aaaGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTa                                                          >  1:5243096/1‑93 (MQ=255)
                                  aaGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTAc                                                         >  1:2044670/1‑93 (MQ=255)
                                  aaGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTAc                                                         <  2:5268687/93‑1 (MQ=255)
                                   aGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACg                                                        >  2:6010627/1‑93 (MQ=255)
                                   aGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACg                                                        >  1:611451/1‑93 (MQ=255)
                                    ggCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGc                                                       >  1:2711697/1‑93 (MQ=255)
                                     gCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCt                                                      >  1:1017444/1‑93 (MQ=255)
                                       gATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCttt                                                    >  1:1973882/1‑93 (MQ=255)
                                        aTGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTa                                                   <  2:1226096/93‑1 (MQ=255)
                                        aTGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTa                                                   <  2:2902604/93‑1 (MQ=255)
                                         tGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTAc                                                  <  1:2336078/93‑1 (MQ=255)
                                          tgCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACt                                                 <  2:2720360/92‑1 (MQ=255)
                                          ggCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACt                                                 <  1:6254015/93‑1 (MQ=255)
                                          ggCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACt                                                 >  2:3462885/1‑93 (MQ=255)
                                          ggCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACt                                                 <  2:3145696/93‑1 (MQ=255)
                                           gCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTg                                                <  1:2445349/93‑1 (MQ=255)
                                             gAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGAc                                              >  1:4674072/1‑93 (MQ=255)
                                              aaCAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACa                                             >  2:7196121/1‑93 (MQ=255)
                                              aaCAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACa                                             >  1:3059019/1‑93 (MQ=255)
                                              aaCAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACa                                             >  1:6532323/1‑93 (MQ=255)
                                               aCAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACAt                                            <  1:1616086/93‑1 (MQ=255)
                                               aCAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACAt                                            <  1:3721416/93‑1 (MQ=255)
                                                cAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATg                                           <  2:1568926/93‑1 (MQ=255)
                                                  ggTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGAt                                         <  2:5204011/93‑1 (MQ=255)
                                                  ggTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGAt                                         <  1:6799699/93‑1 (MQ=255)
                                                  ggTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGAt                                         <  1:805098/93‑1 (MQ=255)
                                                   gTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATc                                        >  2:1067425/1‑93 (MQ=255)
                                                    ttGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCg                                       <  2:6998148/93‑1 (MQ=255)
                                                      gCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGAc                                     <  2:2079782/93‑1 (MQ=255)
                                                       cTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACg                                    <  1:3376625/93‑1 (MQ=255)
                                                        tGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACgg                                   <  1:6459251/93‑1 (MQ=255)
                                                        tGGCGGCGAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACgg                                   <  1:2598227/93‑1 (MQ=255)
                                                          gcggcgAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGcc                                 <  1:5586117/93‑1 (MQ=255)
                                                          gcggcgAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGcc                                 <  1:693460/93‑1 (MQ=255)
                                                           cggcgAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCa                                >  1:120668/1‑93 (MQ=255)
                                                           cggcgAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCa                                <  1:21716/93‑1 (MQ=255)
                                                            ggcgAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAg                               <  1:443821/93‑1 (MQ=255)
                                                            ggcgAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAg                               <  2:6122512/93‑1 (MQ=255)
                                                            ggcgAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGGTATCACCGAAGTAAGCGGGTTTGTACGCTTTACTGACATGATCGACGGCCAg                               >  1:6187650/1‑93 (MQ=255)
                                                             gcgAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGa                              <  2:1865771/93‑1 (MQ=255)
                                                               gAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGAcc                            >  1:3838007/1‑93 (MQ=255)
                                                               gAAACCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGAcc                            >  2:6904272/1‑93 (MQ=255)
                                                                aaaCCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCa                           <  2:5522277/93‑1 (MQ=255)
                                                                aaaCCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCa                           <  2:7287272/93‑1 (MQ=255)
                                                                  aCCGTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCAtt                         <  2:6739374/93‑1 (MQ=255)
                                                                  aCCGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCAtt                         <  2:3562210/93‑1 (MQ=255)
                                                                  aCCGTTCCAAACTGGGACCCGCACACCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCAtt                         >  2:7096921/1‑93 (MQ=255)
                                                                   ccGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTa                        <  1:3493011/93‑1 (MQ=255)
                                                                     gTTGCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTAcg                      >  2:4072655/1‑93 (MQ=255)
                                                                        gCAAACTGGGACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGt                   >  1:6727169/1‑93 (MQ=255)
                                                                               gggACCCGCACACCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACcg            >  1:325316/1‑93 (MQ=255)
                                                                                ggACCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACcga           <  1:1567040/93‑1 (MQ=255)
                                                                                ggACCCGCACACCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACcga           <  2:6796130/93‑1 (MQ=255)
                                                                                  aCCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACcgacg         >  2:1453761/1‑93 (MQ=255)
                                                                                  aCCCGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACcgacg         <  1:1092139/93‑1 (MQ=255)
                                                                                   cccGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACcgacga        >  2:2375360/1‑93 (MQ=255)
                                                                                    ccGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGaa       >  1:1998782/1‑93 (MQ=255)
                                                                                     cGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAAc      <  2:5463636/93‑1 (MQ=255)
                                                                                     cGCACGCCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAAc      >  2:4627537/1‑93 (MQ=255)
                                                                                          |                                                                                           
AGCTACAAAGTACCTTACGGTGCGGTACTGGCGAAAGGCGATGGCGAACAGGTTGCTGGCGGCGAAACCGTTGCAAACTGGGACCCGCACACCATGCCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGAC  >  NC_000913/4185907‑4186088

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑