Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,798,234 A→T 6.6% intergenic (+27/+5) waaL → / ← waaU O‑antigen ligase/lipopolysaccharide core biosynthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,798,2340AT 6.6% 3416.7 / 161.3 960intergenic (+27/+5)waaL/waaUO‑antigen ligase/lipopolysaccharide core biosynthesis
Reads supporting (aligned to +/- strand):  ref base A (458/438);  new base T (32/31);  total (491/469)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-01

Only 100 of 993 total aligned reads displayed.

AATTGAAGCAGGGTCACTGAAAGGTCTGATGGGAATTTTTTCCACACTTTTCCTCTATTTTTCACTATTTTATATAGCATATAAAAAACGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTATGATAAGCCGCAGCGATATGAACAAGAGCTGTATCAACAGAAATGACAAAATCACTATACTTGGTCAA  >  NC_000913/3797970‑3798487
                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                             
                                                                               tataAAAAACGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaaaa                                                                                                                                                                                                                                                            >  3:261937/1‑189 (MQ=255)
                                                                                                                      ggCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaaa                                                                                                                                                                                                                                                             <  7:832884/149‑1 (MQ=255)
                                                                                                                       gCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaaaa                                                                                                                                                                                                                                                            <  6:403870/149‑1 (MQ=255)
                                                                                                                            gtgGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCAc                                                                                                                                                                                                                                                       <  8:781738/149‑1 (MQ=255)
                                                                                                                             tgGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCAct                                                                                                                                                                                                                                                      <  5:824055/149‑1 (MQ=255)
                                                                                                                             tgGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCAct                                                                                                                                                                                                                                                      >  8:1254511/1‑149 (MQ=255)
                                                                                                                               gggATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActtc                                                                                                                                                                                                                                                    <  6:1215360/149‑1 (MQ=255)
                                                                                                                                     ggACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCt                                                                                                                                                                                                                                              <  6:110615/149‑1 (MQ=255)
                                                                                                                                         tCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTATAAGAAGTGAGTTTTAACTCACTTCTTaaacttgt                                                                                                                                                                                                                                          >  7:1405881/1‑142 (MQ=39)
                                                                                                                                           aGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACaa                                                                                                                                                                                                                                        >  7:407332/1‑149 (MQ=255)
                                                                                                                                            gtgatgtgatCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAAc                                                                                                                                                                                                                                       >  8:1559950/1‑149 (MQ=255)
                                                                                                                                                  tgatCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAAtc                                                                                                                                                                                                                                 <  5:192401/149‑1 (MQ=255)
                                                                                                                                                    atCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAAtctc                                                                                                                                                                                                                               <  6:129053/149‑1 (MQ=255)
                                                                                                                                                     tCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTatg                                                                                                                                                  >  1:10613/1‑225 (MQ=255)
                                                                                                                                                      cATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCtt                                                                                                                                                                                                                             <  6:1033105/149‑1 (MQ=255)
                                                                                                                                                       atatGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaa                                                                                                                                                                                                                                                              >  4:443469/1‑115 (MQ=255)
                                                                                                                                                       atatGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCttt                                                                                                                                                                                                                            >  7:570146/1‑149 (MQ=255)
                                                                                                                                                       atatGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCttt                                                                                                                                                                                                                            <  8:407332/149‑1 (MQ=255)
                                                                                                                                                       atatGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTaaa                                                                                                                            >  2:466740/1‑245 (MQ=255)
                                                                                                                                                       atatGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCg                                                                                                                        <  1:8183/249‑1 (MQ=255)
                                                                                                                                                                 gcagcaTTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAg                                                                                                                                                                                                                  >  7:1170224/1‑149 (MQ=255)
                                                                                                                                                                  cagcaTTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGt                                                                                                                                                                                                                 <  5:1370272/149‑1 (MQ=255)
                                                                                                                                                                    gcaTTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTa                                                                                                                                                                                                               >  6:155521/1‑149 (MQ=255)
                                                                                                                                                                    gcaTTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTa                                                                                                                                                                                                               >  5:247048/1‑149 (MQ=255)
                                                                                                                                                                     caTTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTaa                                                                                                                                                                                                              >  5:94000/1‑149 (MQ=255)
                                                                                                                                                                      aTTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAAt                                                                                                                                                                                                                     <  7:147544/141‑1 (MQ=255)
                                                                                                                                                                       ttCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAg                                                                                                                                                                                                            <  7:1028941/149‑1 (MQ=255)
                                                                                                                                                                       ttCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAg                                                                                                                                                                                                            >  8:266934/1‑149 (MQ=255)
                                                                                                                                                                           aaTTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTAc                                                                                                                                                                                                                          >  6:303361/1‑131 (MQ=255)
                                                                                                                                                                           aaTTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTc                                                                                                                                                                                                        <  7:331413/149‑1 (MQ=255)
                                                                                                                                                                           aaTTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTc                                                                                                                                                                                                        <  7:968996/149‑1 (MQ=255)
                                                                                                                                                                            aTTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCa                                                                                                                                                                                                       <  8:943045/149‑1 (MQ=255)
                                                                                                                                                                                  aTTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCg                                                                                                                                                                                                 >  5:1455163/1‑149 (MQ=255)
                                                                                                                                                                                    tataTCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGat                                                                                                                                                                                               <  5:1142696/149‑1 (MQ=255)
                                                                                                                                                                                      tatCCGCTATAGTCCTCCTACTCGTCATTAATAATCGTAACAATACAATTAGTTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGatat                                                                                                                                                                                             <  6:1370975/149‑1 (MQ=255)
                                                                                                                                                                                       atCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATc                                                                                                                                                                                            <  6:893441/149‑1 (MQ=255)
                                                                                                                                                                                        tCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCt                                                                                                                                                                                           >  6:704465/1‑149 (MQ=255)
                                                                                                                                                                                          cGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCttt                                                                                                                                                                                         >  7:642501/1‑149 (MQ=255)
                                                                                                                                                                                            cTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTa                                                                                                                                                                                       <  8:1170224/149‑1 (MQ=255)
                                                                                                                                                                                                  tCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGtat                                                                                                                                                                                 >  7:1520673/1‑149 (MQ=255)
                                                                                                                                                                                                  tCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGtat                                                                                                                                                                                 <  5:367272/149‑1 (MQ=255)
                                                                                                                                                                                                   ccTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGCTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGtata                                                                                                                                                                                <  6:63970/149‑1 (MQ=255)
                                                                                                                                                                                                     tctTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATg                                                                                                                                                                              >  5:969675/1‑149 (MQ=255)
                                                                                                                                                                                                      ctTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAttt                                                                                                                                                               >  2:717891/1‑163 (MQ=255)
                                                                                                                                                                                                       ttACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTaaa                                                                                                                                                                                                             >  8:1590498/1‑116 (MQ=255)
                                                                                                                                                                                                         aCTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAAc                                                                                                                                                                                                                                                           >  6:74617/1‑68 (MQ=255)
                                                                                                                                                                                                         aCTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTgg                                                                                                                                                                          >  5:1062233/1‑149 (MQ=255)
                                                                                                                                                                                                         aCTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTgg                                                                                                                                                                          <  6:1616328/149‑1 (MQ=255)
                                                                                                                                                                                                          cTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAAt                                                                                                                                                                                                                     >  6:932149/1‑105 (MQ=255)
                                                                                                                                                                                                          cTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTggg                                                                                                                                                                         <  5:666660/149‑1 (MQ=255)
                                                                                                                                                                                                           tCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTgggg                                                                                                                                                                        >  5:1670573/1‑149 (MQ=255)
                                                                                                                                                                                                             gTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGAt                                                                                                                                                                      >  5:904611/1‑149 (MQ=255)
                                                                                                                                                                                                              tCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGtt                                                                                                                                                                                                                <  5:293012/106‑1 (MQ=255)
                                                                                                                                                                                                              tCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTggg                                                                                                                                                                         <  7:209499/145‑1 (MQ=255)
                                                                                                                                                                                                              tCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATa                                                                                                                                                                     >  8:1111382/1‑149 (MQ=255)
                                                                                                                                                                                                               cATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATAc                                                                                                                                                                    <  5:715723/149‑1 (MQ=255)
                                                                                                                                                                                                                aTTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCtt                                                                                                                                                                                                                             >  5:1208276/1‑91 (MQ=255)
                                                                                                                                                                                                                aTTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACt                                                                                                                                                                   >  5:1543438/1‑149 (MQ=255)
                                                                                                                                                                                                                 ttAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTa                                                                                                                                                                  >  8:214688/1‑149 (MQ=255)
                                                                                                                                                                                                                   aataatCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAtt                                                                                                                                                                >  7:1126001/1‑149 (MQ=255)
                                                                                                                                                                                                                   aataatCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAtt                                                                                                                                                                >  8:480286/1‑149 (MQ=255)
                                                                                                                                                                                                                    ataatCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAttt                                                                                                                                                               <  6:978626/149‑1 (MQ=255)
                                                                                                                                                                                                                    ataatCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAttt                                                                                                                                                               <  8:667133/149‑1 (MQ=255)
                                                                                                                                                                                                                     taatCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAttt                                                                                                                                                               <  7:259693/148‑1 (MQ=255)
                                                                                                                                                                                                                      aatCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCa                                                                                                                                                                                                       <  8:322137/107‑1 (MQ=255)
                                                                                                                                                                                                                         cGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGtata                                                                                                                                                          <  8:142333/149‑1 (MQ=255)
                                                                                                                                                                                                                         cGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGtata                                                                                                                                                          <  6:1164779/149‑1 (MQ=255)
                                                                                                                                                                                                                            aaCAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGAt                                                                                                                                                       >  7:688445/1‑149 (MQ=255)
                                                                                                                                                                                                                               aatacaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTa                                                                                                                                                    >  8:8982/1‑149 (MQ=255)
                                                                                                                                                                                                                                atacaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTCTTTGTATAGATTTat                                                                                                                                                   >  5:1286522/1‑149 (MQ=255)
                                                                                                                                                                                                                                  acaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTatga                                                                                                                                                 >  7:1219071/1‑149 (MQ=255)
                                                                                                                                                                                                                                  acaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTatga                                                                                                                                                 <  7:130760/149‑1 (MQ=255)
                                                                                                                                                                                                                                         attaaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTa                                                                                                                                                                                       >  5:466078/1‑104 (MQ=255)
                                                                                                                                                                                                                                         attaaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCg                                                                                                                                          >  8:1499135/1‑149 (MQ=255)
                                                                                                                                                                                                                                             aGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGa                                                                                                                                      >  5:1608995/1‑149 (MQ=255)
                                                                                                                                                                                                                                             aGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGa                                                                                                                                      >  7:922267/1‑149 (MQ=255)
                                                                                                                                                                                                                                              aaaTAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGAc                                                                                                                                     <  8:1663723/148‑1 (MQ=255)
                                                                                                                                                                                                                                                    aCAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTg                                                                                                                                                                           >  8:543041/1‑105 (MQ=255)
                                                                                                                                                                                                                                                    aCAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGta                                                                                                                                                            >  5:551702/1‑120 (MQ=255)
                                                                                                                                                                                                                                                    aCAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTat                                                                                                                                                   >  7:513615/1‑129 (MQ=255)
                                                                                                                                                                                                                                                    aCAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAAtt                                                                                                                               >  7:472729/1‑149 (MQ=255)
                                                                                                                                                                                                                                                       aGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTaaa                                                                                                                            <  7:1584786/149‑1 (MQ=255)
                                                                                                                                                                                                                                                        gTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAAt                                                                                                                           <  5:429628/149‑1 (MQ=255)
                                                                                                                                                                                                                                                         tttAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATg                                                                                                                          >  5:766044/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         tttAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATg                                                                                                                          >  8:828416/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         tttAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATg                                                                                                                          >  6:359278/1‑149 (MQ=255)
                                                                                                                                                                                                                                                            aagaagTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGa                                                                                                                       <  6:194170/149‑1 (MQ=255)
                                                                                                                                                                                                                                                              gaagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAgg                                                                                                                     >  6:166907/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               aagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAAATAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCggg                                                                                                                                        >  6:1469922/1‑129 (MQ=255)
                                                                                                                                                                                                                                                               aagTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAggg                                                                                                                    <  7:165343/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                 gTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTTGGGATACTATTTGTATAGATTTa                                                                                                                                                    <  8:104921/115‑1 (MQ=255)
                                                                                                                                                                                                                                                                 gTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAttt                                                                                                                                                               <  6:795726/104‑1 (MQ=255)
                                                                                                                                                                                                                                                                  tGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATc                                                                                                                                                                                            >  4:721982/1‑74 (MQ=255)
                                                                                                                                                                                                                                                                   gAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAg                                                                                                                      <  8:725510/143‑1 (MQ=255)
                                                                                                                                                                                                                                                                   gAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAg                                                                                                                      >  7:725510/1‑143 (MQ=255)
                                                                                                                                                                                                                                                                     gTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATc                                                                                                                                                                                                  <  8:605321/65‑1 (MQ=255)
                                                                                                                                                                                                                                                                     gTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATTc                                                                                                              >  8:1328427/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                      ttAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAttt                                                                                                                                                               >  6:497564/1‑99 (MQ=255)
                                                                                                                                                                                                                                                                      ttAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTa                                                                                                                                                    >  7:770455/1‑110 (MQ=255)
                                                                                                                                                                                                                                                                        aaaaCTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGt                                                                                                                                                                                   <  6:255562/77‑1 (MQ=255)
                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                             
AATTGAAGCAGGGTCACTGAAAGGTCTGATGGGAATTTTTTCCACACTTTTCCTCTATTTTTCACTATTTTATATAGCATATAAAAAACGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTATGATAAGCCGCAGCGATATGAACAAGAGCTGTATCAACAGAAATGACAAAATCACTATACTTGGTCAA  >  NC_000913/3797970‑3798487

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: