Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 63,070 G→A 6.0% D65D (GAC→GAT) ‡ rapA ← RNA polymerase remodeling/recycling factor ATPase; RNA polymerase‑associated, ATP‑dependent RNA translocase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091363,0700GA 6.0% 324.1 / 6.9 102D65D (GAC→GAT) ‡rapARNA polymerase remodeling/recycling factor ATPase; RNA polymerase‑associated, ATP‑dependent RNA translocase
Reads supporting (aligned to +/- strand):  ref base G (47/48);  new base A (4/2);  total (51/51)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

Only 100 of 102 total aligned reads displayed.

ACGGTCCTGCGGTTTGCTGAACACCAGTTTGCTATCAAGGAAAACTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAAAAGTAAAGTGACAGTTCGCGCATCCACCGCGACAACGGTTCCAAGTCCCAA  >  NC_000913/62923‑63216
                                                                                                                                                   |                                                                                                                                                  
aCGGTCCTGCGGTTTGCTGAACACCAGTTTGCTATCAAGGAAAACTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGt                                                                                                                                                   >  2:1117720/1‑149 (MQ=255)
aCGGTCCTGCGGTTTGCTGAACACCAGTTTGCTATCAAGGAAAACTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGt                                                                                                                                                   <  2:17061/149‑1 (MQ=255)
     ccTGCGGTTTGCTGAACACCAGTTTGCTATCAAGGAAAACTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATgg                                                                                                                                              <  1:1413894/149‑1 (MQ=255)
      cTGCGGTTTGCTGAACACCAGTTTGCTATCAAGGAAAACTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGc                                                                                                                                             <  2:606175/149‑1 (MQ=255)
        gCGGTTTGCTGAACACCAGTTTGCTATCAAGGAAAACTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTg                                                                                                                                           >  2:422677/1‑149 (MQ=255)
            tttGCTGAACACCAGTTTGCTATCAAGGAAAACTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTaa                                                                                                                                       >  2:1362391/1‑149 (MQ=255)
            tttGCTGAACACCAGTTTGCTATCAAGGAAAACTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCATGACGGCTggcag                                                                                                                                       <  1:16555/149‑4 (MQ=255)
            gttGCTGAACACCAGTTTGCTATCAAGGAAAACTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCATGACGGCTggcag                                                                                                                                       <  1:148168/148‑4 (MQ=255)
             ttGCTGAACACCAGTTTGCTATCAAGGAAAACTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAAt                                                                                                                                      <  1:409474/149‑1 (MQ=255)
             ttGCTGAACACCAGTTTGCTATCAAGGAAAACTTCACGCAGGGCTACGAAGGAATCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTAATTCTTCGACTTGCATCTGCCAGCCTTCATGGCTGGTAAt                                                                                                                                      <  1:612271/149‑1 (MQ=255)
              tGCTGAACACCAGTTTGCTATCAAGGAAAACTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATg                                                                                                                                     <  1:1311030/149‑1 (MQ=255)
               gCTGAACACCAGTTTGCTATCAAGGAAAACTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATgg                                                                                                                                    <  2:784952/149‑1 (MQ=255)
               gCTGAACACCAGTTTGCTATCAAGGAAAACTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATgg                                                                                                                                    <  2:1301786/149‑1 (MQ=255)
                cTGAACACCAGTTTGCTATCAAGGAAAACTTCACGCTGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGt                                                                                                                                   <  1:649171/149‑1 (MQ=255)
                cTGAACACCAGTTTGCTATCAAGGAAAACTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGt                                                                                                                                   <  1:903158/149‑1 (MQ=255)
                  gAACACCAGTTTGCTATCAAGGAAAACTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTAt                                                                                                                                 <  2:715583/149‑1 (MQ=255)
                  gAACACCAGTTTGCTATCAAGGAAAACTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTATTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTAt                                                                                                                                 >  2:1127379/1‑149 (MQ=255)
                    acacCAGTTTGCTATCAAGGAAAACTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCa                                                                                                                               <  2:111743/149‑1 (MQ=255)
                      acCAGTTTGCTATCAAGGAAAACTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCAcc                                                                                                                             >  2:1530337/1‑149 (MQ=255)
                        cAGTTTGCTATCAAGGAAAACTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAg                                                                                                                           <  1:1416476/149‑1 (MQ=255)
                        cAGTTTGCTATCAAGGAAAACTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAg                                                                                                                           <  1:19794/149‑1 (MQ=255)
                         aGTTTGCTATCAAGGAAAACTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAgg                                                                                                                          <  1:109123/149‑1 (MQ=255)
                          gTTTGCTATCAAGGAAAACTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAggg                                                                                                                         <  1:1476463/149‑1 (MQ=255)
                           tttGCTATCAAGGAAAACTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGt                                                                                                                        >  2:1509636/1‑149 (MQ=255)
                               cTATCAAGGAAAACTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCGGGGTTGaa                                                                                                                    <  1:894001/149‑1 (MQ=255)
                               cTATCAAGGAAAACTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGaa                                                                                                                    <  2:940078/149‑1 (MQ=255)
                                 aTCAAGGAAAACTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACa                                                                                                                  >  1:390499/1‑149 (MQ=255)
                                  tCAAGGAAAACTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACAt                                                                                                                 >  2:581753/1‑149 (MQ=255)
                                  tCAAGGAAAACTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACAt                                                                                                                 >  1:466585/1‑149 (MQ=255)
                                   cAAGGAAAACTTCCCGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATc                                                                                                                <  1:86097/149‑1 (MQ=255)
                                   cAAGGAAAACTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATc                                                                                                                <  1:312528/149‑1 (MQ=255)
                                     aGGAAAACTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATGGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCAc                                                                                                              >  2:1045337/1‑149 (MQ=255)
                                     aGGAAAACTTAACGCAGGGCTACGCCGGACTCTTAAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGTCTGGAAAAGGTATAACCAGGGTTGAAAATCAc                                                                                                              >  1:454535/1‑149 (MQ=255)
                                      ggAAAACTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCAcg                                                                                                             >  2:61051/1‑149 (MQ=255)
                                       gAAAACTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCAcgc                                                                                                            <  1:1024627/149‑1 (MQ=255)
                                        aaaaCTTGACGCAGGGCTACGCCTGACTCTTCAGTATGCAGTCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCAcgcg                                                                                                           <  2:950609/149‑1 (MQ=255)
                                        aaaaCTTCACGCAGTGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCAagcg                                                                                                           <  2:763966/149‑1 (MQ=255)
                                           aCTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACTCGGGt                                                                                                        >  2:1422098/1‑149 (MQ=255)
                                             ttCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCa                                                                                                      >  2:626366/1‑149 (MQ=255)
                                             ttCACGCAGGGCTACGCCGGACTCGTCAGTATCCAGGAGAGTACCGATATAGGACAGCACGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCa                                                                                                      >  1:637280/1‑149 (MQ=255)
                                             ttCACGCAGGGCTAAGACGGACTCATAAGTATCCAGGCGAGTACCGATATAGGACAGCAAGCCATTTTCATCTTTTACTTCTTCGACTTGCATCAGCAAGCCGTCATGGCTGGTAATGGTATAACCAGGGTTGAACATCACCCGGGACa                                                                                                      >  1:1092632/1‑149 (MQ=255)
                                               cACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCCCg                                                                                                    >  2:10827/1‑149 (MQ=255)
                                                          aCGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACt                                                                                         <  2:466585/149‑1 (MQ=255)
                                                             ccGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACtgcg                                                                                      >  1:1336044/1‑149 (MQ=255)
                                                                   tctTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTa                                                                                >  1:1142498/1‑149 (MQ=255)
                                                                     ttCAGTATCCGGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACa                                                                              >  2:754579/1‑149 (MQ=255)
                                                                        aGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGAc                                                                           >  2:531516/1‑149 (MQ=255)
                                                                          tATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACgg                                                                         <  2:212803/149‑1 (MQ=255)
                                                                             ccAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGttt                                                                      >  2:871145/1‑149 (MQ=255)
                                                                              cAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGtttt                                                                     <  2:1579322/149‑1 (MQ=255)
                                                                               aGGCGAGTCCCGATATAGGTCAGCACGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTc                                                                    >  1:449472/1‑149 (MQ=255)
                                                                                 gCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCAc                                                                  >  2:801980/1‑149 (MQ=255)
                                                                                 gCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCAc                                                                  <  2:921504/149‑1 (MQ=255)
                                                                                   gAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCa                                                                <  2:390499/149‑1 (MQ=255)
                                                                                    aGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAg                                                               <  1:1599194/149‑1 (MQ=255)
                                                                                         cGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGAt                                                          >  2:1404040/1‑149 (MQ=255)
                                                                                          gATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATg                                                         <  1:1063784/149‑1 (MQ=255)
                                                                                           atatAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATgg                                                        <  1:27726/149‑1 (MQ=255)
                                                                                           atatAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATgg                                                        >  2:561297/1‑149 (MQ=255)
                                                                                           atatAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATgg                                                        >  2:908667/1‑149 (MQ=255)
                                                                                           atatAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATgg                                                        <  1:1362391/149‑1 (MQ=255)
                                                                                             ataGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGa                                                      >  2:1583423/1‑149 (MQ=255)
                                                                                                  tCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCCCGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAaaagt                                                 >  2:127179/1‑149 (MQ=255)
                                                                                                   cAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAaaagta                                                <  1:331797/149‑1 (MQ=255)
                                                                                                     gCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAaaagtaaa                                              <  2:529189/149‑1 (MQ=255)
                                                                                                     gCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAaaagtaaa                                              <  1:1253941/149‑1 (MQ=255)
                                                                                                     gCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAaaagtaaa                                              <  1:1117720/149‑1 (MQ=255)
                                                                                                     gCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTAGCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAaaagtaaa                                              <  2:304339/149‑1 (MQ=255)
                                                                                                      cAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAaaagtaaag                                             >  2:1409579/1‑149 (MQ=255)
                                                                                                         gCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAAAAGTAAAGTGa                                          <  2:472862/149‑1 (MQ=255)
                                                                                                          ccATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAAAAGTAAAGTGAc                                         <  2:1336044/149‑1 (MQ=255)
                                                                                                           cATTTTCTTCTTTTACTTCTTCGACTTGCAGCTTCCAGCCGTCAGGGCTGGTAATGGGATCACCAGGGGTGAACATCACGCGGGTCACGGGGGAATCACTGCGGGCGTACAGACGGTGTTCACCAGTAGATGGGAAAAGTAAAGTGACa                                        <  1:805957/149‑1 (MQ=255)
                                                                                                           aattggattcttTTACTTCTTCGACTTGCATCTGCCAGCCGTCATTGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAAAAGTAAAGTGACa                                        <  1:811484/142‑1 (MQ=255)
                                                                                                              tttCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAAAAGTAAAGTGACAGtt                                     <  2:162440/149‑1 (MQ=255)
                                                                                                                tcttctTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAAAAGTAAAGTGACAGTTcg                                   <  1:1509636/149‑1 (MQ=255)
                                                                                                                 cttcttTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAAAAGTAAAGTGACAGTTcgc                                  >  1:169967/1‑149 (MQ=255)
                                                                                                                      tttACTTATTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGCGAATCACTGCGGGCGTACAGACGGTTTTCACCAGTAGATGGGAAAAGTAAAGTGACAGTTCGCGCATc                             >  1:989568/1‑149 (MQ=255)
                                                                                                                              ttcGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAAAAGTAAAGTGACAGTTCGCGCATCCACCGCGa                     >  2:799382/1‑149 (MQ=255)
                                                                                                                               tcGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAAAAGTAAAGTGACAGTTCGCGCATCCACCGCGAc                    <  1:552610/149‑1 (MQ=255)
                                                                                                                               tcGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAAAAGTAAAGTGACAGTTCGCGCATCCACCGCGAc                    >  1:1292702/1‑149 (MQ=255)
                                                                                                                                 gACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAAAAGTAAAGTGACAGTTCGCGCATCCACCGCGACaa                  <  1:531516/149‑1 (MQ=255)
                                                                                                                                   cTTGCATCTGCCAGCCGTCATGGATGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGTAAAGTAAAGTGCCAGTTCGCGCATCCACCGCGACAACg                >  2:285756/1‑149 (MQ=255)
                                                                                                                                    ttGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAAAAGTAAAGTGACAGTTCGCGCATCCACCGCGACAACgg               >  1:1308588/1‑149 (MQ=255)
                                                                                                                                    ttGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAAAAGTAAAGTGACAGTTCGCGCATCCACCGCGACAACgg               >  1:1570622/1‑149 (MQ=255)
                                                                                                                                    ttGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAAAAGTAAAGTGACAGTTCGCGCATCCACCGCGACAACgg               >  2:574404/1‑149 (MQ=255)
                                                                                                                                     tGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATTGGAAAAGTAAAGTGACAGTTCGCGCATCCACCGCGACAACGGt              <  1:626366/149‑1 (MQ=255)
                                                                                                                                      gCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAAAAGTAAAGTGACAGTTCGCGCATCCACCGCGACAACGGtt             <  1:581753/149‑1 (MQ=255)
                                                                                                                                       cATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAAAAGTAAAGTGACAGTTCGCGCATCCACCGCGACAACGGTTc            >  1:808045/1‑149 (MQ=255)
                                                                                                                                        cattaccAGCCATGACGGCTGGCAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAAAAGTAAAGTGACAGTTCGCGCATCCACCGCGCCAACGGTTcc           >  2:1255731/6‑149 (MQ=255)
                                                                                                                                        aTCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAAAAGTAAAGTGACAGTTCGCGCATCCACCGCGACAACGGTTcc           >  1:1076854/1‑149 (MQ=255)
                                                                                                                                          ttaccAGCCATGACGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAAAAGTAAAGTGACAGTTCGCGCATCCACCGCGACAACGGTTCCaa         >  2:1178941/4‑149 (MQ=255)
                                                                                                                                          cTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAAAAGTAAAGTGACAGTTCGCGCATCCACCGCGACAACGGTTCCaa         >  2:576284/1‑149 (MQ=255)
                                                                                                                                            accAGCCATGACGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACAAGTAGATGGGAAAAGTAAAGTGACAGTTCGCGCATCCACCGCGACAACGGTTCCAAGt       >  2:1040082/2‑149 (MQ=255)
                                                                                                                                             ccAGCCGTCATGGCTGGTAATGGTATCCCCAGGGTTGAACATCACGCGTGTCACGGGGGAATCACTGCGTGCGAACAGACGGTTATCACGAGTAGATGGGAAAAGTAAAGTGACAGTTCGCGCATCCACCGCGACAACGGATCCAAGTc      >  2:1400137/1‑149 (MQ=255)
                                                                                                                                             ccAGCCATGACGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAAAAGTAAAGTGACAGTTCGCGCATCCACCGCGACAACGGTTCCAAGTc      >  2:233140/1‑149 (MQ=255)
                                                                                                                                              cAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAAAAGTAAAGTGACAGTTCGCGCATCCACCGCGACAACGGTTCCAAGTcc     <  1:871145/149‑1 (MQ=255)
                                                                                                                                               aGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAAAAGTAAAGTGACAGTTCGCGCATCCACCGCGACAACGGTTCCAAGTccc    >  2:891257/1‑149 (MQ=255)
                                                                                                                                                gCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAAAAGTAAAGTGACAGTTCGCGCATCCACCGCGACAACGGTTCCAAGTCCCa   >  2:1000740/1‑149 (MQ=255)
                                                                                                                                                 ccGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAAAAGTAAAGTGACAGTTCGCGCATCCACCGCGACAACGGTTCCAAGTCCCaa  >  2:61200/1‑149 (MQ=255)
                                                                                                                                                 ccGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAAAAGTAAAGTGACAGTTCGCGCATCCACCGCGACAACGGTTCCAAGTCCCaa  >  2:1338056/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ACGGTCCTGCGGTTTGCTGAACACCAGTTTGCTATCAAGGAAAACTTCACGCAGGGCTACGCCGGACTCTTCAGTATCCAGGCGAGTACCGATATAGGTCAGCAAGCCATTTTCTTCTTTTACTTCTTCGACTTGCATCTGCCAGCCGTCATGGCTGGTAATGGTATCACCAGGGTTGAACATCACGCGGGTCACGGGGGAATCACTGCGTGCGTACAGACGGTTTTCACCAGTAGATGGGAAAAGTAAAGTGACAGTTCGCGCATCCACCGCGACAACGGTTCCAAGTCCCAA  >  NC_000913/62923‑63216

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