Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,186,577 G→C D410H (GAT→CAT)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,186,5770GC98.8% 305.3 / NA 83D410H (GAT→CAT) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (40/42);  minor base A (0/1);  total (40/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATA  >  NC_000913/4186436‑4186721
                                                                                                                                             |                                                                                                                                                
cTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCt                                                                                                                                           <  6:274202/149‑1 (MQ=255)
  gCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTg                                                                                                                                          >  3:49031/1‑148 (MQ=255)
  gCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTg                                                                                                                                          >  8:39236/1‑148 (MQ=255)
   cATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTgg                                                                                                                                         <  7:310172/148‑1 (MQ=255)
   cATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGa                                                                                                                                        >  2:412623/1‑149 (MQ=255)
              tCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAg                                                                                                                                   >  5:352141/1‑143 (MQ=255)
              tCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAg                                                                                                                                   <  6:352141/143‑1 (MQ=255)
                tGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGACCTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCcgc                                                                                                                           >  1:106184/1‑149 (MQ=255)
                tGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTCCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCcgc                                                                                                                           >  6:351514/1‑149 (MQ=255)
                  ccGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTCCGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCcgc                                                                                                                           <  4:245536/147‑1 (MQ=255)
                  ccGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCcgc                                                                                                                           <  3:117512/147‑1 (MQ=255)
                     aagaaAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGACGAAGCCGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAaca                                                                                                                      <  5:170313/149‑1 (MQ=255)
                     aagaaAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAaca                                                                                                                      >  2:19845/1‑149 (MQ=255)
                        aaaaTGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGAAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACAccc                                                                                                                   <  2:49788/149‑1 (MQ=255)
                            tGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTa                                                                                                               <  5:361249/149‑1 (MQ=255)
                              gCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTAct                                                                                                             >  4:346336/1‑149 (MQ=255)
                                aCTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTTGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTActgc                                                                                                           <  3:227703/149‑1 (MQ=255)
                                aCTGGAGCAGTTCAAACCGTGCATCTACCGCAAGCTGGAACTGCGTGGGCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTActgc                                                                                                           <  2:421552/149‑1 (MQ=255)
                                  tGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTActgct                                                                                                          >  8:178446/1‑148 (MQ=255)
                                          ttCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAg                                                                                                                                   <  5:341868/115‑1 (MQ=255)
                                          ttCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAg                                                                                                                                   >  6:341868/1‑115 (MQ=255)
                                          ttCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGt                                                                                                   <  8:41516/147‑1 (MQ=255)
                                            cAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTActgc                                                                                                           <  7:383710/137‑1 (MQ=255)
                                            cAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGAGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTActgc                                                                                                           >  8:383710/1‑137 (MQ=255)
                                             aaaCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAAcc                                                                                                     <  1:123540/142‑1 (MQ=255)
                                             aaaCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAAcc                                                                                                     >  2:123540/1‑142 (MQ=255)
                                             aaaCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCAcc                                                                                              >  2:328546/1‑149 (MQ=255)
                                             aaaCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCAcc                                                                                              <  3:346336/149‑1 (MQ=255)
                                             aaaCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGAAcc                                                                                              >  4:241731/1‑149 (MQ=255)
                                               aCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGTGCGAAGAAGCTGTCGTTTGGCATATc                                                                                                                                             <  2:70851/100‑1 (MQ=255)
                                               aCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGTGCGAAGAAGCTGTCGTTTGGCATATc                                                                                                                                             >  1:70851/1‑100 (MQ=255)
                                               aCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATc                                                                                                                                             <  8:45695/100‑1 (MQ=255)
                                               aCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATc                                                                                                                                             >  7:45695/1‑100 (MQ=255)
                                                ccGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGAc                                                                                           >  1:116515/1‑149 (MQ=255)
                                                   ttCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACAccc                                                                                                                   <  7:336637/122‑1 (MQ=255)
                                                   ttCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACAccc                                                                                                                   >  8:336637/1‑122 (MQ=255)
                                                           cGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCAc                                                                                    <  7:66161/145‑1 (MQ=255)
                                                            ggCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCt                                                                               >  5:144904/1‑149 (MQ=255)
                                                            ggCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCt                                                                               >  4:8483/1‑149 (MQ=255)
                                                             gCAAGCTGTAACTGCTTGGTATTGCTAACAACATTAAAGCTGCGAATAAAATGGTTGAGCGCTAAGAAGCTGTCGTTTGGAATATCCTGGAAGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTg                                                                              <  4:49031/149‑1 (MQ=255)
                                                             gCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTg                                                                              <  7:39236/149‑1 (MQ=255)
                                                                  cTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTAt                                                                         <  1:328546/149‑1 (MQ=255)
                                                                    ggAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCt                                                                               >  6:272758/1‑141 (MQ=255)
                                                                    ggAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCt                                                                               <  5:272758/141‑1 (MQ=255)
                                                                     gAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATcc                                                                       <  2:116515/148‑1 (MQ=255)
                                                                       aCTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAgg                                                                    <  1:412623/149‑1 (MQ=255)
                                                                          gCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCAt                                                                 >  2:306601/1‑149 (MQ=255)
                                                                           cGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGGCTGGGTATCCAGGCAtt                                                                <  2:106184/149‑1 (MQ=255)
                                                                              ggTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCAttt                                                               >  7:350366/1‑147 (MQ=255)
                                                                              ggTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCAttt                                                               <  8:350366/147‑1 (MQ=255)
                                                                               gTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAg                                                                     >  4:51049/1‑140 (MQ=255)
                                                                               gTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAg                                                                     <  3:51049/140‑1 (MQ=255)
                                                                               gTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGaa                                                            >  3:350071/1‑149 (MQ=255)
                                                                                  ttGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCg                                                         >  1:122855/1‑149 (MQ=255)
                                                                                    gCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCCCCGTCTGGGTATCCAGGCATTTGAACCgg                                                        >  2:313988/1‑148 (MQ=255)
                                                                                    gCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCgg                                                        >  2:196259/1‑148 (MQ=255)
                                                                                     cTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGt                                                       <  1:19845/148‑1 (MQ=255)
                                                                                        ccaccaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTAc                                                     <  4:335868/147‑1 (MQ=255)
                                                                                         caccaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGa                                                  <  2:122855/149‑1 (MQ=255)
                                                                                              ttAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTActg                                                                                                            >  4:387275/1‑86 (MQ=255)
                                                                                              ttAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTActg                                                                                                            <  3:387275/86‑1 (MQ=255)
                                                                                              ttAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTc                                                                                <  8:287737/114‑1 (MQ=255)
                                                                                              ttAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTc                                                                                >  7:287737/1‑114 (MQ=255)
                                                                                              ttAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCg                                               >  7:157287/1‑147 (MQ=255)
                                                                                              ttAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCg                                               >  4:344327/1‑147 (MQ=255)
                                                                                              ttAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGaa                                             <  3:8483/149‑1 (MQ=255)
                                                                                              ttAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGaa                                             <  4:308161/149‑1 (MQ=255)
                                                                                                aaaGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAg                                            <  1:306601/148‑1 (MQ=255)
                                                                                                        gaagaaAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTa                                   <  5:351514/149‑1 (MQ=255)
                                                                                                        gaagaaAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGCGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTa                                   >  2:362253/1‑149 (MQ=255)
                                                                                                              aaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCa                                                                      <  5:211109/108‑1 (MQ=255)
                                                                                                              aaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCa                                                                      >  6:211109/1‑108 (MQ=255)
                                                                                                                 ggTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTg                           >  7:32801/1‑148 (MQ=255)
                                                                                                                    tGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGActct                                                                                        <  3:108580/84‑1 (MQ=255)
                                                                                                                    tGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGActct                                                                                        >  4:108580/1‑84 (MQ=255)
                                                                                                                         gcgAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTg                  <  6:144904/149‑1 (MQ=255)
                                                                                                                           gaagaaGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGt                >  7:260588/1‑149 (MQ=255)
                                                                                                                                 gCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGc          >  4:26325/1‑149 (MQ=255)
                                                                                                                                 gCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGc          <  1:196259/149‑1 (MQ=255)
                                                                                                                                   tGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAAACGGTACTGATCGAAGGTAGAGCTATCCAGCTGCACCCGCTGGTTTGTGCgg        >  5:114930/1‑149 (MQ=255)
                                                                                                                                       gTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCata    <  8:260588/149‑1 (MQ=255)
                                                                                                                                        tttGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCata    >  8:167010/1‑148 (MQ=255)
                                                                                                                                          tGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAATCGGTACTGATCGAAGGTAGAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCAtata  <  6:114930/148‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                
CTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATA  >  NC_000913/4186436‑4186721

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