Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,249,437 G→T 5.6% intergenic (+126/‑117) lamB → / → malM maltose outer membrane channel/phage lambda receptor protein/maltose regulon periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,249,4370GT 5.6% 718.6 / 20.6 215intergenic (+126/‑117)lamB/malMmaltose outer membrane channel/phage lambda receptor protein/maltose regulon periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (91/110);  new base T (6/6);  total (98/117)
Fisher's exact test for biased strand distribution p-value = 7.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.44e-01

Only 100 of 231 total aligned reads displayed.

GCCCAGATGGAAATCTGGTGGTAATAGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTG  >  NC_000913/4249288‑4249585
                                                                                                                                                     |                                                                                                                                                    
tcccAGATGGAAATCTGGTGGTAATAGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTAcgccg                                                                                                                                                      >  2:1075253/2‑150 (MQ=255)
 cccAGATGGAAATCTGGTGGTAATAGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTAcgccgc                                                                                                                                                     >  2:1897377/1‑150 (MQ=255)
 cccAGATGGAAATCTGGTGGTAATAGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTAcgccgc                                                                                                                                                     >  2:382091/1‑150 (MQ=255)
   cAGATGGAAATCTGGTGGTAATAGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCAt                                                                                                                                                   <  2:324341/150‑1 (MQ=255)
     gATGGAAATCTGGTGGTAATAGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATcc                                                                                                                                                 <  2:1199955/150‑1 (MQ=255)
     gATGGAAATCTGGTGGTAATAGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATcc                                                                                                                                                 >  2:1400251/1‑150 (MQ=255)
        ggAAATCTGGTGGTAATAGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                              <  1:562876/150‑1 (MQ=255)
                ggtggtAATAGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTCGGCCGGATGCGGCGTAAACGCCTTATCCGGCCTacattt                                                                                                                                      <  2:926691/150‑6 (MQ=21)
                ggtggtAATAGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAc                                                                                                                                      >  1:1628391/1‑150 (MQ=255)
                ggtggtAATAGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCATATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                             >  1:151071/1‑142 (MQ=25)
                   ggtAATAGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAgg                                                                                                                                   >  1:1297387/1‑150 (MQ=255)
                   ggtAATAGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAgg                                                                                                                                   >  1:287113/1‑150 (MQ=255)
                    gtAATAGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGt                                                                                                                                  <  1:1141482/150‑1 (MQ=255)
                    gtAATAGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTACCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGt                                                                                                                                  >  1:1900683/1‑150 (MQ=255)
                       aTAGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGt                                                                                                                               <  2:1971508/150‑1 (MQ=255)
                         aGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCa                                                                                                                             <  1:1481430/150‑1 (MQ=255)
                         aGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCa                                                                                                                             >  2:1828806/1‑150 (MQ=255)
                         aGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCa                                                                                                                             <  2:1166454/150‑1 (MQ=255)
                          gCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCAt                                                                                                                            <  2:762036/150‑1 (MQ=255)
                          gCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCAt                                                                                                                            <  1:979175/150‑1 (MQ=255)
                          gCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCAt                                                                                                                            <  2:870047/150‑1 (MQ=255)
                            aaaaCCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTa                                                                                                                          <  1:35451/150‑1 (MQ=255)
                               aCCTTGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCt                                                                                                                       <  1:1122635/150‑1 (MQ=255)
                                ccTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTg                                                                                                                      >  2:1277332/1‑150 (MQ=255)
                                   gggCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGaaa                                                                                                                   >  1:1068676/1‑150 (MQ=255)
                                   gggCCGGATAAGGCGTTTACGCCGCATGCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGaaa                                                                                                                   >  1:875446/1‑150 (MQ=255)
                                     gCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAgg                                                                                                                 >  2:599387/1‑150 (MQ=255)
                                             aGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAg                                                                                                         <  2:1381219/150‑1 (MQ=255)
                                              ggCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGc                                                                                                        >  1:228583/1‑150 (MQ=255)
                                              ggCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGc                                                                                                        <  1:382091/150‑1 (MQ=255)
                                               gCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGcc                                                                                                       <  2:1453478/150‑1 (MQ=255)
                                                      cgccgcATTCGGCAACCGACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGa                                                                                                <  1:1809388/150‑1 (MQ=255)
                                                      cgccgcATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGa                                                                                                >  1:1660637/1‑150 (MQ=255)
                                                      cgccgcATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGa                                                                                                >  1:657794/1‑150 (MQ=255)
                                                       gccgcATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCGTTATCTGAAAGGGGCGAAAGCCCCTCTGAt                                                                                               >  1:1017549/1‑150 (MQ=255)
                                                        ccgcATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGAtt                                                                                              >  2:640250/1‑150 (MQ=255)
                                                          gcATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTAt                                                                                            >  2:1998435/1‑150 (MQ=255)
                                                           cATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATc                                                                                           >  1:1202456/1‑150 (MQ=255)
                                                              tCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCggg                                                                                        <  2:1628391/150‑1 (MQ=255)
                                                              tCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCggg                                                                                        <  2:1202456/150‑1 (MQ=255)
                                                               cGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGt                                                                                       >  2:1782140/1‑150 (MQ=255)
                                                               cGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGt                                                                                       <  2:895551/150‑1 (MQ=255)
                                                               cGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGt                                                                                       >  1:1325954/1‑150 (MQ=255)
                                                               cGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGt                                                                                       >  1:421261/1‑150 (MQ=255)
                                                                ggCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGtt                                                                                      >  2:671733/1‑150 (MQ=255)
                                                                 gCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGttt                                                                                     <  2:1248789/150‑1 (MQ=255)
                                                                  caaccaacGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTa                                                                                    >  1:229854/1‑150 (MQ=255)
                                                                   aaccaacGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAg                                                                                   <  2:1172517/150‑1 (MQ=255)
                                                                         cGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCt                                                                             <  1:1558063/150‑1 (MQ=255)
                                                                          gCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTa                                                                            >  2:1105126/1‑150 (MQ=255)
                                                                          gCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCACATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTa                                                                            >  1:1742752/1‑150 (MQ=255)
                                                                             tGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTg                                                                         <  2:481767/150‑1 (MQ=255)
                                                                                 gCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTg                                                                     >  1:1233953/1‑150 (MQ=255)
                                                                                  cGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAg                                                                                   >  1:816208/1‑135 (MQ=255)
                                                                                  cGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTgg                                                                    <  1:1732926/150‑1 (MQ=255)
                                                                                   gACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGc                                                                   >  2:1345695/1‑150 (MQ=255)
                                                                                   gACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGc                                                                   >  2:1238756/1‑150 (MQ=255)
                                                                                       cTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTAcc                                                               <  1:1405829/150‑1 (MQ=255)
                                                                                            gcgTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGa                                                          <  2:269062/150‑1 (MQ=255)
                                                                                             cgTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAg                                                         <  1:508774/150‑1 (MQ=255)
                                                                                                 ttATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTcc                                                     <  2:1286056/150‑1 (MQ=255)
                                                                                                  tATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTAGTCATTATCTGAAAGGGTCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCa                                                    >  2:1013130/1‑150 (MQ=255)
                                                                                                   aTCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAg                                                   <  2:1660637/150‑1 (MQ=255)
                                                                                                   aTCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAg                                                   <  2:1714564/150‑1 (MQ=255)
                                                                                                    tCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGa                                                  >  1:579302/1‑150 (MQ=255)
                                                                                                       ggCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGAttt                                               >  1:636087/1‑150 (MQ=255)
                                                                                                           tACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAg                                           <  1:1782140/150‑1 (MQ=255)
                                                                                                           tACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAg                                           >  2:1520836/1‑150 (MQ=255)
                                                                                                           tACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAg                                           <  2:1015029/150‑1 (MQ=255)
                                                                                                            aCAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAgg                                          <  1:1881629/150‑1 (MQ=255)
                                                                                                              aaCGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTg                                        <  1:671733/150‑1 (MQ=255)
                                                                                                                 ggCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGaaa                                     <  1:617075/150‑1 (MQ=255)
                                                                                                                    tGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACa                                  <  2:664721/150‑1 (MQ=255)
                                                                                                                    tGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACa                                  <  2:1760363/150‑1 (MQ=255)
                                                                                                                       cAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCt                                                                  >  2:251508/1‑115 (MQ=255)
                                                                                                                       cAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACaatg                               <  2:172783/150‑1 (MQ=255)
                                                                                                                             tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACaatgaaaatg                         <  2:1959800/150‑1 (MQ=255)
                                                                                                                              atgCCGGATGCGGCGTAAACGCCTTATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACaatgaaaatga                        >  2:655498/3‑150 (MQ=255)
                                                                                                                              aGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACaatgaaaatga                        <  1:1985489/150‑1 (MQ=255)
                                                                                                                               tgCCGGATGCGGCGTAAACGCCTTATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACaatgaaaatgaa                       >  2:1918429/2‑150 (MQ=255)
                                                                                                                               ggCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGTCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACaatgaaaatgaa                       >  1:1098835/1‑150 (MQ=255)
                                                                                                                               ggCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACaatgaaaatgaa                       >  1:907098/1‑150 (MQ=255)
                                                                                                                                gCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGTCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAAt                      <  2:1098835/150‑1 (MQ=255)
                                                                                                                                gCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAAt                      <  2:907098/150‑1 (MQ=255)
                                                                                                                                 ccGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATa                     >  2:1307187/1‑150 (MQ=255)
                                                                                                                                       aaGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGtc               >  2:461846/1‑150 (MQ=255)
                                                                                                                                         ggCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGtctc             >  2:1473783/1‑150 (MQ=255)
                                                                                                                                          gCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCa            >  2:800368/1‑150 (MQ=255)
                                                                                                                                           cgtaaaCGCCTTATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGa                                       >  1:1267312/6‑122 (MQ=255)
                                                                                                                                           cGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCAt           >  2:1145502/1‑150 (MQ=255)
                                                                                                                                             tttACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCg         <  1:1956714/150‑1 (MQ=255)
                                                                                                                                              ttACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGt        <  1:458429/150‑1 (MQ=255)
                                                                                                                                               tACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGAtt                                                                                              >  1:55602/1‑63 (MQ=255)
                                                                                                                                               tACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGAtt                                                                                              <  2:55602/63‑1 (MQ=255)
                                                                                                                                                  gccgcATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCctc    <  2:229854/150‑1 (MQ=255)
                                                                                                                                                  gccgcATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCctc    >  1:104611/1‑150 (MQ=255)
                                                                                                                                                   ccgcATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCt                                                                  >  1:1554884/1‑87 (MQ=255)
                                                                                                                                                   ccgcATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCt                                                                  <  2:1554884/87‑1 (MQ=255)
                                                                                                                                                   ccgcATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCctct   <  1:1008806/150‑1 (MQ=255)
                                                                                                                                                    cgcATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTg  <  1:1678435/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCCCAGATGGAAATCTGGTGGTAATAGCAAAACCTGGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACAACGGCTGTCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTG  >  NC_000913/4249288‑4249585

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