Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,790,303 A→T 5.4% intergenic (‑222/+17) waaH ← / ← tdh LPS(HepIII)‑glucuronic acid glycosyltransferase/L‑threonine 3‑dehydrogenase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,790,3030AT 5.4% 4012.2 / 142.1 1114intergenic (‑222/+17)waaH/tdhLPS(HepIII)‑glucuronic acid glycosyltransferase/L‑threonine 3‑dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (542/510);  new base T (28/33);  total (570/544)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.99e-02

Only 100 of 1121 total aligned reads displayed.

TATAACGGAATAATAACACTAAGTTTATTGGTGCTGTTCATCATAAATTACCGACGACCTCTGAAGTTATCCTGACGTTTTACAGGCTTTGTATTAAGAAACTATTAAGCCTGTGAGATATTATTTTTACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGG  >  NC_000913/3790038‑3790575
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                
   aaCGGAATATTTACTCTAAGCTTTTTGGTACTTTTCATCATAAATTCCCGTCGCCCTCTGTAGTTATCCTGTCGTTGTACAGGCCTTGTAGTAAGAAACAATTAAGCCGTTGATCTTTTACTTTTACATGTGAGCGCCTGACTTCAGGGTAGCGAGATAATGGGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTATATTACTCTGTCATGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATTCCAGCCCTTGtg                                                                                                                                                                                                                                                                      <  4:255744/275‑2 (MQ=255)
               acacTAAGTTTATTGGTGCTGTTCATCATAAATTACCGACGACCTCTGAAGTTATCCTGACGTTTTACAGGCTTTGTATTAAGAAACTATTAAGCCTGTGAGATATTATTTTGACAGTTTAGAGAGTCACTTCAGGGTAGATAGACAATGAGAGCAGCTATATAAGAGATAAGCACAATTACGAACTCTTTATTACACTGACAGGTATTTACAGATTTAACAGATTACCCTCAGATCAAAGGACCGGAAAACCAGCCCTTTCTCTTTTTAATcgc                                                                                                                                                                                                                                                          >  4:551910/1‑273 (MQ=255)
                                                                        tGACGTTTTACAGGCTTTGTATTAAGAAACTATTAAGCCTGTGAGATATTATTTTTACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAg                                                                                                                                                                                                                                                        <  4:471176/220‑1 (MQ=255)
                                                                                                             ccTGTGAGATATTATTTTTTCAGGGGGGAGAGTAACTTAAGGGTAGGGGGATAATGAGAGCAGCTATATAAGAGATAAGCACAATTACGAACACTTTATTACAATGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGa                                                                                                                                                            <  4:165731/275‑1 (MQ=255)
                                                                                                                                  aGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTCCGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTc                                                                                                                                                                                                                                                                     <  8:956275/149‑1 (MQ=255)
                                                                                                                                      gagagaGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtgt                                                                                                                                                                                                                                                                 <  8:1376067/149‑1 (MQ=255)
                                                                                                                                          gagTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAAt                                                                                                                                                                                                                                                             >  7:1154093/1‑149 (MQ=255)
                                                                                                                                              cACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCa                                                                                                                                                                                                                                                         >  5:621878/1‑149 (MQ=255)
                                                                                                                                                 ttCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCt                                                                                                                                                                                                                                                      >  6:1714672/1‑149 (MQ=255)
                                                                                                                                                   cAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCa                                                                                                                                                                                                                                                    <  5:456229/149‑1 (MQ=255)
                                                                                                                                                      ggTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGaa                                                                                                                                                                                                                                                 <  7:1052875/149‑1 (MQ=255)
                                                                                                                                                      ggTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGaa                                                                                                                                                                                                                                                 <  6:871568/149‑1 (MQ=255)
                                                                                                                                                                  aTGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCgg                                                                                                                                                                                                                                     >  7:794735/1‑149 (MQ=255)
                                                                                                                                                                   tGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGa                                                                                                                                                                                                                                    >  7:715828/1‑149 (MQ=255)
                                                                                                                                                                    gagaGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatc                                                                                                                                                                                                                                                                    <  7:1343494/116‑6 (MQ=255)
                                                                                                                                                                     agagCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACt                                                                                                                                                                                                                                  >  6:650375/1‑149 (MQ=255)
                                                                                                                                                                      gagCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTg                                                                                                                                                                                                                                 >  6:785593/1‑149 (MQ=255)
                                                                                                                                                                          agcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGccc                                                                                                                                                                                                                             >  5:940979/1‑149 (MQ=255)
                                                                                                                                                                             tatataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGaa                                                                                                                                                                                                                          <  6:1485592/149‑1 (MQ=255)
                                                                                                                                                                             tagataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAAAAGATTATCCGCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCAGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACg                                                                                                <  2:667090/268‑1 (MQ=255)
                                                                                                                                                                              atataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAAc                                                                                                                                                                                                                         <  8:160531/149‑1 (MQ=255)
                                                                                                                                                                                  aaGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCAt                                                                                                                                                                                                                     >  5:1626582/1‑149 (MQ=255)
                                                                                                                                                                                  aaGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGATAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCAt                                                                                                                                                                                                                     >  7:1443116/1‑149 (MQ=255)
                                                                                                                                                                                   agagATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCg                                                                                                                                                                                                                                      <  5:1157024/131‑1 (MQ=255)
                                                                                                                                                                                    gagaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGc                                                                                                                                                                                                                               <  8:1303285/137‑1 (MQ=255)
                                                                                                                                                                                       aTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtgt                                                                                                                                                                                                                                                                 >  8:620746/1‑100 (MQ=255)
                                                                                                                                                                                       aTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTTTACCAGGTTTCAAACATCCCACGCCCGTAAATAccc                                                                                   >  4:64759/1‑273 (MQ=255)
                                                                                                                                                                                         cAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTa                                                                                >  2:527099/1‑275 (MQ=255)
                                                                                                                                                                                         cAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCAAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCATAAATACCTTTAAtca                                                                            >  1:93214/1‑277 (MQ=255)
                                                                                                                                                                                          aGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCaa                                                                                                                                                                                                             >  5:223065/1‑149 (MQ=255)
                                                                                                                                                                                           gCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCaaa                                                                                                                                                                                                            <  5:1251190/149‑1 (MQ=255)
                                                                                                                                                                                                aTTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACa                                                                         <  2:325169/275‑1 (MQ=255)
                                                                                                                                                                                                 ttACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGccc                                                                                                                                                                                                                                                                           >  1:286182/1‑80 (MQ=255)
                                                                                                                                                                                                 ttACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCtt                                                                                                                                                                                                      >  5:292561/1‑149 (MQ=255)
                                                                                                                                                                                                   aCGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCt                                                                                                                                                                                                    <  7:219225/149‑1 (MQ=255)
                                                                                                                                                                                                     gAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATaa                                                                                                                                                                                                                                              >  6:1332385/1‑105 (MQ=255)
                                                                                                                                                                                                        ctctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaaa                                                                                                                                                                                               <  7:1534094/149‑1 (MQ=255)
                                                                                                                                                                                                         tctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAat                                                                                                                                                                                              >  6:589082/1‑149 (MQ=255)
                                                                                                                                                                                                           tttATTACACTGACCGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCaaa                                                                                                                                                                                                            >  7:119864/1‑133 (MQ=255)
                                                                                                                                                                                                           tttATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtgt                                                                                                                                                                                                                                                                 >  5:1023501/1‑80 (MQ=255)
                                                                                                                                                                                                            ttattaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatctt                                                                                                                                                                                           <  8:1157422/149‑3 (MQ=255)
                                                                                                                                                                                                            ttattaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcat                                                                                                                                                                                           <  5:1687637/149‑1 (MQ=255)
                                                                                                                                                                                                               ttaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGaa                                                                                                                                                                                                                                                 <  6:1550115/92‑1 (MQ=255)
                                                                                                                                                                                                               ttaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGAc                                                                                                                                                                                                                                   <  4:409740/106‑1 (MQ=255)
                                                                                                                                                                                                                   acTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATAc                                                                                                                                                                                                                                                                                 >  4:772389/1‑56 (MQ=255)
                                                                                                                                                                                                                   acTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATaga                                                                                                                                                                                    >  8:754767/1‑149 (MQ=255)
                                                                                                                                                                                                                    tnaaaagaaTTATACATATTTAACAGTTTAACCNNAGATAAAAGGGCTGGAATACCGGCCCTTTTNCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGAATGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCANCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTcc                                                 <  3:112593/271‑1 (MQ=255)
                                                                                                                                                                                                                     tGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGcc                                                                                                                                                                                                         <  1:110299/126‑1 (MQ=255)
                                                                                                                                                                                                                     tGACAGGATTTTACATATTTAACAGATTATCCGCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTTGAAat                                                                                                                                                                                              <  1:147965/137‑1 (MQ=255)
                                                                                                                                                                                                                       aCAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaaa                                                                                                                                                                                               >  2:750421/1‑134 (MQ=255)
                                                                                                                                                                                                                       aCAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaa                                                                                                                                                                                                >  7:979949/1‑133 (MQ=255)
                                                                                                                                                                                                                       aCAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcatc                                                                                                                                                                                          >  5:75745/1‑139 (MQ=255)
                                                                                                                                                                                                                       aCAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGaaa                                                                                                                                                                                <  6:701780/149‑1 (MQ=255)
                                                                                                                                                                                                                       aCAGGATTTTAAATATTTAACAGATTAGCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTc                                                  <  2:173241/275‑1 (MQ=255)
                                                                                                                                                                                                                        cAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAAc                                                                                                                                                                               >  7:278578/1‑149 (MQ=255)
                                                                                                                                                                                                                         aGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTAcc                                                                                                                   <  4:346860/208‑1 (MQ=255)
                                                                                                                                                                                                                           gATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGAt                                                                                                                                                                            <  7:428711/149‑1 (MQ=255)
                                                                                                                                                                                                                            aTTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATg                                                                                                                                                                           <  6:432792/149‑1 (MQ=255)
                                                                                                                                                                                                                             ttttACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATgg                                                                                                                                                                          <  6:231616/149‑1 (MQ=255)
                                                                                                                                                                                                                                  cATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAAt                                                                                                                                                                     >  7:571317/1‑149 (MQ=255)
                                                                                                                                                                                                                                   atatTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTc                                                                                                                                                                                                                                         >  8:345969/1‑80 (MQ=255)
                                                                                                                                                                                                                                   atatTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATg                                                                                                                                                                    >  6:1432138/1‑149 (MQ=255)
                                                                                                                                                                                                                                   atatTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATg                                                                                                                                                                    >  7:138652/1‑149 (MQ=255)
                                                                                                                                                                                                                                     atTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCCTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCCGCGCCGCCATCTTGTACCAGGTTTCa                                                                                                           <  2:248827/204‑1 (MQ=255)
                                                                                                                                                                                                                                      tttAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATc                                                                                                                                                                 >  7:1276273/1‑149 (MQ=255)
                                                                                                                                                                                                                                         aaCAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAat                                                                                                                                                                                              <  8:872209/117‑1 (MQ=255)
                                                                                                                                                                                                                                          aCAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCg                                                                                                                                                             >  8:700994/1‑149 (MQ=255)
                                                                                                                                                                                                                                          aCAGATTATACTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcatc                                                                                                                                                                                          >  5:1477000/1‑120 (MQ=255)
                                                                                                                                                                                                                                             gATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGaaa                                                                                                                                                          <  6:1279503/149‑1 (MQ=255)
                                                                                                                                                                                                                                             gATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGaaa                                                                                                                                                          >  5:862855/1‑149 (MQ=255)
                                                                                                                                                                                                                                              aTTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTc                                                                                                                                                                                                                                         >  6:108795/1‑69 (MQ=255)
                                                                                                                                                                                                                                              aTTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTc                                                                                                                                                                                                                                         <  5:108795/69‑1 (MQ=255)
                                                                                                                                                                                                                                                 aTCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATaa                                                                                                                                                                                                                                              <  3:88309/61‑1 (MQ=255)
                                                                                                                                                                                                                                                   ccTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGa                                                                                                                                                    <  6:1012513/149‑1 (MQ=255)
                                                                                                                                                                                                                                                   ccTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGa                                                                                                                                                    <  6:1186491/149‑1 (MQ=255)
                                                                                                                                                                                                                                                   ccTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGa                                                                                                                                                    <  8:284842/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    cTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAg                                                                                                                                                   <  6:504972/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    cTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAg                                                                                                                                                   <  5:1368258/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    cTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAg                                                                                                                                                   >  5:677440/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     tCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAAc                                                                                                                                                                                                                                             >  2:45053/1‑58 (MQ=255)
                                                                                                                                                                                                                                                     tCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaa                                                                                                                                                                                                <  7:1651767/103‑1 (MQ=255)
                                                                                                                                                                                                                                                     tCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGa                                                                                                                                                                                        <  8:1051104/111‑1 (MQ=255)
                                                                                                                                                                                                                                                     tCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGa                                                                                                                                                    <  5:540772/147‑1 (MQ=255)
                                                                                                                                                                                                                                                       aGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGcc                                                                                                                                                <  7:901994/149‑1 (MQ=255)
                                                                                                                                                                                                                                                        gATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCt                                                                                                                                                                                                    <  8:21231/96‑1 (MQ=255)
                                                                                                                                                                                                                                                         gaacaaGGGCTGGTATTCCAGCCCTTTTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAg                                                                                                                                              >  8:1418660/5‑149 (MQ=255)
                                                                                                                                                                                                                                                           aaaaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGAc                                                                                                                                            <  5:1547719/149‑1 (MQ=255)
                                                                                                                                                                                                                                                           aaaaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGAc                                                                                                                                            <  6:1646922/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                 gCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATc                                                                                                                                      <  7:651280/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                 gCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATg                                                                             >  2:582378/1‑206 (MQ=255)
                                                                                                                                                                                                                                                                   tGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATgg                                                                                                                                                                          <  6:982111/111‑1 (MQ=255)
                                                                                                                                                                                                                                                                   tGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTaaa                                                                                        >  3:225691/1‑193 (MQ=255)
                                                                                                                                                                                                                                                                   tGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTgg                                                    >  2:377415/1‑229 (MQ=255)
                                                                                                                                                                                                                                                                   tGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGGCGATAGACATATCAAACGGCGGAATACCCAGACTCGCAATACgg  >  1:217494/1‑279 (MQ=255)
                                                                                                                                                                                                                                                                    ggAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAgc                                                                                                                                   <  8:456359/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                    ggAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAgc                                                                                                                                   <  7:350793/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                     gAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAgcg                                                                                                                                  >  8:125035/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                      aaTACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAgcgc                                                                                                                                 <  8:1513567/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                         aCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGa                                                                                                                                                    <  6:34951/127‑1 (MQ=255)
                                                                                                                                                                                                                                                                         aCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGa                                                                                                                                                    >  5:34951/1‑127 (MQ=255)
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                
TATAACGGAATAATAACACTAAGTTTATTGGTGCTGTTCATCATAAATTACCGACGACCTCTGAAGTTATCCTGACGTTTTACAGGCTTTGTATTAAGAAACTATTAAGCCTGTGAGATATTATTTTTACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAATACGG  >  NC_000913/3790038‑3790575

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