Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,798,235 A→T 5.4% 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 5.4% 1154.8 / 37.8 318intergenic (+28/+4)waaL/waaUO‑antigen ligase/lipopolysaccharide core biosynthesis
Reads supporting (aligned to +/- strand):  ref base A (158/142);  new base T (11/6);  total (170/148)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.68e-01

Only 100 of 332 total aligned reads displayed.

TGAGATAATTGAAGCAGGGTCACTGAAAGGTCTGATGGGAATTTTTTCCACACTTTTCCTCTATTTTTCACTATTTTATATAGCATATAAAAAACGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTATGATAAGCCGCAGCGATATGAACAAGAGCTGTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAA  >  NC_000913/3797964‑3798501
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                          
tGAGATAATTGAAGCAGGGTCACTGAAAGGTCTGATGGGAATTTTTTCCACACTTTTCCTCTATTTTTCACTATTTTATATAGCATATAAAAAACGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAAAATAAACAAGTTTAATAATTGATTTAAAAc                                                                                                                                                                                                                                                                         >  4:609152/1‑275 (MQ=255)
                                     ggAATTTTTTCCACACTTTTCCTCTATTTTTCACTATTTTATATAGCATATAAAAAACGAGCTTTTGGTTTGTTGATATTAACGCTTGGCATTTTGGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAGTTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGaa                                                                                                                                                                                                                                    <  2:363458/275‑1 (MQ=255)
                                                                    cACTATTTTATATAGCATATAAAAAACGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTa                                                                                                                                                                                                     <  2:609154/275‑1 (MQ=255)
                                                                                   cATATAAAAAACGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTgggg                                                                                                                                                                                      >  2:10492/1‑275 (MQ=255)
                                                                                                         tttGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaa                                                                                                                                                                                                                                                                            >  4:362345/1‑167 (MQ=255)
                                                                                                                             gCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaaaa                                                                                                                                                                                                                                                                          <  8:888747/149‑1 (MQ=255)
                                                                                                                             gCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaaa                                                                                                                                                                                                                                                                           <  6:401211/148‑1 (MQ=255)
                                                                                                                             ccATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCAGTATTTCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTaaa                                                                                                                                                                                                                                                                           <  6:655346/147‑1 (MQ=255)
                                                                                                                                 tgtgGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTc                                                                                                                                                                                                                                                                       <  5:876942/148‑1 (MQ=255)
                                                                                                                                  gtgGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCAc                                                                                                                                                                                                                                                                     <  8:113382/149‑1 (MQ=255)
                                                                                                                                   tgGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCAct                                                                                                                                                                                                                                                                    >  8:234658/1‑149 (MQ=255)
                                                                                                                                     gggATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCActtc                                                                                                                                                                                                                                                                  <  8:390754/149‑1 (MQ=255)
                                                                                                                                         ttGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTa                                                                                                                                                                                                                                                               >  6:515823/1‑148 (MQ=255)
                                                                                                                                              cTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTTTAACTCAc                                                                                                                                                                                                                                                                     >  5:824164/1‑137 (MQ=255)
                                                                                                                                              cTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAAc                                                                                                                                                                                                                                                                         <  6:746221/133‑1 (MQ=255)
                                                                                                                                               tCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAAc                                                                                                                                                                                                                                                                         >  7:620209/1‑132 (MQ=255)
                                                                                                                                               tCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATAc                                                                                                                                                                                                                                                        >  7:421237/1‑149 (MQ=255)
                                                                                                                                                cAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTAt                                                                                                                                                                                                                                                              <  5:646034/142‑1 (MQ=255)
                                                                                                                                                cAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATAc                                                                                                                                                                                                                                                        >  6:126011/1‑148 (MQ=255)
                                                                                                                                                        tgatCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTaa                                                                                                                                                                                                                                                 >  5:377704/1‑147 (MQ=255)
                                                                                                                                                         gatCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATc                                                                                                                                                                                                                                                             >  7:557674/1‑134 (MQ=255)
                                                                                                                                                         gatCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAAtc                                                                                                                                                                                                                                               >  6:72653/1‑148 (MQ=255)
                                                                                                                                                         gatCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAAt                                                                                                                                                                                                                                                >  7:739312/1‑147 (MQ=255)
                                                                                                                                                             atatGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTc                                                                                                                                                                                                                                                                       >  3:49185/1‑120 (MQ=255)
                                                                                                                                                             atatGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCttt                                                                                                                                                                                                                                          >  7:722115/1‑149 (MQ=255)
                                                                                                                                                              tatGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActtctt                                                                                                                                                                                                                                                                <  5:593230/126‑1 (MQ=255)
                                                                                                                                                              tatGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCttt                                                                                                                                                                                                                                          >  6:189104/1‑148 (MQ=255)
                                                                                                                                                              tatGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTa                                                                                                                                                                                                                                         <  8:590216/149‑1 (MQ=255)
                                                                                                                                                                   gCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGa                                                                                                                                                                                                                                     >  6:893246/1‑148 (MQ=255)
                                                                                                                                                                        cagcaTTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAg                                                                                                                                                                                                                                >  5:669748/1‑148 (MQ=255)
                                                                                                                                                                                 aaTTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTaa                                                                                                                                                                                                                                                 <  6:782285/122‑1 (MQ=255)
                                                                                                                                                                                 aaTTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTaa                                                                                                                                                                                                                                                 >  5:782285/1‑122 (MQ=255)
                                                                                                                                                                                  aTTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTc                                                                                                                                                                                                                      <  5:126011/148‑1 (MQ=255)
                                                                                                                                                                                   ttATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAg                                                                                                                                                                                                                    >  7:602867/1‑149 (MQ=255)
                                                                                                                                                                                     aTCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCa                                                                                                                                                                                                                     <  7:813336/146‑1 (MQ=255)
                                                                                                                                                                                     aTCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCa                                                                                                                                                                                                                     >  8:813336/1‑146 (MQ=255)
                                                                                                                                                                                     aTCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAAGTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTa                                                                                                                                                                                                                  >  8:482081/1‑149 (MQ=255)
                                                                                                                                                                                              tCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAg                                                                                                                                                                                                                                >  7:336792/1‑126 (MQ=255)
                                                                                                                                                                                                  cTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTa                                                                                                                                                                                                     <  8:702217/149‑1 (MQ=255)
                                                                                                                                                                                                   tataGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTACAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTAAGTATCGATATCTTTTAc                                                                                                                                                                                                    >  8:534544/1‑149 (MQ=255)
                                                                                                                                                                                                   tataGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActtc                                                                                                                                                                                                                                                                  >  2:746778/1‑87 (MQ=255)
                                                                                                                                                                                                     taGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAAtt                                                                                                                                                                                                                                  >  6:538804/1‑117 (MQ=255)
                                                                                                                                                                                                        tCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAg                                                                                                                                                                                                                    <  8:603683/128‑1 (MQ=255)
                                                                                                                                                                                                        tCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGta                                                                                                                                                                                                >  6:133212/1‑148 (MQ=255)
                                                                                                                                                                                                          ctctTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTatat                                                                                                                                                                                             >  8:30807/1‑149 (MQ=255)
                                                                                                                                                                                                          ctctTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTGTCAGTATCGATATCTTTTACCTTatat                                                                                                                                                                                             >  8:781973/1‑149 (MQ=255)
                                                                                                                                                                                                             ttACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTc                                                                                                                                                                                                                                                                       >  4:549053/1‑72 (MQ=255)
                                                                                                                                                                                                              tACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGCGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGtt                                                                                                                                                                                          >  5:744038/1‑148 (MQ=255)
                                                                                                                                                                                                               aCTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTgg                                                                                                                                                                                        <  7:232403/149‑1 (MQ=255)
                                                                                                                                                                                                                cTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTggg                                                                                                                                                                                       <  8:210826/149‑1 (MQ=255)
                                                                                                                                                                                                                 tCGTCATTAATAATCGTAACGATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCActt                                                                                                                                                                                                                                                                   <  2:549055/72‑1 (MQ=255)
                                                                                                                                                                                                                 tCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTgggg                                                                                                                                                                                      >  7:853660/1‑149 (MQ=255)
                                                                                                                                                                                                                    tCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGAt                                                                                                                                                                                    >  6:563916/1‑148 (MQ=255)
                                                                                                                                                                                                                       ttAATAATCGTAACAATACAATTAATAAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTa                                                                                                                                                                                >  7:57970/1‑149 (MQ=255)
                                                                                                                                                                                                                        taataaTCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAt                                                                                                                                                                               <  8:722115/149‑1 (MQ=255)
                                                                                                                                                                                                                         aataatCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAt                                                                                                                                                                               <  5:588956/148‑1 (MQ=255)
                                                                                                                                                                                                                         aataatCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAt                                                                                                                                                                               <  6:742812/148‑1 (MQ=255)
                                                                                                                                                                                                                          ataatCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAttt                                                                                                                                                                             <  7:181757/149‑1 (MQ=255)
                                                                                                                                                                                                                             atCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGta                                                                                                                                                                          <  7:92272/149‑1 (MQ=255)
                                                                                                                                                                                                                               cGTAACAATACAATTTATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTAACGTATATGTTGGGGATACTATTTGTATAGATTTATCATGATTCGGTTACCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTGGGGTAAATTGCAAGCGTTGGTTTATTATAACCCGAAGCGATATGAACAAGAGATGTATAAACAGaa                                          >  4:671854/1‑275 (MQ=255)
                                                                                                                                                                                                                                gTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTCGTATAg                                                                                                                                                                       >  7:163408/1‑149 (MQ=255)
                                                                                                                                                                                                                                      atacaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTCTTTGTATAGATTTat                                                                                                                                                                 <  8:668702/149‑1 (MQ=255)
                                                                                                                                                                                                                                      atacaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTa                                                                                                                                                                  >  5:276503/1‑148 (MQ=255)
                                                                                                                                                                                                                                        acaatTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTatg                                                                                                                                                                <  6:235014/148‑1 (MQ=255)
                                                                                                                                                                                                                                           attaattaAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTAc                                                                                                                                                                                                    >  5:343485/1‑109 (MQ=255)
                                                                                                                                                                                                                                           attaattaAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTatgatg                                                                                                                                                             <  6:690090/148‑1 (MQ=255)
                                                                                                                                                                                                                                             taattaaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATgatgat                                                                                                                                                           <  5:239759/148‑1 (MQ=255)
                                                                                                                                                                                                                                                ttaaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCgg                                                                                                                                                       <  8:172871/149‑1 (MQ=255)
                                                                                                                                                                                                                                                 taaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCgg                                                                                                                                                       >  6:60667/1‑148 (MQ=255)
                                                                                                                                                                                                                                                 taaGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCg                                                                                                                                                        >  6:177906/1‑147 (MQ=255)
                                                                                                                                                                                                                                                     aaTAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAg                                                                                                                                                                       <  5:419275/128‑1 (MQ=255)
                                                                                                                                                                                                                                                     aaTAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAg                                                                                                                                                                       >  6:419275/1‑128 (MQ=255)
                                                                                                                                                                                                                                                      aTAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTc                                                                                                                                                                                                                      <  6:814517/80‑1 (MQ=255)
                                                                                                                                                                                                                                                      aTAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTc                                                                                                                                                                                                                      >  5:814517/1‑80 (MQ=255)
                                                                                                                                                                                                                                                       tAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCa                                                                                                                                                 >  6:552019/1‑148 (MQ=255)
                                                                                                                                                                                                                                                        aaaCAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCaa                                                                                                                                                >  5:909430/1‑148 (MQ=255)
                                                                                                                                                                                                                                                         aaCAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCaaa                                                                                                                                               >  5:182773/1‑148 (MQ=255)
                                                                                                                                                                                                                                                               tttAAGAAGTGAGTTAAAACTCTCTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCgg                                                                                                                                                       <  7:442638/134‑1 (MQ=255)
                                                                                                                                                                                                                                                               tttAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTg                                                                                                                                                                                         <  7:296621/100‑1 (MQ=255)
                                                                                                                                                                                                                                                                ttAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCg                                                                                                                                                        <  6:374766/132‑1 (MQ=255)
                                                                                                                                                                                                                                                                ttAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATg                                                                                                                                        >  5:533548/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                  aagaagTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCg                                                                                                                                      <  5:742470/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                  aagaagTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGa                                                                                                                                     <  8:437238/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                   agaagTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGatat                                                                                                                              >  1:170326/1‑155 (MQ=255)
                                                                                                                                                                                                                                                                    gaagTGAGTTTTAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAg                                                                                                                                    >  5:555089/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                    gaagTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAg                                                                                                                                    <  5:543021/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aagTGAGTTAAAACTCACTTCTTATCTATACAACTTAGTCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAgg                                                                                                                                   >  6:282064/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                          aGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTAtt                                                                                                                                                                              >  1:36273/1‑100 (MQ=255)
                                                                                                                                                                                                                                                                            ttAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGt                                                                                                                                                                                                 <  3:174296/79‑1 (MQ=255)
                                                                                                                                                                                                                                                                            ttAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCggg                                                                                                                                                      <  8:28331/122‑1 (MQ=255)
                                                                                                                                                                                                                                                                            ttAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCggg                                                                                                                                                      >  7:28331/1‑122 (MQ=255)
                                                                                                                                                                                                                                                                            ttAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCggg                                                                                                                                                      >  6:318832/1‑122 (MQ=255)
                                                                                                                                                                                                                                                                            ttAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCggg                                                                                                                                                      <  5:318832/122‑1 (MQ=255)
                                                                                                                                                                                                                                                                            ttAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATTcc                                                                                                                           >  7:681600/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                             tAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATTcc                                                                                                                           <  6:268268/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                             tAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATTcc                                                                                                                           <  5:512773/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                               aaaCTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATTCCgg                                                                                                                         <  5:705524/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                               aaaCTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATTCCgg                                                                                                                         <  5:149344/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                               aaaCTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATTCCgg                                                                                                                         <  5:133212/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                               aaaCTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATTCCGGa                                                                                                                        >  8:352114/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                          
TGAGATAATTGAAGCAGGGTCACTGAAAGGTCTGATGGGAATTTTTTCCACACTTTTCCTCTATTTTTCACTATTTTATATAGCATATAAAAAACGAGCTTTGGGTTTGTTGATATTAACGCTTGGCATTGTGGGGATTGGACTCAGTGATGTGATCATATGGGCACGCAGCATTCCAATTATCATTATATCCGCTATAGTCCTCTTACTCGTCATTAATAATCGTAACAATACAATTAATTAAGAATAAACAAGTTTAAGAAGTGAGTTAAAACTCACTTCTTATCTATACAACTTAATCTCTTTACTGAATTAGTTAAAGTTTCAGTATCGATATCTTTTACCGTATATGTTGGGGATACTATTTGTATAGATTTATGATGATTCGGGGACCAAATTAAATGCGAGGGATATTCCGGAGTACGTGAATTAGGGTAAAATGCAAGCGTTGGTTTATGATAAGCCGCAGCGATATGAACAAGAGCTGTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAA  >  NC_000913/3797964‑3798501

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