Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,798,235 A→T 6.7% intergenic (+28/+4) 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,2350AT 6.7% 3397.5 / 162.5 960intergenic (+28/+4)waaL/waaUO‑antigen ligase/lipopolysaccharide core biosynthesis
Reads supporting (aligned to +/- strand):  ref base A (458/438);  new base T (33/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 = 2.40e-01

Only 100 of 1002 total aligned reads displayed.

AATTGAAGCAGGGTCACTGAAAGGTCTGATGGGAATTTTTTCCACACTTTTCCTCTATTTTTCACTATTTTATATAGCATATAAAAAACGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTATGATAAGCCGCAGCGATATGAACAAGAGCTGTATCAACAGAAATGACAAAATCACTATACTTGGTCAA  >  NC_000913/3797970‑3798487
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                            
                                                                             cATATAAAAAACGAGCTTTGGGTTTGTTGATATTAACGCTTGGCCTTGTGGGGCTTGGACTCAGTGATGTGATCATATGGGAACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATAGTAACAATACACTTAATTAATAATAAAGAAGTTTAAAAAGTGAGTTAAAACCCAATTTTTATCTATACAAATTAATCTCTTTAATGAATTAATTAAAGTTTCAATTACGAAATATTTTACCTGATATGTTgggg                                                                                                                                                                        >  4:445517/1‑275 (MQ=255)
                                                                                                      gttgATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGATATAGTCATCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATTATGATTc                                                                                                                                           >  1:83089/1‑279 (MQ=255)
                                                                                                                       gCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaaaa                                                                                                                                                                                                                                                            <  8:319009/149‑1 (MQ=255)
                                                                                                                           tgtgGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCa                                                                                                                                                                                                                                                        <  8:1255149/149‑1 (MQ=255)
                                                                                                                           tgtgGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCa                                                                                                                                                                                                                                                        <  7:728493/149‑1 (MQ=255)
                                                                                                                               gggATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCActtc                                                                                                                                                                                                                                                    <  8:939609/149‑1 (MQ=255)
                                                                                                                                        cTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCtata                                                                                                                                                                                                                                           <  6:692280/149‑1 (MQ=255)
                                                                                                                                         tCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATAc                                                                                                                                                                                                                                          >  5:1375647/1‑149 (MQ=255)
                                                                                                                                           aGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACaa                                                                                                                                                                                                                                        >  7:407332/1‑149 (MQ=255)
                                                                                                                                            gtgatgtgatCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAAc                                                                                                                                                                                                                                       >  8:1559950/1‑149 (MQ=255)
                                                                                                                                                  tgatCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAAtc                                                                                                                                                                                                                                 <  6:628541/149‑1 (MQ=255)
                                                                                                                                                  tgatCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAAtc                                                                                                                                                                                                                                 <  5:192401/149‑1 (MQ=255)
                                                                                                                                                     tCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACaa                                                                                                                                                                                                                                        <  8:1050910/139‑1 (MQ=255)
                                                                                                                                                     tCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTatg                                                                                                                                                  >  1:10613/1‑225 (MQ=255)
                                                                                                                                                      aatatGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCtt                                                                                                                                                                                                                             <  6:299185/148‑1 (MQ=255)
                                                                                                                                                       atatGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCttt                                                                                                                                                                                                                            <  8:407332/149‑1 (MQ=255)
                                                                                                                                                       atatGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCttt                                                                                                                                                                                                                            >  7:570146/1‑149 (MQ=255)
                                                                                                                                                        tatGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTa                                                                                                                                                                                                                           <  8:775710/149‑1 (MQ=255)
                                                                                                                                                        tatGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTa                                                                                                                                                                                                                           >  6:745584/1‑149 (MQ=255)
                                                                                                                                                           gggCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTg                                                                                                                                                                                                                        >  6:400986/1‑149 (MQ=255)
                                                                                                                                                                cGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTa                                                                                                                                                                                                                   <  7:851690/149‑1 (MQ=255)
                                                                                                                                                                 gcagcaTTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAg                                                                                                                                                                                                                  >  7:1170224/1‑149 (MQ=255)
                                                                                                                                                                 gcagcaTTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAg                                                                                                                                                                                                                  <  8:82105/149‑1 (MQ=255)
                                                                                                                                                                    gcaTTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTa                                                                                                                                                                                                               >  5:247048/1‑149 (MQ=255)
                                                                                                                                                                       ttCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAg                                                                                                                                                                                                            >  8:266934/1‑149 (MQ=255)
                                                                                                                                                                       ttCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAg                                                                                                                                                                                                            <  7:1028941/149‑1 (MQ=255)
                                                                                                                                                                            aTTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCa                                                                                                                                                                                                       <  6:1477812/149‑1 (MQ=255)
                                                                                                                                                                                  aTTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCg                                                                                                                                                                                                 >  5:1455163/1‑149 (MQ=255)
                                                                                                                                                                                       atCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATc                                                                                                                                                                                            <  6:893441/149‑1 (MQ=255)
                                                                                                                                                                                          cGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCttt                                                                                                                                                                                         >  7:642501/1‑149 (MQ=255)
                                                                                                                                                                                               taGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCtttt                                                                                                                                                                                        <  8:1170767/145‑1 (MQ=255)
                                                                                                                                                                                                  tCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGtat                                                                                                                                                                                 <  5:367272/149‑1 (MQ=255)
                                                                                                                                                                                                      ctTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGAt                                                                                                                                                                      <  4:773422/156‑1 (MQ=255)
                                                                                                                                                                                                       ttACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTaaa                                                                                                                                                                                                             >  8:1590498/1‑116 (MQ=255)
                                                                                                                                                                                                       ttACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGtata                                                                                                                                                                                >  2:492487/1‑145 (MQ=255)
                                                                                                                                                                                                        tACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTg                                                                                                                                                                           >  7:1412106/1‑149 (MQ=255)
                                                                                                                                                                                                         aCTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTgg                                                                                                                                                                          <  5:1143655/149‑1 (MQ=255)
                                                                                                                                                                                                          cTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAAt                                                                                                                                                                                                                     >  6:932149/1‑105 (MQ=255)
                                                                                                                                                                                                          cTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTggg                                                                                                                                                                         <  8:718301/149‑1 (MQ=255)
                                                                                                                                                                                                          cTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTggg                                                                                                                                                                         <  7:97366/149‑1 (MQ=255)
                                                                                                                                                                                                          cTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTggg                                                                                                                                                                         <  5:666660/149‑1 (MQ=255)
                                                                                                                                                                                                           tCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCa                                                                                                                                                                                                       >  8:1115452/1‑118 (MQ=255)
                                                                                                                                                                                                           tCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTgggg                                                                                                                                                                        >  7:206532/1‑149 (MQ=255)
                                                                                                                                                                                                           tCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTgggg                                                                                                                                                                        >  5:1670573/1‑149 (MQ=255)
                                                                                                                                                                                                            cGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACt                                                                                                                                                                                                                                                          <  3:36110/66‑1 (MQ=255)
                                                                                                                                                                                                             gTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGCAt                                                                                                                                                                      >  7:1488144/1‑149 (MQ=255)
                                                                                                                                                                                                              tCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGtt                                                                                                                                                                                                                <  5:293012/106‑1 (MQ=255)
                                                                                                                                                                                                              tCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTaaa                                                                                                                                                                                                             <  6:790952/109‑1 (MQ=255)
                                                                                                                                                                                                               cATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATAc                                                                                                                                                                    <  5:715723/149‑1 (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)
                                                                                                                                                                                                                      aatCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGt                                                                                                                                                             <  7:774292/149‑1 (MQ=255)
                                                                                                                                                                                                                        tCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGt                                                                                                                                                             >  6:1633358/1‑147 (MQ=255)
                                                                                                                                                                                                                         cGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGtata                                                                                                                                                          <  8:142333/149‑1 (MQ=255)
                                                                                                                                                                                                                          gTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAg                                                                                                                                                         >  5:905004/1‑149 (MQ=255)
                                                                                                                                                                                                                           tAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTgg                                                                                                                                                                          >  6:1678469/1‑131 (MQ=255)
                                                                                                                                                                                                                           tAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAttt                                                                                                                                                               <  6:1381049/142‑1 (MQ=255)
                                                                                                                                                                                                                            aaCAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGAt                                                                                                                                                       >  7:688445/1‑149 (MQ=255)
                                                                                                                                                                                                                               aatacaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTa                                                                                                                                                    <  6:1322406/149‑1 (MQ=255)
                                                                                                                                                                                                                                atacaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTat                                                                                                                                                   <  5:944480/149‑1 (MQ=255)
                                                                                                                                                                                                                                  acaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTatga                                                                                                                                                 <  7:130760/149‑1 (MQ=255)
                                                                                                                                                                                                                                  acaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTatga                                                                                                                                                 >  7:1219071/1‑149 (MQ=255)
                                                                                                                                                                                                                                  acaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTatga                                                                                                                                                 <  6:1559978/149‑1 (MQ=255)
                                                                                                                                                                                                                                   caatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTatgat                                                                                                                                                <  5:1208877/149‑1 (MQ=255)
                                                                                                                                                                                                                                    aattaattAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGt                                                                                                                                                                                                                 >  5:979288/1‑83 (MQ=255)
                                                                                                                                                                                                                                    aattaattAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTatgatg                                                                                                                                               >  5:1084738/1‑149 (MQ=255)
                                                                                                                                                                                                                                    aattaattAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTatgatg                                                                                                                                               >  6:1228519/1‑149 (MQ=255)
                                                                                                                                                                                                                                    aattaattAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTatgatg                                                                                                                                               >  6:1044729/1‑149 (MQ=255)
                                                                                                                                                                                                                                        aattaaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTc                                                                                                                                           <  6:1123754/149‑1 (MQ=255)
                                                                                                                                                                                                                                         attaaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTa                                                                                                                                                                                       >  5:466078/1‑104 (MQ=255)
                                                                                                                                                                                                                                         attaaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCg                                                                                                                                          >  8:1499135/1‑149 (MQ=255)
                                                                                                                                                                                                                                              aaaTAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGAc                                                                                                                                     <  8:1663723/148‑1 (MQ=255)
                                                                                                                                                                                                                                               aaTAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGAcc                                                                                                                                    >  7:922981/1‑149 (MQ=255)
                                                                                                                                                                                                                                                  aaaCAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCaaa                                                                                                                                 <  6:34530/149‑1 (MQ=255)
                                                                                                                                                                                                                                                   aaCAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAAt                                                                                                                                <  5:1601005/149‑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                                                                                                                            <  6:1181311/149‑1 (MQ=255)
                                                                                                                                                                                                                                                       aGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTaaa                                                                                                                            <  7:1584786/149‑1 (MQ=255)
                                                                                                                                                                                                                                                        gTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAAt                                                                                                                           <  5:429628/149‑1 (MQ=255)
                                                                                                                                                                                                                                                         tttAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATg                                                                                                                          >  5:766044/1‑149 (MQ=255)
                                                                                                                                                                                                                                                            aagaagTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGa                                                                                                                       <  6:194170/149‑1 (MQ=255)
                                                                                                                                                                                                                                                            aagaagTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGa                                                                                                                       <  6:252261/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               aagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTat                                                                                                                                                   >  5:840580/1‑118 (MQ=255)
                                                                                                                                                                                                                                                               aagTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCaaa                                                                                                                                 <  6:466346/136‑1 (MQ=255)
                                                                                                                                                                                                                                                                agTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGt                                                                                                                                                                                                     <  8:893772/67‑1 (MQ=255)
                                                                                                                                                                                                                                                                  tGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGata                                                                                                                 <  7:1522992/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                     gTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATc                                                                                                                                                                                                  <  8:605321/65‑1 (MQ=255)
                                                                                                                                                                                                                                                                     gTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATTc                                                                                                              >  8:1328427/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                      ttAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCtttt                                                                                                                                                                                        <  2:721987/74‑1 (MQ=255)
                                                                                                                                                                                                                                                                         aaaCTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAAt                                                                                                                           <  5:1294238/132‑1 (MQ=255)
                                                                                                                                                                                                                                                                         aaaCTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATTCCGGa                                                                                                          >  7:324307/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         aaaCTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATTCCGGa                                                                                                          <  7:1647428/149‑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: