Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,709,445 T→C 5.1% R190G (AGG→GGG)  rpoE ← RNA polymerase sigma E factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,709,4450TC 5.1% 252.9 / 2.8 80R190G (AGG→GGG) rpoERNA polymerase sigma E factor
Reads supporting (aligned to +/- strand):  ref base T (38/37);  new base C (2/2);  total (41/39)
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.36e-01

GCTTTCCCAGGTTTTCTGCATTTCTGGGTTATGAGCCAGTTCGTTAAGCAGCTCACTATCCAGCGTTTCGCCATCCATTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTATCTCTTCATAGCTCAGGCCATCCAGCTCCCGCAAGGTTATTGCCATGCGTAAATCT  >  NC_000913/2709300‑2709593
                                                                                                                                                 |                                                                                                                                                    
gCTTTCCCAGGTTTTCTGCATTTCTGGGTTATGAGCCAGTTCGTTAAGCAGCTCACTATCCAGCGTTTCGCCATCCATTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGAt                                                                                                                                                   <  5:211933/149‑1 (MQ=255)
 cTTTCCCAGGTTTTCTGCATTTCTGGGTTATGAGCCAGTTCGTTAAGCAGCTCACTATCCAGCGTTTCGCCATCCATTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATCCCCTTATCCAGTATCCCGCTATCGTCAACGCCTGAt                                                                                                                                                   >  2:394983/1‑148 (MQ=255)
      ccAGGTTTTCTGCATTTCTGGGTTATGAGCCAGTTCGTTAAGCAGCTCACTATCCAGCGTTTCGCCATCCATTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCgg                                                                                                                                             >  7:374789/1‑149 (MQ=255)
          gTTTTCTGCATTTCTGGGTTATGAGCCAGTTCGTTAAGCAGCTCACTATCCAGCGTTTCGCCATCCATTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTg                                                                                                                                          >  8:102033/1‑148 (MQ=255)
             ttCTGCATTTCTGGGTTATGAGCCAGTTCGTTAAGCAGCTCACTATCCAGCGTTTCGCCATCCATTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACt                                                                                                                                      >  8:51162/1‑149 (MQ=255)
                  cATTTCTGGGTTATGAGCCAGTTCGTTAAGCAGCTCACTATCCAGCGTTTCGCCATCCATTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACttt                                                                                                                                    <  1:217452/146‑1 (MQ=255)
                  cATTTCTGGGTTATGAGCCAGTTCGTTAAGCAGCTCACTATCCAGCGTTTCGCCATCCATTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACttt                                                                                                                                    <  8:257704/146‑1 (MQ=255)
                  cATTTCTGGGTTATGAGCCAGTTCGTTAAGCAGCTCACTATCCAGCGTTTCGCCATCCATTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACttt                                                                                                                                    <  4:82908/146‑1 (MQ=255)
                  cATTTCTGGGTTATGAGCCAGTTCGTTAAGCAGCTCACTATCCAGCGTTTCGCCATCCATTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACttt                                                                                                                                    <  3:389849/146‑1 (MQ=255)
                  cATTTCTGGGTTATGAGCCAGTTCGTTAAGCAGCTCACTATCCAGCGTTTCGCCATCCATTAAAGCGGAAAGTTGTGCTCTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACttt                                                                                                                                    <  5:323923/146‑1 (MQ=255)
                        tGGGTTATGAGCCAGTTCGTTAAGCAGCTCACTATCCAGCGTTTCGCCATCCATTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCaa                                                                                                                            >  6:138654/1‑148 (MQ=255)
                        tGGGTTATGAGCCAGTTCGTTAAGCAGCTCACTATCCAGCGTTTCGCCATCCATTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAAt                                                                                                                           <  4:226963/149‑1 (MQ=255)
                         gggTTATGAGCCAGTTCGTTAAGCATCTCACTATCCAGCGTTTCGCCATCCATTAAAGCGGAAACTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCaa                                                                                                                            >  8:112891/1‑147 (MQ=255)
                          ggTTATGAGCCAGTTCGTTAAGCAGCTCACTATCCAGCGTTTCGCCATCCATTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAg                                                                                                                         >  8:145962/1‑149 (MQ=255)
                            ttATGAGCCAGTTCGTTAAGCAGCTCACTATCCAGCGTTTCGCCATCCATTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCt                                                                                                                       <  7:51162/149‑1 (MQ=255)
                               tGAGCCAGTTCGTTAAGCAGCTCACTATCCAGCGTTTCGCCATCCATTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATa                                                                                                                          <  6:124641/143‑1 (MQ=255)
                               tGAGCCAGTTCGTTAAGCAGCTCACTATCCAGCGTTTCGCCATCCATTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATa                                                                                                                          >  5:124641/1‑143 (MQ=255)
                               tGAGCCAGTTCGTTAAGCAGCTCACTATCCAGCGTTTCGCCATCCATTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTcc                                                                                                                    >  7:360529/1‑149 (MQ=255)
                                  gCCAGTTCGTTAAGCAGCTCACTATCCAGCGTTTCGCCATCCATTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCctc                                                                                                                 >  1:161632/1‑149 (MQ=255)
                                  gCCAGTTCGTTAAGCAGCTCACTATCCAGCGTTTCGCCATCCATTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTGTGTTAGCAATAGCTTCCctc                                                                                                                 >  4:216684/1‑149 (MQ=255)
                                     aGTTCGTTAAGCAGCTCACTATCCAGCGTTTCGCCATCCATTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCctcgct                                                                                                              >  4:155158/1‑149 (MQ=255)
                                      gTTCGTTAAGCAGCTCACTATCCAGCGTTTCGCCATCCATTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCctcgctc                                                                                                             <  7:102033/149‑1 (MQ=255)
                                      gTTCGTTAAGCAGCTCACTATCCAGCGTTTCGCCATCCATTAAAGCGGAAAGTGGTTCTTTCTGCATGCCTAATACGCTTATCCAGTATCCCGCTATCGTCAACGGCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCctcgctc                                                                                                             <  1:97495/149‑1 (MQ=255)
                                             aaGCAGCTCACTATCCAGCGTTTCGCCATCCATTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGa                                                                                                       <  8:374789/148‑1 (MQ=255)
                                              agcagcTCACTATCCAGCGTTTCGCCATCCATTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGAt                                                                                                      <  3:155158/148‑1 (MQ=255)
                                                cagcTCACTATCCAGCGTTTCGCCATCCATTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCctc                                                                                                                 <  7:319564/135‑1 (MQ=255)
                                                cagcTCACTATCCAGCGTTTCGCCATCCATTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCctc                                                                                                                 >  8:319564/1‑135 (MQ=255)
                                                 agcTCACTATCCAGCGTTTCGCCATCCATTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGt                                                                                                  >  4:214514/1‑149 (MQ=255)
                                                  gcTCACTATCCAGCGTTTCGCCATCCATTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGt                                                                                                  >  5:217366/1‑148 (MQ=255)
                                                  gcTCACTATCCAGCGTTTCGCCATCCATTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGt                                                                                                  >  3:286623/1‑148 (MQ=255)
                                                     cACTATCCAGCGTTTCGCCATCCATTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAAc                                                                                              >  1:43275/1‑149 (MQ=255)
                                                          tcaaGCGTTTCGCCATCCATTAAAGCCGAAAGTTGTTCTTTCTGCATGCCAAATACCCTTATCCAGTATCCCGCTTTCGTCAACGCCTCATAAGCCGTTGAACTTTCTTATCAATAGCTTCCCTAGCTCGGAAGATACCTGAACGCAcc                                                                                         >  5:283737/4‑149 (MQ=255)
                                                             aGCGTTTCGCCATCCATTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGta                                                                                      <  8:360529/149‑1 (MQ=255)
                                                              gCGTTTCGCCATCCATTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGGCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGtac                                                                                     <  2:114466/149‑1 (MQ=255)
                                                                 tttCGCCATCCATTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCAGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGtaccta                                                                                  >  1:20941/1‑149 (MQ=255)
                                                                      ccatccatTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGa                                                                             <  7:145962/149‑1 (MQ=255)
                                                                        atccatTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCTACGCCTGATAAGCGGTTGAACTTTGTTATCATTAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGAc                                                                            >  4:10493/1‑148 (MQ=255)
                                                                         tccatTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACaa                                                                          >  7:293385/1‑149 (MQ=255)
                                                                            atTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATcc                                                                       <  5:190576/149‑1 (MQ=255)
                                                                            atTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATc                                                                        >  3:138990/1‑148 (MQ=255)
                                                                               aaaGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCAt                                                                     <  8:404559/148‑1 (MQ=255)
                                                                                   cGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAg                                                                >  6:382301/1‑149 (MQ=255)
                                                                                          ttgttCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTAt                                                         <  5:382301/149‑1 (MQ=255)
                                                                                             ttctttctGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTAtctc                                                      >  8:411378/1‑149 (MQ=255)
                                                                                             ttctttctGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAAAAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTAtct                                                       <  8:338534/148‑1 (MQ=255)
                                                                                               ctttctGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTATctct                                                     <  8:259621/148‑1 (MQ=255)
                                                                                                      cATGCGTAATACCCTTTTCCAGTATCCCGCTATATTTAACGGCTGATAAGCGGTTGAACTTTGTTATTTATAGCTTCCCTCGCTCGGAATATATGTGAACGGAGCGTACCTACCGGACAATCGATGATAGCGGTTATTTCTTCATa                                                <  6:283737/146‑1 (MQ=255)
                                                                                                      cATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTATCTCTTCATa                                                <  2:43275/146‑1 (MQ=255)
                                                                                                        tGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTATCTCTTCATAGCt                                             >  8:384526/1‑147 (MQ=255)
                                                                                                        tGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTATCTCTTCATAGCt                                             >  7:81846/1‑147 (MQ=255)
                                                                                                          ccTAATACGCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTAGCTACCGGACAATCCATGATAGCGGCTATCTCTTCATAGCTCAg                                          <  7:112891/148‑1 (MQ=255)
                                                                                                           cTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTATCTCTTCATAGCTCAGGc                                        >  3:123050/1‑149 (MQ=255)
                                                                                                            tAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTATCTCTTCATAGCTCAGGc                                        >  1:5252/1‑148 (MQ=255)
                                                                                                            tAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCCGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTATCTCTTCATAGCTCAg                                          <  2:209126/146‑1 (MQ=255)
                                                                                                            tAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCCGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTATCTCTTCATAGCTCAg                                          >  3:67693/1‑146 (MQ=255)
                                                                                                            tAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCCGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTATCTCTTCATAGCTCAg                                          >  1:209126/1‑146 (MQ=255)
                                                                                                            tAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCCGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTATCTCTTCATAGCTCAg                                          <  4:67693/146‑1 (MQ=255)
                                                                                                               tACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTATCTCTTCATAGCTCAGGCCATc                                    <  2:20941/149‑1 (MQ=255)
                                                                                                                 cccTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTATCTCTTCATAGCTCAGGCCATCCa                                  <  4:286623/149‑1 (MQ=255)
                                                                                                                 cccTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTATCTCTTCATAGCTCAGGCCATCCa                                  <  6:217366/149‑1 (MQ=255)
                                                                                                                  ccTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGa                                                                   <  8:58639/115‑1 (MQ=255)
                                                                                                                  ccTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGa                                                                   >  7:58639/1‑115 (MQ=255)
                                                                                                                  ccTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTATCTCTTCATAGCTCAGGCCATcccg                                 >  1:103545/1‑147 (MQ=255)
                                                                                                                  ccTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTATCTCTTCATAGCTCAGGCCATCCAg                                 >  2:415281/1‑149 (MQ=255)
                                                                                                                  ccTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTATCTCTTCATAGCTCAGGCCATCCAg                                 >  6:202524/1‑149 (MQ=255)
                                                                                                                        tcAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCTGACAATCCATGATATCGGCTATCTCTTCATAGCTCAGGCCATCCAGCTccc                            <  1:219182/147‑1 (MQ=255)
                                                                                                                         cAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAACAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTATCTCTTCATAGCTCAGGCCATCCAGCTCCCg                           >  2:198767/1‑148 (MQ=255)
                                                                                                                          aGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTATCTCTTCATAGCTCAGGCCATCCAGCTccc                            >  2:219182/1‑146 (MQ=255)
                                                                                                                           gTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTATCTCTTCATAGCTCAGGCCATCCAGCTCCCGCaa                        <  7:239917/149‑1 (MQ=255)
                                                                                                                            tATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTATCTCTTCATAGCTCAGGCCATCCAGCTCCCGCAAg                       <  1:88991/149‑1 (MQ=255)
                                                                                                                              tCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTATCTCTTCATAGCTCAGGCCATCCAGCTCCCGCAAgg                      <  8:293385/148‑1 (MQ=255)
                                                                                                                                ccGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTATCTCTTCATAGCTCAGGCCATCCAGCTCCCGCAAGGTTa                   <  6:310595/149‑1 (MQ=255)
                                                                                                                                 cGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTATCTCTTCATAGCTCAGGCCATCCAGCTCCCGCAAGGTTAt                  <  1:394983/149‑1 (MQ=255)
                                                                                                                                   cTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTATCTCTTCATAGCTCAGGCCATCCAGCTCCCGCAAGGTTATTg                >  2:349432/1‑149 (MQ=255)
                                                                                                                                     aTCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTATCTCTTCATAGCTCAGGCCATCCAGCTCCCGCAAGGTTATTGc               >  8:350109/1‑148 (MQ=255)
                                                                                                                                      tCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTATCTCTTCATAGCTCAGGCCATCCAGCTCCCGCAAGGTTATTGcc              <  1:178608/148‑1 (MQ=255)
                                                                                                                                         tCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTATCTCTTCATAGCTCAGGCCATCCAGCTCCCGCAAGGTTATTGCCATGc          >  1:327757/1‑149 (MQ=255)
                                                                                                                                           aaCGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTATCTCTTCATAGCTCAGGCCATCCAGCTCCCGCAAGGTTATTGCCATGCg         >  2:170699/1‑148 (MQ=255)
                                                                                                                                             cGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTATCTCTTCATAGCTCAGGCCATCCAGCTCCCGCAAGGTTATTGCCATGCGTaa      <  5:138654/149‑1 (MQ=255)
                                                                                                                                              gCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTATCTCTTCATAGCTCAGGCCATCCAGCTCCCGCAAGGTTATTGCCATGCGTaaa     >  5:258802/1‑149 (MQ=255)
                                                                                                                                                 tGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACACACCGTACCTACCGGACAATCCATGATAGCGGCTATCTCTTCATAGCTCAGGCCATCCATCTCCCGCAAGGTTATTGCCATGCGTAAATCt  >  6:29026/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GCTTTCCCAGGTTTTCTGCATTTCTGGGTTATGAGCCAGTTCGTTAAGCAGCTCACTATCCAGCGTTTCGCCATCCATTAAAGCGGAAAGTTGTTCTTTCTGCATGCCTAATACCCTTATCCAGTATCCCGCTATCGTCAACGCCTGATAAGCGGTTGAACTTTGTTATCAATAGCTTCCCTCGCTCGGAAGATACGTGAACGCACCGTACCTACCGGACAATCCATGATAGCGGCTATCTCTTCATAGCTCAGGCCATCCAGCTCCCGCAAGGTTATTGCCATGCGTAAATCT  >  NC_000913/2709300‑2709593

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