Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,188,572 C→T 19.2% R1075C (CGT→TGT) ‡ 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,5720CT19.2% 180.7 / 35.0 84R1075C (CGT→TGT) ‡rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base C (35/32);  new base T (8/8);  total (44/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.64e-01

CCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCG  >  NC_000913/4188425‑4188717
                                                                                                                                                   |                                                                                                                                                 
ccGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGct                                                                                                                                                  <  5:74009/149‑2 (MQ=255)
ccGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGct                                                                                                                                                  <  1:385329/149‑2 (MQ=255)
  ggTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTc                                                                                                                                                <  3:193386/149‑1 (MQ=255)
   gTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTcc                                                                                                                                               >  6:208646/1‑149 (MQ=255)
       tCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCa                                                                                                                                           >  8:221210/1‑149 (MQ=255)
        cACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACAGACGAACTGACCGGTCTGGCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGAAGGTGGTAAAGATCTGATTCCGGCAc                                                                                                                                          >  5:168955/1‑149 (MQ=255)
         aCCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACt                                                                                                                                         >  4:199109/1‑149 (MQ=255)
               gTAAGCGGGTTTGGACGCTTTACTGACATTAGAGACGGCCAGACCATTACGCGTCAGACCGACGTACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGAGCTGCTTCCGGCACTGAAAAt                                                                                                                                   <  8:434994/149‑1 (MQ=255)
                  aGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGt                                                                                                                                >  3:253346/1‑149 (MQ=255)
                   gCGTTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGtg                                                                                                                               <  4:326331/149‑2 (MQ=255)
                   gCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGtt                                                                                                                               <  6:163697/149‑1 (MQ=255)
                      gTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGAt                                                                                                                            <  3:199109/149‑1 (MQ=255)
                            tACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAg                                                                                                                      >  8:324418/1‑149 (MQ=255)
                            tACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAg                                                                                                                      >  7:146470/1‑149 (MQ=255)
                              cGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAggg                                                                                                                    >  2:355520/1‑149 (MQ=255)
                              cGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAgg                                                                                                                     <  6:143495/148‑1 (MQ=255)
                              cGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAgg                                                                                                                     >  1:350704/1‑148 (MQ=255)
                              cGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAgg                                                                                                                     <  2:350704/148‑1 (MQ=255)
                                   tACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTaa                                                                                                                 <  2:217379/147‑1 (MQ=255)
                                      tGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAacaacg                                                                                                            >  7:89279/1‑149 (MQ=255)
                                       gACATGATCGACGGCCAGACCATTACGTGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGAGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGt                                                                                                           >  2:435351/1‑149 (MQ=255)
                                        aCATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGtt                                                                                                          >  2:119568/1‑149 (MQ=255)
                                          aTGATCGACGGCCAGACCATTACGCGTCAGCCCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCt                                                                                                        >  7:183926/1‑149 (MQ=255)
                                          aTGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCt                                                                                                        <  7:354739/149‑1 (MQ=255)
                                           tGATCGACTGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTg                                                                                                       <  5:326716/149‑1 (MQ=255)
                                                  cGGCCAGACCCTTACGCGGCAGACCGACGAACTGACCGGGCTGGCTTCGGGGGTGGGTCTGGATTCCGCAGAACGTACCGGCGGGGGGAAAGCTCTGTGGCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAg                                                                                                <  1:283456/149‑1 (MQ=255)
                                                     ccAGACCATTACCCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTa                                                                                             >  7:363628/1‑149 (MQ=255)
                                                      cAGACCATTACGCGTCAGACCGACGAACTGCCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACGGAAAATCGTTGATGCTCAGGGTAACGACGGTCTGATCCCCGGTAc                                                                                            >  5:340053/1‑149 (MQ=255)
                                                           cATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGata                                                                                       >  5:136473/1‑149 (MQ=255)
                                                             ttACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCTGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATg                                                                                     <  2:66386/149‑1 (MQ=255)
                                                             ttACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATg                                                                                     >  4:204506/1‑149 (MQ=255)
                                                              tACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGtt                                                                                                                               >  6:6672/1‑106 (MQ=255)
                                                              tACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGtt                                                                                                                               <  5:6672/106‑1 (MQ=255)
                                                               aCGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGc                                                                                    >  8:144760/1‑148 (MQ=255)
                                                               aCGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGc                                                                                    <  1:355520/148‑1 (MQ=255)
                                                                cgcgTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGc                                                                                    >  4:47890/1‑147 (MQ=255)
                                                                cgcgTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCa                                                                                  <  7:144760/149‑1 (MQ=255)
                                                                cgcgTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACAACGTTCTGATCCCAGGTACCGATATGCCa                                                                                  <  8:89279/149‑1 (MQ=255)
                                                                       gACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGt                                                                           <  8:58725/149‑1 (MQ=255)
                                                                        aCCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTa                                                                          >  4:127191/1‑149 (MQ=255)
                                                                          cgacgaACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACt                                                                        >  7:231601/1‑149 (MQ=255)
                                                                           gacgaACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACtt                                                                       >  4:89865/1‑149 (MQ=255)
                                                                              gaACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCt                                                                    <  3:127191/149‑1 (MQ=255)
                                                                                    aCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCgg                                                               <  5:208646/148‑1 (MQ=255)
                                                                                    aCCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCgg                                                               <  3:47890/148‑1 (MQ=255)
                                                                                      cGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGt                                                             >  7:410401/1‑148 (MQ=255)
                                                                                      cGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGt                                                             >  8:424172/1‑148 (MQ=255)
                                                                                           tgtctTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGc                                                        >  3:79882/1‑148 (MQ=255)
                                                                                           tgtctTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGc                                                        >  4:220688/1‑148 (MQ=255)
                                                                                            gtctTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGa                                                      >  4:117934/1‑149 (MQ=255)
                                                                                                    tggtggTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAgg                                                                                                                     <  5:421429/78‑1 (MQ=255)
                                                                                                    tggtggTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAgg                                                                                                                     >  6:421429/1‑78 (MQ=255)
                                                                                                    tggtggTTCCGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGAGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAg                                              >  7:59494/1‑149 (MQ=255)
                                                                                                     ggtggtTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAgg                                                                                               <  5:47467/99‑1 (MQ=255)
                                                                                                     ggtggtTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAgg                                                                                               >  6:47467/1‑99 (MQ=255)
                                                                                                      gtggtTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCt                                            >  2:268334/1‑149 (MQ=255)
                                                                                                      gtggtTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGc                                             <  4:253346/148‑1 (MQ=255)
                                                                                                       tggtTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTg                                           >  1:400886/1‑149 (MQ=255)
                                                                                                       tggtTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTg                                           >  6:409142/1‑149 (MQ=255)
                                                                                                         gtTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATACCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGa                                         <  7:221210/149‑1 (MQ=255)
                                                                                                              ggATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATg                                    >  7:192207/1‑149 (MQ=255)
                                                                                                               gATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATgg                                   >  3:265277/1‑149 (MQ=255)
                                                                                                               gATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAg                                                                                                <  8:56803/88‑1 (MQ=255)
                                                                                                               gATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAg                                                                                                >  7:56803/1‑88 (MQ=255)
                                                                                                               gATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATg                                    <  5:188303/148‑1 (MQ=255)
                                                                                                               gATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATg                                    <  2:400886/148‑1 (MQ=255)
                                                                                                                 ttCCGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCg                                 <  8:192207/149‑1 (MQ=255)
                                                                                                                    cGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGc                                             <  7:195634/134‑1 (MQ=255)
                                                                                                                    cGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTAc                              <  4:79882/149‑1 (MQ=255)
                                                                                                                    cGCAGAACGTACCGCAGGTGGTAAAGATCTGTGTCAGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGc                                             >  8:195634/1‑134 (MQ=255)
                                                                                                                    cGCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTa                               <  3:89865/148‑1 (MQ=255)
                                                                                                                     gCAGAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTa                               >  4:138971/1‑147 (MQ=255)
                                                                                                                        gAACGTACCGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGAt                          <  7:424172/149‑1 (MQ=255)
                                                                                                                         aaCGTACCGCAGGTGGTAAAGATCTGTGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATc                         <  6:136473/149‑1 (MQ=255)
                                                                                                                                cGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTATAGCGATTGTTCTGCTGGAAGATGGCGTACAGATCAGCTCTg                  <  8:146470/149‑1 (MQ=255)
                                                                                                                                cGCAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTg                  <  5:409142/149‑1 (MQ=255)
                                                                                                                                  cAGGTGGTAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGt                <  8:231601/149‑1 (MQ=255)
                                                                                                                                         tAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGTAAGATGGCGTACAGATCAGCTCTGGTGACAccc         <  8:59494/149‑1 (MQ=255)
                                                                                                                                         tAAAGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACAccc         >  7:378646/1‑149 (MQ=255)
                                                                                                                                           aaGATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCt        <  3:117934/148‑1 (MQ=255)
                                                                                                                                             gATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGgc     >  4:430850/1‑149 (MQ=255)
                                                                                                                                             gATCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGgc     <  7:324418/149‑1 (MQ=255)
                                                                                                                                               tCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGtt                                                 <  8:420926/103‑1 (MQ=255)
                                                                                                                                               tCTGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGtt                                                 >  7:420926/1‑103 (MQ=255)
                                                                                                                                                 tGCTTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGcgcg  >  4:89659/1‑148 (MQ=255)
                                                                                                                                                 tGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTGAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGcgcg  >  4:321595/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCGGTTATCACCGAAGTAAGCGGTTTTGTACGCTTTACTGACATGATCGACGGCCAGACCATTACGCGTCAGACCGACGAACTGACCGGTCTGTCTTCGCTGGTGGTTCTGGATTCCGCAGAACGTACCGCAGGTGGTAAAGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCG  >  NC_000913/4188425‑4188717

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: ‑