Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,790,303 A→T 10.1% intergenic (‑222/+17) waaH ← / ← tdh UDP‑glucuronate:LPS(HepIII) glycosyltransferase/threonine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,790,3030AT10.1% 338.9 / 23.1 119intergenic (‑222/+17)waaH/tdhUDP‑glucuronate:LPS(HepIII) glycosyltransferase/threonine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (50/57);  new base T (7/5);  total (57/62)
Fisher's exact test for biased strand distribution p-value = 5.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.62e-01

Only 100 of 121 total aligned reads displayed.

GATATTATTTTTACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGC  >  NC_000913/3790154‑3790451
                                                                                                                                                     |                                                                                                                                                    
tatatTATTTTTACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTAGATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATa                                                                                                                                                      <  2:949186/149‑1 (MQ=255)
 atatTATTTTTACAGGTGAGAGAGACAATTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAAATACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATACTCAGATAAAATTGCTGGAATaa                                                                                                                                                     >  2:601220/1‑149 (MQ=255)
  tattatTTTTACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATAcc                                                                                                                                                    >  1:341473/1‑150 (MQ=255)
   attattTTTACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCa                                                                                                                                                   >  1:202999/1‑150 (MQ=255)
   attattTTTACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCa                                                                                                                                                   >  2:575627/1‑150 (MQ=255)
     tattTTTACAGTTGAGAGAGTCTCTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGc                                                                                                                                                 <  2:592337/150‑1 (MQ=255)
     tattTTTACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAAGTGCGCACTCTTTATTACACTGAGAGGATTTTACAGTTTTACCAGATAATCCTCAGATAAAAGGGATGGAATACCAGc                                                                                                                                                 >  2:649182/1‑150 (MQ=255)
      attTTTACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGcc                                                                                                                                                >  1:897508/1‑150 (MQ=255)
        ttttACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCt                                                                                                                                              <  2:957770/150‑1 (MQ=255)
        ttttACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCt                                                                                                                                              <  1:820546/150‑1 (MQ=255)
        ttttACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCt                                                                                                                                              <  2:202999/150‑1 (MQ=255)
          ttACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttt                                                                                                                                            <  1:405197/150‑2 (MQ=255)
          ttACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttt                                                                                                                                            <  2:803112/150‑2 (MQ=255)
          ttACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTg                                                                                                                                            <  1:614552/150‑1 (MQ=255)
          ttACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTg                                                                                                                                            <  1:155483/150‑1 (MQ=255)
             cAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttat                                                                                                                                         >  1:1572340/1‑146 (MQ=255)
             cAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTc                                                                                                                                         >  2:179282/1‑150 (MQ=255)
             cAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTc                                                                                                                                         >  1:895328/1‑150 (MQ=255)
              aGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatc                                                                                                                                        <  2:44432/150‑6 (MQ=255)
               ggTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgt                                                                                                                                       >  1:535478/1‑150 (MQ=255)
                gTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtg                                                                                                                                      <  2:1043017/150‑1 (MQ=255)
                gTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtg                                                                                                                                      <  1:1065724/150‑1 (MQ=255)
                 tGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtgt                                                                                                                                     <  2:287393/150‑1 (MQ=255)
                 tGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtgt                                                                                                                                     <  2:1371990/150‑1 (MQ=255)
                 tGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtgt                                                                                                                                     <  2:1373791/150‑1 (MQ=255)
                        gTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATccc                                                                                                                              >  1:1577988/1‑150 (MQ=255)
                         tCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCa                                                                                                                             >  2:705312/1‑150 (MQ=255)
                               cAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAg                                                                                                                       <  2:323956/150‑1 (MQ=255)
                                  ggTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAAt                                                                                                                    <  1:433736/150‑1 (MQ=255)
                                     agagagATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAAc                                                                                                                 <  2:510564/150‑1 (MQ=255)
                                     agagagATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAAc                                                                                                                 >  1:1342114/1‑150 (MQ=255)
                                     agagagATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAAc                                                                                                                 <  1:709172/150‑1 (MQ=255)
                                      gagagaTAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACt                                                                                                                >  1:578443/1‑150 (MQ=255)
                                           aTAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTccc                                                                                                           >  1:713299/1‑150 (MQ=255)
                                            tAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCg                                                                                                          <  2:1055584/150‑1 (MQ=255)
                                            tAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCg                                                                                                          <  2:1050141/150‑1 (MQ=255)
                                             aaTGAGACCAGCTATATAAGAGATCAGCACAATTACCAACTCTATATTACACTGACAGGATTTTACATATTTAACACATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCgg                                                                                                         <  2:1136861/150‑1 (MQ=255)
                                               tGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGAc                                                                                                       <  2:1136824/150‑1 (MQ=255)
                                                  gagCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTgg                                                                                                    <  1:439536/150‑1 (MQ=255)
                                                    gcagcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGcc                                                                                                  <  2:1568038/150‑1 (MQ=255)
                                                      agcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCg                                                                                                <  2:347366/150‑1 (MQ=255)
                                                         tatataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAAc                                                                                             >  2:795850/1‑150 (MQ=255)
                                                         tatataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTAAATATTTAACAGATTATCCTCAGATAAAAGTGCTGGAATACCAGCCCTTTTTCGTGTTAATAAAAGCTAATAATAACTTTCACGGACTGGACCGaaa                                                                                             >  1:445815/1‑149 (MQ=255)
                                                          atattAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACg                                                                                            <  1:1043199/150‑1 (MQ=255)
                                                           tataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGc                                                                                           <  2:1299769/150‑1 (MQ=255)
                                                           tataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGc                                                                                           <  1:268404/150‑1 (MQ=255)
                                                           tataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGc                                                                                           >  2:1114801/1‑150 (MQ=255)
                                                               agagATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAg                                                                                       >  1:1397247/1‑150 (MQ=255)
                                                                   aTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTc                                                                                   <  1:1387573/150‑1 (MQ=255)
                                                                   aTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTc                                                                                   <  1:705312/150‑1 (MQ=255)
                                                                      aGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCaaa                                                                                >  2:1602967/1‑150 (MQ=255)
                                                                          caATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGccc                                                                            >  1:1220085/1‑150 (MQ=255)
                                                                             ttACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTc                                                                         <  1:575627/150‑1 (MQ=255)
                                                                                cGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTgg                                                                      <  1:834997/150‑1 (MQ=255)
                                                                                 gAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGa                                                                     <  2:1584325/150‑1 (MQ=255)
                                                                                 gAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGa                                                                     <  1:1467501/150‑1 (MQ=255)
                                                                                    atctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAat                                                                  <  2:1570938/149‑1 (MQ=255)
                                                                                         tattaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaaa                                                                   >  2:1069585/1‑144 (MQ=255)
                                                                                         tattaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaaa                                                                   <  1:1069585/144‑1 (MQ=255)
                                                                                         tattaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCg                                                             >  1:1209756/1‑150 (MQ=255)
                                                                                          attaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGa                                                            >  1:542804/1‑150 (MQ=255)
                                                                                          attaCACTGACAGGAGTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGa                                                            >  2:139793/1‑150 (MQ=255)
                                                                                              cacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGTGTCAAAGCCCTTCTGGAAATCATCGATaga                                                        <  2:1315183/150‑1 (MQ=255)
                                                                                              cacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATaga                                                        <  2:753032/150‑1 (MQ=255)
                                                                                              cacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATaga                                                        <  1:1114801/150‑1 (MQ=255)
                                                                                              cacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATaga                                                        >  1:199808/1‑150 (MQ=255)
                                                                                                     aGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGa                                                 >  1:108651/1‑150 (MQ=255)
                                                                                                              cATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATg                                        <  1:179282/150‑1 (MQ=255)
                                                                                                               atatTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGa                                       <  2:535478/150‑1 (MQ=255)
                                                                                                                tatTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGAt                                      >  1:1370135/1‑150 (MQ=255)
                                                                                                                   ttAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCgg                                   >  2:494035/1‑150 (MQ=255)
                                                                                                                    tAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGc                                  >  1:4269/1‑150 (MQ=255)
                                                                                                                        aGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGt                                            <  1:282716/136‑1 (MQ=255)
                                                                                                                        aGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGt                                            >  2:282716/1‑136 (MQ=255)
                                                                                                                         taTTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAg                             <  2:1566736/149‑1 (MQ=255)
                                                                                                                         gATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAg                             <  2:1382118/150‑1 (MQ=255)
                                                                                                                         gATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAg                             >  1:1350754/1‑150 (MQ=255)
                                                                                                                                cTCATATAAAATGGATGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAgg                      <  2:895328/150‑1 (MQ=255)
                                                                                                                                cTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAg                             >  1:1568973/1‑143 (MQ=255)
                                                                                                                                cTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAg                             <  2:1568973/143‑1 (MQ=255)
                                                                                                                                cTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAgg                      <  2:199808/150‑1 (MQ=255)
                                                                                                                                cTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAgg                      >  1:735186/1‑150 (MQ=255)
                                                                                                                                cTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAgg                      <  2:578443/150‑1 (MQ=255)
                                                                                                                                      aacaaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGAc                >  2:854454/4‑150 (MQ=255)
                                                                                                                                       acaaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACt               >  1:248636/3‑150 (MQ=255)
                                                                                                                                       aaaaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACt               <  2:676409/150‑1 (MQ=255)
                                                                                                                                       aaaaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACt               <  1:1272372/150‑1 (MQ=255)
                                                                                                                                        caaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTg              >  2:324555/2‑150 (MQ=255)
                                                                                                                                        caaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTg              >  1:77584/2‑150 (MQ=255)
                                                                                                                                        aaaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTg              <  1:243294/150‑1 (MQ=255)
                                                                                                                                        aaaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACACATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACtc              >  1:199131/1‑149 (MQ=255)
                                                                                                                                         aaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGa             >  1:1060406/1‑150 (MQ=255)
                                                                                                                                         aaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGa             >  2:1136903/1‑150 (MQ=255)
                                                                                                                                         aaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGa             >  1:794786/1‑150 (MQ=255)
                                                                                                                                          aGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGaa            >  2:1447031/1‑150 (MQ=255)
                                                                                                                                            ggCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATc          <  2:166552/150‑1 (MQ=255)
                                                                                                                                             gCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCa         >  2:1330731/1‑150 (MQ=255)
                                                                                                                                             gCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCa         >  2:1327980/1‑150 (MQ=255)
                                                                                                                                               tGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAgc       <  1:711300/150‑1 (MQ=255)
                                                                                                                                                    tACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGcgccgc  >  1:1151122/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GATATTATTTTTACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGC  >  NC_000913/3790154‑3790451

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