Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,094,906 T→C 5.3% I59V (ATA→GTA)  yggR ← putative PilT family AAA+ ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,094,9060TC 5.3% 230.8 / 2.9 75I59V (ATA→GTA) yggRputative PilT family AAA+ ATPase
Reads supporting (aligned to +/- strand):  ref base T (36/35);  new base C (2/2);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.81e-01

GCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACAC  >  NC_000913/3094760‑3095052
                                                                                                                                                  |                                                                                                                                                  
gcCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATAtt                                                                                                                                                   <  1:356720/148‑1 (MQ=255)
  cAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGcc                                                                                                                                                >  7:205151/1‑149 (MQ=255)
           cGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGGTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCGAGCGACACGGCAAAATCCAGCTGTCCATTCTCCAGCAATATTGCCCGCTGATc                                                                                                                                        <  3:137046/148‑1 (MQ=255)
             aGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTc                                                                                                                                     <  6:241578/149‑1 (MQ=255)
               cTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCAt                                                                                                                                   >  3:309001/1‑149 (MQ=255)
                  cGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATc                                                                                                                                  >  3:307516/1‑147 (MQ=255)
                   gggCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGGTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATAGCCCGCTGATCGTCATCCAg                                                                                                                               <  2:345559/149‑1 (MQ=255)
                            gAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCg                                                                                                                      >  7:275308/1‑149 (MQ=255)
                            gAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTGTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCg                                                                                                                      <  5:170205/149‑1 (MQ=255)
                               ggTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAg                                                                                                                   <  8:205151/149‑1 (MQ=255)
                               ggTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAg                                                                                                                   <  5:65207/149‑1 (MQ=255)
                                gTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAg                                                                                                                   >  5:15604/1‑148 (MQ=255)
                                gTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATACTGCCCGCTGATCGTCATCCAGc                                                                                                                              <  5:192788/137‑1 (MQ=255)
                                gTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATACTGCCCGCTGATCGTCATCCAGc                                                                                                                              >  6:197732/1‑137 (MQ=255)
                                gTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATACTGCCCGCTGATCGTCATCCAGc                                                                                                                              >  6:192788/1‑137 (MQ=255)
                                gTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATACTGCCCGCTGATCGTCATCCAGc                                                                                                                              <  5:197732/137‑1 (MQ=255)
                                 tAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTa                                                                                                                 <  8:203454/149‑1 (MQ=255)
                                   aCAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAg                                                                                                                <  4:105010/148‑1 (MQ=255)
                                         gTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTc                                                                                                          >  7:135195/1‑148 (MQ=255)
                                              gCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTc                                                                                                    >  6:299149/1‑149 (MQ=255)
                                              gCCAGCGAAATGCCATGCCGTTGCGCGAAAGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTc                                                                                                    >  8:119217/1‑149 (MQ=255)
                                                     aaaTGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCg                                                                                                         >  4:86749/1‑137 (MQ=255)
                                                            atgccGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCa                                                                                                                            >  6:221555/1‑111 (MQ=255)
                                                               ccGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGgc                                                                                    <  4:330000/148‑1 (MQ=255)
                                                                    gcgcGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGc                                                                               <  5:299149/148‑1 (MQ=255)
                                                                      gcgAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATc                                                                                                                                        <  2:78134/89‑1 (MQ=255)
                                                                      gcgAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATc                                                                                                                                        >  1:78134/1‑89 (MQ=255)
                                                                      gcgAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTc                                                                                                                                     <  6:155604/92‑1 (MQ=255)
                                                                      gcgAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTc                                                                                                                                     >  5:155604/1‑92 (MQ=255)
                                                                         aaTGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCAt                                                                         <  6:305808/149‑1 (MQ=255)
                                                                             cgcTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTg                                                                     >  6:173687/1‑149 (MQ=255)
                                                                                   cACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCcgcg                                                                >  7:339230/1‑148 (MQ=255)
                                                                                    aCGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGaa                                                              >  3:325065/1‑149 (MQ=255)
                                                                                     cGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAAt                                                             <  2:354999/149‑1 (MQ=255)
                                                                                     cGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAAt                                                             <  8:228282/149‑1 (MQ=255)
                                                                                            gCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGCCGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCAt                                                      >  5:323114/1‑149 (MQ=255)
                                                                                            gCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCAt                                                      >  4:94481/1‑149 (MQ=255)
                                                                                            gCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCAt                                                      >  8:73146/1‑149 (MQ=255)
                                                                                             cTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACTGCGCAGCTTCCATTCTGCCGCGAATGCGCAATc                                                     >  6:62303/1‑149 (MQ=255)
                                                                                                      cGCCAGCGACACGTCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCGGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGc                                             <  6:302295/148‑1 (MQ=255)
                                                                                                      cGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTg                                                  <  6:317796/143‑1 (MQ=255)
                                                                                                      cGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTg                                                  >  5:317796/1‑143 (MQ=255)
                                                                                                        cgaGCGACACGGCAAAATCCTGCTGACCTTTGTCCAGCTATATTGCCGGCTGTTCGTCTTCCAGCCACTCCCGCAGTTGCTCTTCGACGTCCTGCGGATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAg                                          <  6:323114/147‑1 (MQ=255)
                                                                                                          aGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGgc                                        >  3:188192/1‑149 (MQ=255)
                                                                                                                  ggCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGca                                  >  5:333054/1‑147 (MQ=255)
                                                                                                                   gCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAg                               <  8:339230/149‑1 (MQ=255)
                                                                                                                   gCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAg                               <  5:173687/149‑1 (MQ=255)
                                                                                                                     aaaaTCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCa                                           <  5:210925/135‑1 (MQ=255)
                                                                                                                     aaaaTCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCa                                           >  6:210925/1‑135 (MQ=255)
                                                                                                                     aaaaTCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGgt                             <  4:309001/149‑1 (MQ=255)
                                                                                                                     aaaaTCCAGCTGACCATTCTCCAGCAATATTGCCCGCTCATCGTCATCCAGCCACTCCCCCAGAAGCTCTACGACGTCCCGCGTATCAAACGGCGCACCTTCCATTCTGCCGCGAATGCGCCATCGTACGGGCAAGGCGCTGAACAGgt                             >  6:182817/1‑149 (MQ=255)
                                                                                                                      aaaTCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGgtg                            <  3:94481/149‑1 (MQ=255)
                                                                                                                       aaTCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGgtgt                           <  4:307516/149‑1 (MQ=255)
                                                                                                                       aaTCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGgtg                            >  2:79478/1‑148 (MQ=255)
                                                                                                                        aTCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGgtg                            >  3:301045/1‑147 (MQ=255)
                                                                                                                             gCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATcc                     <  4:188192/149‑1 (MQ=255)
                                                                                                                                 aCCTTTCTCCAGCAATATTGCCCGCTGATAGTCATCCAGGCACTCCCGAAGTAGCTCTTCGACGTCAGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCga                   <  6:352138/147‑1 (MQ=255)
                                                                                                                                 aCCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGgt                             <  3:121626/137‑1 (MQ=255)
                                                                                                                                 aCCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGgt                             >  4:121626/1‑137 (MQ=255)
                                                                                                                                 aCCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCgaga                 >  5:297579/1‑149 (MQ=255)
                                                                                                                                 aCCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCga                   <  4:325065/147‑1 (MQ=255)
                                                                                                                                 aCCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCga                   <  5:182817/147‑1 (MQ=255)
                                                                                                                                   cATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAg                               >  3:21173/1‑133 (MQ=255)
                                                                                                                                   cATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAg                               <  4:21173/133‑1 (MQ=255)
                                                                                                                                   cATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACg               >  3:270963/1‑149 (MQ=255)
                                                                                                                                     ttCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGt              >  3:207951/1‑148 (MQ=255)
                                                                                                                                      tctcCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTGTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGcaca                                <  4:176503/129‑1 (MQ=255)
                                                                                                                                      tctcCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGcaca                                >  3:176503/1‑129 (MQ=255)
                                                                                                                                        tcCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAg                         >  4:257071/1‑134 (MQ=255)
                                                                                                                                        tcCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAg                         <  3:257071/134‑1 (MQ=255)
                                                                                                                                        tcCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATg          >  3:14598/1‑149 (MQ=255)
                                                                                                                                           aGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCgg                                               <  6:189029/109‑1 (MQ=255)
                                                                                                                                           aGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCgg                                               >  5:189029/1‑109 (MQ=255)
                                                                                                                                             cAATATTGCCCGCTGAGCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTa     <  3:258185/149‑1 (MQ=255)
                                                                                                                                              aaTATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGgt                             <  7:216395/124‑1 (MQ=255)
                                                                                                                                              aaTATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGgt                             >  8:216395/1‑124 (MQ=255)
                                                                                                                                                tatTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGTCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTacac  <  6:15604/149‑1 (MQ=255)
                                                                                                                                                tatTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCcgc                                                                 >  7:73829/1‑86 (MQ=255)
                                                                                                                                                tatTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCcgc                                                                 <  8:73829/86‑1 (MQ=255)
                                                                                                                                                tatTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCATTacac  >  1:326596/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCCAAGCTGTTCGAGCTGCGGGCAGTGCGAAGGTAACAGCCGTAACGCCAGCGAAATGCCATGCCGTTGCGCGAATGCGCTGCCACGCAATCGCTGGTTTTCCGCCAGCGACACGGCAAAATCCAGCTGACCATTCTCCAGCAATATTGCCCGCTGATCGTCATCCAGCCACTCCCGCAGTAGCTCTTCGACGTCCGGCGTATCAAACGGCGCAGCTTCCATTCTGCCGCGAATGCGCCATCGTGCGGGCCAGGCGCTGCACAGGTGTAGATCCGAGACGTTATGCTTTACAC  >  NC_000913/3094760‑3095052

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