Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,808,813 C→T 6.2% R92C (CGC→TGC)  waaA → 3‑deoxy‑D‑manno‑octulosonic‑acid transferase (KDO transferase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,808,8130CT 6.2% 320.1 / 7.3 97R92C (CGC→TGC) waaA3‑deoxy‑D‑manno‑octulosonic‑acid transferase (KDO transferase)
Reads supporting (aligned to +/- strand):  ref base C (45/46);  new base T (3/3);  total (48/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

CCATCCGCTAAAACCAGGCGGCATTATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGA  >  NC_000913/3808668‑3808951
                                                                                                                                                 |                                                                                                                                          
ccATCCGCTAAAACCAGGCGGCATTATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAgcgc                                                                                                                                          <  8:241766/148‑1 (MQ=255)
   tCCGCTAAAACCAGGCGGCATTATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTa                                                                                                                                       <  8:109147/148‑1 (MQ=255)
    ccGCTAAAACCAGGCGGCATTATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTa                                                                                                                                       >  1:115441/1‑147 (MQ=255)
     cGCTAAAACCAGGCGGCATTATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTAc                                                                                                                                      <  5:285356/147‑1 (MQ=255)
      gCTAAAACCAGGCGGCATTATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACaa                                                                                                                                    >  1:270270/1‑148 (MQ=255)
       cTAAAACCAGGCGGCATTATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAAt                                                                                                                                   <  2:238305/148‑1 (MQ=255)
         aaaaCCAGGCGGCATTATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCg                                                                                                                                 >  3:331022/1‑148 (MQ=255)
          aaaCCAGGCGGCATTATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCAGCGTTATCCTGATTTACCGATTACCGTAACACCCATGCCGCCAACCGGTTCGGAGCGCGTACAATCgg                                                                                                                                >  8:385797/1‑148 (MQ=255)
              cAGGCGGCATTATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCttt                                                                                                                            <  4:86149/148‑1 (MQ=255)
              cAGGCGGCATTATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCttt                                                                                                                            <  2:115441/148‑1 (MQ=255)
                  cggcATTATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCATCCGGTTCGGAGCGCGTACATTCGGCTTTCgg                                                                                                                         >  5:128442/1‑147 (MQ=255)
                  cggcATTATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGt                                                                                                                                        >  2:233384/1‑132 (MQ=255)
                  cggcATTATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGt                                                                                                                                        <  1:233384/132‑1 (MQ=255)
                  cggcATTATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCgg                                                                                                                         >  5:198215/1‑147 (MQ=255)
                     cATTATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAg                                                                                                                     <  6:198215/148‑1 (MQ=255)
                      aTTATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAgg                                                                                                                    <  2:270270/148‑1 (MQ=255)
                        tATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAg                                                                                                                     >  7:338909/1‑145 (MQ=255)
                        tATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAg                                                                                                                     <  8:338909/145‑1 (MQ=255)
                        tATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGAt                                                                                                                  >  4:117743/1‑148 (MQ=255)
                        tATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGAt                                                                                                                  >  3:69743/1‑148 (MQ=255)
                        tATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGAt                                                                                                                  >  1:20509/1‑148 (MQ=255)
                            ctgcACTCCGTCTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGtt                                                                                                               <  4:331022/147‑1 (MQ=255)
                            ctgcACTCCGTCTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGtt                                                                                                               >  5:65978/1‑147 (MQ=255)
                            ctgcACTCCGTCTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGtt                                                                                                               <  6:65978/147‑1 (MQ=255)
                            ctgcACTCCGTCTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTGCGGGAAGGa                                                                                                                   <  4:254088/143‑1 (MQ=255)
                            ctgcACTCCGTCTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTACGGGAAGGa                                                                                                                   >  3:254088/1‑143 (MQ=255)
                                 ctccggctccggCGGGGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCggg                                                                                                                        <  6:205806/133‑1 (MQ=255)
                                  tccgtctccgtcGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTAc                                                                                                                                      <  7:401609/118‑1 (MQ=255)
                                  tccgtctccgtcGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTAc                                                                                                                                      >  8:401609/1‑118 (MQ=255)
                                  tccgtctccgtcGGTGAAACTCTGGCGGCAATCCCGTAGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAAACGGCTTTCGGGa                                                                                                                       >  5:205806/1‑133 (MQ=255)
                                  tccgtctccgtcGGTGAAACTCTGGCGCCAATCCAGTTGGTTCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAAACATGACGCCAACCGGATCGCACCGCCTACAATCGGCTTTCGGGAAGGATGTTCAGCACg                                                                                                        >  3:341800/1‑148 (MQ=255)
                                   ccgtctccgtcGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGt                                                                                                       >  1:62607/1‑148 (MQ=255)
                                    aggctccgtcGGTGTAACTCTGGCGGAAATGCCGTTGGTGCGCGCGCGGCTTTATCGTTATCCTGATTTACCGTTTACCGTAAGAACCATGACGCAATCCTGTTCGGAGCGCGTACAATCGTCTTTCTGGAAGGATGTTCTGGACGt                                                                                                       <  2:150778/144‑1 (MQ=255)
                                          cgtcGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTAt                                                                                                   <  2:20509/145‑1 (MQ=255)
                                          cgtcGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTAACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTAt                                                                                                   <  3:117743/145‑1 (MQ=255)
                                            tcGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGTGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGcc                                                                                              >  1:187619/1‑148 (MQ=255)
                                            tcGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGGTATCCTGATTTACCGATTACCGTAACAACCATGACGCCCACCGGTTCGGAGTGCGAACAATCGGCTTACGGGAAGCATGTTCAGCACGTTTATCTGcc                                                                                              >  1:193236/1‑148 (MQ=255)
                                             cGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCg                                                                                             >  8:20050/1‑148 (MQ=255)
                                             cGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCg                                                                                             >  2:77074/1‑148 (MQ=255)
                                             cGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACATTTATCTGCCg                                                                                             >  4:165249/1‑148 (MQ=255)
                                             cGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCG‑AACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTa                                                                                                    >  7:45326/1‑140 (MQ=255)
                                             cGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCG‑AACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTa                                                                                                    <  8:45326/140‑1 (MQ=255)
                                                  aaaCTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGa                                                                                        <  7:20050/148‑1 (MQ=255)
                                                  aaaCTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACATTTATCTGCCGTATGa                                                                                        <  3:165249/148‑1 (MQ=255)
                                                   aaCTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATc                                                                                                                                  >  4:118663/1‑105 (MQ=255)
                                                   aaCTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATc                                                                                                                                  <  3:118663/105‑1 (MQ=255)
                                                   aaCTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGAt                                                                                       >  5:297394/1‑148 (MQ=255)
                                                    aCTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTa                                                                                                    <  8:280159/134‑1 (MQ=255)
                                                    aCTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTa                                                                                                    >  7:280159/1‑134 (MQ=255)
                                                       ctGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTa                                                                                                    >  6:36936/1‑131 (MQ=255)
                                                       ctGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTa                                                                                                    <  5:36936/131‑1 (MQ=255)
                                                         ggcggcAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCa                                                                                                             <  2:76424/120‑1 (MQ=255)
                                                         ggcggcAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCa                                                                                                             >  1:76424/1‑120 (MQ=255)
                                                          gcggcAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGTGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCg                                                                                <  2:187619/148‑1 (MQ=255)
                                                          gcggcAATCCCGGTGGTGCGCGCGCGGCGGCATCGGTATCCTGATTTACCGATTACCGTAACATCCTTGACGCCTACCGTTTCGGAGTGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATGTGCCGTATGATCTGCCCg                                                                                <  2:193236/148‑1 (MQ=255)
                                                            ggcAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGAt                                                                              <  7:165006/148‑1 (MQ=255)
                                                             gcAATCCCGTTGGTGCGCGCGCTGGGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGtt                                                                                                      <  1:229443/123‑1 (MQ=255)
                                                             gcAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGtt                                                                                                      >  2:229443/1‑123 (MQ=255)
                                                                 tCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGTGCGTACAATCGGCtt                                                                                                                             >  7:268651/1‑96 (MQ=255)
                                                                 tCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGTGCGTACAATCGGCtt                                                                                                                             <  8:268651/96‑1 (MQ=255)
                                                                 tCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAgg                                                                                                                    <  5:381761/105‑1 (MQ=255)
                                                                 tCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAgg                                                                                                                    >  6:381761/1‑105 (MQ=255)
                                                                 tCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGa                                                                                        >  7:202495/1‑133 (MQ=255)
                                                                 tCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGa                                                                                        <  8:202495/133‑1 (MQ=255)
                                                                  cccGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACt                                                                         >  2:125186/1‑147 (MQ=255)
                                                                  cccGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACt                                                                         >  5:144049/1‑147 (MQ=255)
                                                                      ttGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAAcc                                                                    >  6:254541/1‑148 (MQ=255)
                                                                      ttGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAAcc                                                                    >  3:214829/1‑148 (MQ=255)
                                                                          tGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACCGttt                                                                <  6:144049/148‑1 (MQ=255)
                                                                                     cGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACCGTTTCCTGAATaa                                                       <  2:62607/146‑1 (MQ=255)
                                                                                      gTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACCGTTTCCTGAATAAAGt                                                    >  2:286755/1‑148 (MQ=255)
                                                                                      gTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTGCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACCGTTTCCTGAATAAAGt                                                    <  4:69743/148‑1 (MQ=255)
                                                                                        cATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGc                                                                                                           >  4:14576/1‑91 (MQ=255)
                                                                                        cATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGc                                                                                                           <  3:14576/91‑1 (MQ=255)
                                                                                        cATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGACCGATGCACTCAACCGTTTCCTGAATAAAGTc                                                   >  8:166362/1‑147 (MQ=255)
                                                                                          tCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACCGTTTCCTGAATAAAGTCGAc                                                <  5:254541/148‑1 (MQ=255)
                                                                                                tCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACCGTTTCCTGTATAAAGTCGACCCTaa                                           <  6:297394/147‑1 (MQ=255)
                                                                                                tCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACCGTTTCCTGAATAAAGTCGACCCTaa                                           <  1:286755/147‑1 (MQ=255)
                                                                                                   tGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGAt                                                                                       <  5:162150/100‑1 (MQ=255)
                                                                                                   tGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGAt                                                                                       >  6:162150/1‑100 (MQ=255)
                                                                                                     aTTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTg                                       <  1:125186/146‑1 (MQ=255)
                                                                                                      tttACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGgtg                                    <  1:77074/148‑1 (MQ=255)
                                                                                                            gATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGAtt                              <  7:385797/148‑1 (MQ=255)
                                                                                                              ttACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTg                                 <  7:172279/143‑1 (MQ=255)
                                                                                                               tACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATTAACTCAACCGTTTAATGAATAAAGTCGACCCTAAACTGGTGTTGa                                >  8:172279/1‑143 (MQ=255)
                                                                                                                   gTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCTGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACCGTTTCCTGa                                                           >  6:31513/1‑112 (MQ=255)
                                                                                                                   gTAACAACCATGACGCCAACCGGTGCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACCGTTTCCTGa                                                           <  5:31513/112‑1 (MQ=255)
                                                                                                                    tAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAAc                      >  6:158548/1‑148 (MQ=255)
                                                                                                                     aacaacCATGACGCCAACCGGTTCGTAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAAcc                     >  8:387273/1‑148 (MQ=255)
                                                                                                                     aacaacCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAAcc                     >  7:180813/1‑148 (MQ=255)
                                                                                                                       caacCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGa                   <  7:387273/148‑1 (MQ=255)
                                                                                                                              gACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACCGTTTCCTg                                                            <  6:358223/100‑1 (MQ=255)
                                                                                                                              gACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACCGTTTCCTg                                                            >  5:358223/1‑100 (MQ=255)
                                                                                                                              gACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATgg            >  8:394814/1‑148 (MQ=255)
                                                                                                                                 gCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATgg                          <  4:185314/131‑1 (MQ=255)
                                                                                                                                 gCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATgg                          >  3:185314/1‑131 (MQ=255)
                                                                                                                                        ggTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGa  <  5:158548/148‑1 (MQ=255)
                                                                                                                                        ggTTCGGAGCGCGTACAATCGGCTTGCGGGAAGGACGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATTCACTGAACCGTTTCATTAATTAAGTCGTCCGTAATCTGGTGTTGATTATGGATACCGAACTCTGGCCTAACCTGa  <  4:341800/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                          
CCATCCGCTAAAACCAGGCGGCATTATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGA  >  NC_000913/3808668‑3808951

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