Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,769,727 Δ1 bp 5.5% coding (654/1113 nt) ydiK → UPF0118 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,769,7270G. 5.5% 250.5 / 7.0 73coding (654/1113 nt)ydiKUPF0118 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (38/31);  new base . (2/2);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.71e-01

TGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGA  >  NC_000913/1769579‑1769869
                                                                                                                                                    |                                                                                                                                              
tGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGtgtg                                                                                                                                                >  7:42816/1‑149 (MQ=255)
        ctctTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCTTTCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGAc                                                                                                                                        >  8:294595/1‑149 (MQ=255)
           ttCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGc                                                                                                                                     >  2:165442/1‑149 (MQ=255)
                tGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAg                                                                                                                                <  3:216649/149‑1 (MQ=255)
                  cccTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCGATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGGGGTGGTGACGGCGgtagta                                                                                                                              >  4:234617/1‑149 (MQ=255)
                  cccTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTa                                                                                                                              <  8:42816/149‑1 (MQ=255)
                  cccTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTa                                                                                                                              <  3:173415/149‑1 (MQ=255)
                        tgTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAgg                                                                                                                          >  1:125105/1‑147 (MQ=255)
                                 gcgGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTgtgg                                                                                                                                               >  4:130401/1‑117 (MQ=255)
                                 gcgGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTgtgg                                                                                                                                               <  3:130401/117‑1 (MQ=255)
                                     tGAACAGGTGGCACAAGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTa                                                                                                           >  2:190057/1‑149 (MQ=255)
                                        aCAGGTGGCACAAGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCg                                                                                                        <  1:317375/149‑1 (MQ=255)
                                          aGGTGGCACAAGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGc                                                                                                      <  2:390133/149‑1 (MQ=255)
                                           ggTGGCACAAGGCATTCGCCATTTTGCAACCCGTCTGGGAGGCGTTCGCTGAGTTGCCGGCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCTCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGc                                                                                                      <  4:78172/148‑1 (MQ=255)
                                           ggTGGCACAAGGCATTCGCAATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTTTCCGCGCGGTGGCGCTGGTGGGGGTGGTGACGGCGTTAGTACAGGCAGTGGTTGGCGGTATCGccc                                                                                                     >  3:78172/1‑149 (MQ=255)
                                            gTGGCACAAGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCt                                                                                                    >  4:337801/1‑149 (MQ=255)
                                            gTGGCACAAGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATAGGCCt                                                                                                    >  7:167875/1‑149 (MQ=255)
                                                cacaAGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGcc                                                                                                <  3:337801/149‑1 (MQ=255)
                                                 acaAGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCg                                                                                               >  2:41427/1‑149 (MQ=255)
                                                  caAGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCg                                                                                               >  6:217405/1‑148 (MQ=255)
                                                        aTTCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCgg                                                                                        >  2:293441/1‑149 (MQ=255)
                                                            gacATTTTGCAACCCGTCAGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTAAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTa                                                                                    >  8:327903/3‑149 (MQ=255)
                                                            gCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAg                                                                                                                                >  6:310351/1‑105 (MQ=255)
                                                            gCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAg                                                                                                                                <  5:310351/105‑1 (MQ=255)
                                                            gCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTa                                                                                    >  8:231649/1‑149 (MQ=255)
                                                             ccATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTAc                                                                                   <  5:359206/149‑1 (MQ=255)
                                                               aTTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCt                                                                                 >  3:180001/1‑149 (MQ=255)
                                                                ttttGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCt                                                                                 >  1:401027/1‑148 (MQ=255)
                                                                           gTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTg                                                                     >  6:314961/1‑149 (MQ=255)
                                                                               ggCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTaa                                                                 >  1:349189/1‑149 (MQ=255)
                                                                                gCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTa                                                                  <  7:231649/147‑1 (MQ=255)
                                                                                 cAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACg                                                               >  5:359903/1‑149 (MQ=255)
                                                                                  aGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCt                                                                                                                   <  2:78044/96‑1 (MQ=255)
                                                                                  aGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCt                                                                                                                   >  1:78044/1‑96 (MQ=255)
                                                                                    gCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGgtg                                                            >  5:331278/1‑149 (MQ=255)
                                                                                     cGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGgtgt                                                           <  7:294595/149‑1 (MQ=255)
                                                                                        tCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCaa                                                                          >  1:116717/1‑131 (MQ=255)
                                                                                        tCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCaa                                                                          <  2:116717/131‑1 (MQ=255)
                                                                                         cgcgGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGGGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTaa                                                        <  6:130239/148‑1 (MQ=255)
                                                                                          gcgGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAAt                                                       >  1:353650/1‑148 (MQ=255)
                                                                                          gcgGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGT‑GTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATg                                                      >  2:286669/1‑148 (MQ=255)
                                                                                          gcgGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGT‑GTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGa                                                     >  2:8607/1‑149 (MQ=255)
                                                                                           aggAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGa                                                     <  5:314961/148‑1 (MQ=255)
                                                                                            ggAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGAt                                                    <  1:293441/149‑1 (MQ=255)
                                                                                               gaTGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCct                                                 <  8:167875/149‑1 (MQ=255)
                                                                                                aTGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCctc                                                <  4:106635/149‑1 (MQ=255)
                                                                                                aTGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCct                                                 >  6:161203/1‑148 (MQ=255)
                                                                                                 tGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGAt                                                    >  5:214936/1‑144 (MQ=255)
                                                                                                 tGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGGGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGAt                                                    <  6:214936/144‑1 (MQ=255)
                                                                                                  gccgccGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTaa                                                        <  8:149942/139‑1 (MQ=255)
                                                                                                  gccgccGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTaa                                                        >  7:149942/1‑139 (MQ=255)
                                                                                                    cgccgTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCt                                            >  5:95117/1‑149 (MQ=255)
                                                                                                       cgTGCTGCTGGCGGCTCTGGCTTTCCGGGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGcc                                         <  4:47960/149‑1 (MQ=255)
                                                                                                        gTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCt                                        <  6:331278/149‑1 (MQ=255)
                                                                                                            tgctgGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTc                                    >  4:202507/1‑149 (MQ=255)
                                                                                                              ctgGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTc                                    <  5:161203/147‑1 (MQ=255)
                                                                                                                   ggcACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTg                             <  2:353650/149‑1 (MQ=255)
                                                                                                                   ggcACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTg                             <  3:202507/149‑1 (MQ=255)
                                                                                                                   ggcACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGGGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTg                             <  1:165442/149‑1 (MQ=255)
                                                                                                                    gcACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTgg                            >  5:266108/1‑149 (MQ=255)
                                                                                                                     cacaGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGc                           >  8:293070/1‑149 (MQ=255)
                                                                                                                      acaGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGcc                          >  3:229996/1‑149 (MQ=255)
                                                                                                                       caGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGccc                         >  1:413681/1‑149 (MQ=255)
                                                                                                                         ggCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGt                       >  7:246162/1‑149 (MQ=255)
                                                                                                                            tATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATACTCTCCTTCCTTGTCCATCTTGGCCCGTTGc                    >  4:233579/1‑149 (MQ=255)
                                                                                                                              tCCGCGCGGTGGCGCTGGGTGT‑GTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCg                  <  1:286669/148‑1 (MQ=255)
                                                                                                                              tCCGCGCGGTGGCGCTGGGTGT‑GTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCg                  <  1:8607/148‑1 (MQ=255)
                                                                                                                               ccGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTTCCgg                 >  5:149605/1‑149 (MQ=255)
                                                                                                                                cgcgcgGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATc                                                   <  7:147184/114‑1 (MQ=255)
                                                                                                                                cgcgcgGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATc                                                   >  8:147184/1‑114 (MQ=255)
                                                                                                                                cgcgcgGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCgg                 >  6:173029/1‑148 (MQ=255)
                                                                                                                                  cgcgGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTAc              <  2:125105/149‑1 (MQ=255)
                                                                                                                                              gggTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGa  <  6:95117/149‑1 (MQ=255)
                                                                                                                                               ggTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGa  <  3:234617/148‑1 (MQ=255)
                                                                                                                                                gtgtGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTAc              >  6:40714/1‑135 (MQ=255)
                                                                                                                                                gtgtGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTAc              <  5:40714/135‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
TGATGCTTCTCTTCAGTGCCCTGCTGTACTGGCGCGGTGAACAGGTGGCACAAGGCATTCGCCATTTTGCAACCCGTCTGGCAGGCGTTCGCGGAGATGCCGCCGTGCTGCTGGCGGCACAGGCTATCCGCGCGGTGGCGCTGGGTGTGGTGGTGACGGCGTTAGTACAGGCAGTGCTTGGCGGTATCGGCCTCGCCGTATCCGGCGTACCTTATGCAACGTTGCTAACGGTGTTAATGATCCTCTCCTGCCTTGTCCAGCTTGGCCCGTTGCCGGTACTGATTCCGGCGA  >  NC_000913/1769579‑1769869

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