Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,188,572 C→A R1075S (CGT→AGT)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,188,5720CA83.6% 271.6 / 49.5 116R1075S (CGT→AGT) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base C (8/11);  new base A (50/47);  total (58/58)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

Only 100 of 116 total aligned reads displayed.

GTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGC  >  NC_000913/4188428‑4188716
                                                                                                                                                |                                                                                                                                                
gTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTcc                                                                                                                                              >  7:279917/1‑149 (MQ=255)
      aCCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACt                                                                                                                                        <  4:208872/149‑1 (MQ=255)
      aCCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACt                                                                                                                                        >  5:272049/1‑149 (MQ=255)
      aCCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACt                                                                                                                                        >  1:57398/1‑149 (MQ=255)
       ccGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTg                                                                                                                                       <  3:196369/149‑1 (MQ=255)
          aaGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGaa                                                                                                                                     >  1:292547/1‑148 (MQ=255)
           aGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGaaaa                                                                                                                                   <  6:272049/149‑1 (MQ=255)
               aGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGt                                                                                                                               <  3:102735/149‑1 (MQ=255)
                gtggTTTTTTACGGTTTTCTGAAATGATCGTCGGCCAGACCATTACGGGTCATACCGACGAACTGACCGGTCTGTCTTAGTTGTTGGTTCTGGATTTCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGATAATCGt                                                                                                                               <  4:212106/146‑1 (MQ=255)
                gCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGtt                                                                                                                              <  7:29451/149‑1 (MQ=255)
                 cGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTg                                                                                                                             <  7:196201/149‑1 (MQ=255)
                  ggTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGa                                                                                                                            <  2:106019/149‑1 (MQ=255)
                        gTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGa                                                                                                                            >  7:117097/1‑143 (MQ=255)
                        gTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGa                                                                                                                            <  8:117097/143‑1 (MQ=255)
                        gTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCt                                                                                                                        >  4:245827/1‑147 (MQ=255)
                         tACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAg                                                                                                                     >  7:266881/1‑149 (MQ=255)
                           cGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAggg                                                                                                                   >  6:131376/1‑149 (MQ=255)
                           cGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGGCCGGCACTGAAAATCGTTGATGCTCAgg                                                                                                                    <  1:105869/148‑1 (MQ=255)
                                tACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAac                                                                                                               <  3:300400/148‑1 (MQ=255)
                                     aCATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAacgacg                                                                                                           <  5:244650/147‑1 (MQ=255)
                                     aCATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGtt                                                                                                         <  1:255400/149‑1 (MQ=255)
                                      cATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTc                                                                                                        >  1:236497/1‑149 (MQ=255)
                                       aTGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCt                                                                                                       >  3:83995/1‑149 (MQ=255)
                                       aTGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCt                                                                                                       <  6:11085/149‑1 (MQ=255)
                                         gATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGa                                                                                                     >  2:299703/1‑149 (MQ=255)
                                         gATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGa                                                                                                     >  4:131615/1‑149 (MQ=255)
                                         gATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGa                                                                                                     >  3:192379/1‑149 (MQ=255)
                                           tCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGAt                                                                                                    >  8:133696/1‑148 (MQ=255)
                                           tCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGAt                                                                                                    >  2:77076/1‑148 (MQ=255)
                                           tCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGAt                                                                                                    <  1:77076/148‑1 (MQ=255)
                                              aCGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCa                                                                                                <  6:86621/149‑1 (MQ=255)
                                              aCGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCa                                                                                                <  1:218703/149‑1 (MQ=255)
                                               cGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAg                                                                                               <  5:131376/149‑1 (MQ=255)
                                                  ccAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTa                                                                                            >  6:43845/1‑149 (MQ=255)
                                                  ccAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTa                                                                                            >  2:326172/1‑149 (MQ=255)
                                                  ccAGACAATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTa                                                                                            >  3:281234/1‑149 (MQ=255)
                                                   cagaggaTTACGCGTCAGACAGACGAACTGACTGGTCTGTCTTCGCTGGTTGTTCTGGATTCCGCAGAACGTACCTCAGGTGGTAATGATCTGAGTCCTGCACTTAAATTCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTAc                                                                                           <  8:266881/143‑1 (MQ=255)
                                                   cAGACAATTACGCTTCAGAGCGACGAACTGTCCGGTCTGTCTGCGGTGGTGGTTCTGGATTCCTTAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGGAACGACGTTCTGATCCCAGGTAc                                                                                           <  6:263767/149‑1 (MQ=255)
                                                      aCCATTACGCGTCAGCCCGACGTACTGACCGGTCTGTCTTGGCTGGGGGTGCTGGATGCCGCAGAACGGACCGGAGGTGGTAAAGATCTGAGTCCGGCACTGAAACTCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGa                                                                                        <  6:91/149‑1 (MQ=255)
                                                      aCCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGa                                                                                        >  1:342000/1‑149 (MQ=255)
                                                        cATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGata                                                                                      >  4:30455/1‑149 (MQ=255)
                                                            aCGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCAATGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGc                                                                                   <  2:132555/148‑1 (MQ=255)
                                                             cgcgTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCa                                                                                 <  8:31404/149‑1 (MQ=255)
                                                             cgcgTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCa                                                                                 <  3:131615/149‑1 (MQ=255)
                                                                gTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAgcg                                                                              >  7:129469/1‑149 (MQ=255)
                                                                 tCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAgcg                                                                              >  1:140056/1‑148 (MQ=255)
                                                                   aGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTAc                                                                                           <  5:294452/133‑1 (MQ=255)
                                                                   aGACCGACGAACTGACAGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTAc                                                                                           >  6:294452/1‑133 (MQ=255)
                                                                    gACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGt                                                                          <  2:292547/149‑1 (MQ=255)
                                                                     aCCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTa                                                                         >  3:98694/1‑149 (MQ=255)
                                                                      ccGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTAc                                                                        <  4:98694/149‑1 (MQ=255)
                                                                       cgacgaACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACt                                                                       <  6:115984/149‑1 (MQ=255)
                                                                        gacgaACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACtt                                                                      >  3:13389/1‑149 (MQ=255)
                                                                             aCTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGtt                                                                                                         <  2:257087/109‑1 (MQ=255)
                                                                             aCTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTACGCAGCACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGGTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGAGCAGTACTTCCTGc                                                                 >  5:101401/1‑149 (MQ=255)
                                                                                 aCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCggg                                                             >  4:62665/1‑149 (MQ=255)
                                                                                 aCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCggg                                                             <  7:132807/149‑1 (MQ=255)
                                                                                   cGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGt                                                            >  5:277316/1‑148 (MQ=255)
                                                                                          tctTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGa                                                     <  4:281234/148‑1 (MQ=255)
                                                                                          tctTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGAt                                                    <  3:30455/149‑1 (MQ=255)
                                                                                          tctTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGGAAAGATCTGAGTCCGGCACTGAAAATAGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACAGATATGCCAGCGCAGTACTTACTGCCGGGTAAATCGAt                                                    >  2:316501/1‑149 (MQ=255)
                                                                                            ttCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTg                                                  >  1:199761/1‑149 (MQ=255)
                                                                                            ttCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGata                                                                                      <  6:286996/113‑1 (MQ=255)
                                                                                            ttCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGata                                                                                      >  5:286996/1‑113 (MQ=255)
                                                                                            ttCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTg                                                  >  3:8189/1‑149 (MQ=255)
                                                                                               gCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGtt                                                <  1:299703/148‑1 (MQ=255)
                                                                                               gCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGtt                                                <  8:279917/148‑1 (MQ=255)
                                                                                               gCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGtt                                                <  1:160554/148‑1 (MQ=255)
                                                                                               gCTGGTGGGTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGtt                                                <  8:56318/148‑1 (MQ=255)
                                                                                                   gtggtTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCt                                           >  5:10799/1‑149 (MQ=255)
                                                                                                   gtggtTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTAAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGa                                                     >  1:328187/1‑139 (MQ=255)
                                                                                                   gtggtTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTAAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGa                                                     <  2:328187/139‑1 (MQ=255)
                                                                                                      gtTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGAtt                                                   >  5:77759/1‑138 (MQ=255)
                                                                                                       ttCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTgg                                         >  6:299718/1‑147 (MQ=255)
                                                                                                          tGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGa                                                                                                                            <  8:38785/61‑1 (MQ=255)
                                                                                                          tGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGAt                                    >  3:117453/1‑149 (MQ=255)
                                                                                                           ggATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATg                                   >  2:192201/1‑149 (MQ=255)
                                                                                                              ttCCGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTAc                                                                                           >  6:325075/1‑90 (MQ=255)
                                                                                                                ccGCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGt                               <  3:156257/148‑1 (MQ=255)
                                                                                                                  gCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTa                              >  7:335258/1‑147 (MQ=255)
                                                                                                                  gCAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACa                            <  5:339289/149‑1 (MQ=255)
                                                                                                                   cAGAACGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTa                              >  1:335229/1‑146 (MQ=255)
                                                                                                                       aCGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCa                       <  1:30094/149‑1 (MQ=255)
                                                                                                                       aCGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCa                       <  4:8189/149‑1 (MQ=255)
                                                                                                                        cGTACCGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAg                      <  2:57398/149‑1 (MQ=255)
                                                                                                                         gTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGc                     >  2:166643/1‑149 (MQ=255)
                                                                                                                            ccGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCATGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGctc                   <  3:62665/148‑1 (MQ=255)
                                                                                                                             cGCAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTg                 >  3:301851/1‑149 (MQ=255)
                                                                                                                              gCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTgg                <  1:166643/149‑1 (MQ=255)
                                                                                                                               cAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGt               <  2:236497/149‑1 (MQ=255)
                                                                                                                               cAGGTGGTAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGt               <  6:10799/149‑1 (MQ=255)
                                                                                                                                     gtAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACAcc         >  1:230800/1‑149 (MQ=255)
                                                                                                                                     gtAAAGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACAcc         >  3:222265/1‑149 (MQ=255)
                                                                                                                                       aaaGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCt       >  5:134517/1‑149 (MQ=255)
                                                                                                                                        aaGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGAACAGCTCTGGTGACACCCTg      >  5:182511/1‑149 (MQ=255)
                                                                                                                                        aaGATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCt       <  3:245827/148‑1 (MQ=255)
                                                                                                                                          gATCTTAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGgc    <  6:196529/149‑1 (MQ=255)
                                                                                                                                          gATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGa                                                     <  3:93687/100‑1 (MQ=255)
                                                                                                                                          gATCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGgc    >  2:59806/1‑149 (MQ=255)
                                                                                                                                            tCTGAGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGgcgc  >  7:100345/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                
GTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGC  >  NC_000913/4188428‑4188716

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: