Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,790,303 A→T 5.8% 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.8% 2189.9 / 81.2 604intergenic (‑222/+17)waaH/tdhLPS(HepIII)‑glucuronic acid glycosyltransferase/L‑threonine 3‑dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (313/254);  new base T (21/14);  total (335/269)
Fisher's exact test for biased strand distribution p-value = 6.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.32e-02

Only 100 of 609 total aligned reads displayed.

TAACGGAATAATAACACTAAGTTTATTGGTGCTGTTCATCATAAATTACCGACGACCTCTGAAGTTATCCTGACGTTTTACAGGCTTTGTATTAAGAAACTATTAAGCCTGTGAGATATTATTTTTACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAAT  >  NC_000913/3790040‑3790571
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                            
        taataaCACTAAGTTTATTGGTGCTGTTCATCATAAATTACCGACGACCTCTGAAGTTATCCTGACGTTTTACAGGCTTTGTATTAAGAAACTATTAAGCCTGTGAGATATTATTTTTACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTa                                                                                                                                                                                                                                                           >  1:230071/1‑275 (MQ=255)
                                                                  aTCCTGACGTTTTACAGGCTTTGTATTAAGAAACTATTAAGCCTGTGAGGTATTATTTTTACAGGTGAGGGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTc                                                                                                                                                                                                 >  2:126652/1‑275 (MQ=255)
                                                                                                   cTATTAAGCCTGTGAGATATTATTTTTACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTccc                                                                                                                                                                                                                                   >  2:532710/1‑208 (MQ=255)
                                                                                                       tAAGCCTGTGAGATATTATTTTTACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGAc                                                                                                                                                                                                                               <  4:532708/208‑1 (MQ=255)
                                                                                                                  gaTATTATTTTTACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATcc                                                                                                                                                                                                                                                       <  3:108630/173‑1 (MQ=255)
                                                                                                                      ttattTTTACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATAcc                                                                                                                                                                                                                                                                            >  6:717867/1‑148 (MQ=255)
                                                                                                                           tttACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGccc                                                                                                                                                                                                                                                                       <  6:508250/148‑1 (MQ=255)
                                                                                                                           tttACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCt                                                                                                                                                                                                                                                                      <  7:149439/149‑1 (MQ=255)
                                                                                                                            ttACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCtt                                                                                                                                                                                                                                                                     <  8:598234/149‑1 (MQ=255)
                                                                                                                                    gagagaGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtg                                                                                                                                                                                                                                                              <  6:823409/148‑1 (MQ=255)
                                                                                                                                      gagagTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTa                                                                                                                                                                                                                                                           >  8:398564/1‑149 (MQ=255)
                                                                                                                                       agagTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGTCTGGTATTCCAg                                                                                                                                                                                                                                                                          >  5:305383/1‑133 (MQ=255)
                                                                                                                                       agagTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCa                                                                                                                                                                                                                                                >  2:716225/1‑159 (MQ=255)
                                                                                                                                        gagTCACCTCAGGGTAGAGAGATAGTGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACATTGACAGGACTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTaa                                                                                                                                                                                                                                                          >  5:216730/1‑148 (MQ=255)
                                                                                                                                           tCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATcc                                                                                                                                                                                                                                                       >  5:139633/1‑148 (MQ=255)
                                                                                                                                                 cAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCGCa                                                                                                                                                                                                                                                <  8:72984/149‑1 (MQ=255)
                                                                                                                                                    ggTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAAc                                                                                                                                                                                                                                         >  3:225840/1‑153 (MQ=255)
                                                                                                                                                       agagagATAATGAGAGCAGCTATATAAGAGATCAGCACAAAGACGAACGCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCGGCCCTTGTTCGTGTTAATCCCAGCTCAGAATa                                                                                                                                                                                                                                           <  6:216730/148‑1 (MQ=255)
                                                                                                                                                         agagaTAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGa                                                                                                                                                                                                                                              >  2:61967/1‑143 (MQ=255)
                                                                                                                                                             aTAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATaa                                                                                                                                                                                                                                          <  4:61966/143‑1 (MQ=255)
                                                                                                                                                               taagagacaggcgatataAGAGATCAGCACAATTAAGAAAGCTTTATTACACTGACAGGATGTGACATATTTAACAGATTAGCCACAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCGCAGAAGAACATTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCaaa                                                                                                     <  1:208032/261‑1 (MQ=255)
                                                                                                                                                                  gagaGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGAc                                                                                                                                                                                                                               >  7:21607/1‑149 (MQ=255)
                                                                                                                                                                     agcagcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACg                                                                                                                                                                                                                    >  4:503005/1‑157 (MQ=255)
                                                                                                                                                                      gcagcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCg                                                                                                                                                                                                                                                                >  6:545036/1‑112 (MQ=255)
                                                                                                                                                                       cagcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATAc                                                                                                                                                                                                                                                                             <  7:359674/98‑1 (MQ=255)
                                                                                                                                                                       cagcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATAc                                                                                                                                                                                                                                                                             >  8:359674/1‑98 (MQ=255)
                                                                                                                                                                               taAGAGTTAAGCACAATTACGAACCCTTGTTTACACGGACAGGATTTTACATATTTAAAACATTATCCCCAGTTAAAAGGGCAGGAATACCAGCCCTTGTTCTTGTTAATCCCAGCTCAGAATAACTTGCCCGTACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTaaa                                                                                    <  4:598926/275‑1 (MQ=255)
                                                                                                                                                                               taAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCa                                                                                                                                                                                                                  >  8:77502/1‑149 (MQ=255)
                                                                                                                                                                               taAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCa                                                                                                                                                                                                                  >  7:597256/1‑149 (MQ=255)
                                                                                                                                                                                  gagaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCg                                                                                                                                                                                                                                                                >  6:487394/1‑100 (MQ=255)
                                                                                                                                                                                   agaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGa                                                                                                                                                                                    <  2:411631/175‑1 (MQ=255)
                                                                                                                                                                                    gaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCg                                                                                                                                                                                                                                                                >  4:390647/1‑98 (MQ=255)
                                                                                                                                                                                     aaaacaaataTTACAAACTTATNTNTACACTGAAAGATATTTACATAGTTAACAGATTATCCCCAGATAANAGGGCTGGAATCACAGCCCTTGTTCGTGTTAATCCCAGCACCGAATAACTATCCCGGACTGGCCCGAACNCATAGCGTCAAAGCCCTTCTGTAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTaaa                                                                                    <  3:83644/260‑1 (MQ=255)
                                                                                                                                                                                      tCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATaa                                                                                                                                                                                                                                          <  4:316319/118‑1 (MQ=255)
                                                                                                                                                                                         gCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCaaa                                                                                                                                                                                                        <  8:438108/149‑1 (MQ=255)
                                                                                                                                                                                         gCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCaaa                                                                                                                                                                                                        >  8:655837/1‑149 (MQ=255)
                                                                                                                                                                                              aTTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGccc                                                                                                                                                                                                    >  6:369245/1‑148 (MQ=255)
                                                                                                                                                                                               ttACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCtt                                                                                                                                                                                                  <  8:301449/149‑1 (MQ=255)
                                                                                                                                                                                                      ctctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaaa                                                                                                                                                                                           >  8:163800/1‑149 (MQ=255)
                                                                                                                                                                                                       tctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCCCACGACCGTAAATACCTTTAATGCACAAGCCTTTa                                                             >  2:586892/1‑274 (MQ=255)
                                                                                                                                                                                                        ctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatc                                                                                                                                                                                         >  8:100814/1‑149 (MQ=255)
                                                                                                                                                                                                        ctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAat                                                                                                                                                                                          >  5:501530/1‑148 (MQ=255)
                                                                                                                                                                                                         tttATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGttt                                                                                                         >  1:65049/1‑228 (MQ=255)
                                                                                                                                                                                                           tattaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATa                                                                                                                                                                                                                <  3:203991/123‑1 (MQ=255)
                                                                                                                                                                                                           tattaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcat                                                                                                                                                                                       <  5:803306/148‑1 (MQ=255)
                                                                                                                                                                                                            attaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcatc                                                                                                                                                                                      <  6:168576/148‑1 (MQ=255)
                                                                                                                                                                                                             ttaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCg                                                                                                                                                                                     <  6:403678/148‑1 (MQ=255)
                                                                                                                                                                                                             ttaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGa                                                                                                                                                                                    >  7:72224/1‑149 (MQ=255)
                                                                                                                                                                                                             ttaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGa                                                                                                                                                                                    >  7:855716/1‑149 (MQ=255)
                                                                                                                                                                                                               acacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATa                                                                                                                                                                                  >  7:426965/1‑149 (MQ=255)
                                                                                                                                                                                                                 acTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATaga                                                                                                                                                                                >  8:47604/1‑149 (MQ=255)
                                                                                                                                                                                                                 acTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAgaga                                                                                                                                                                              >  2:782139/1‑151 (MQ=255)
                                                                                                                                                                                                                   tGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACg                                                                                                                                                                                                                    >  4:294703/1‑111 (MQ=255)
                                                                                                                                                                                                                   tGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCGTAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAAAGATGGGTAATGATCGGCGAAAGATCGAGGCAAGACTGAATCAGCGCCGCCATCTTGTACCAGGTATCGAACATCCCACGACCGTAAATACCTTTAATGAACAAACCTTTAAAGACCACCTTgg                                                >  2:324943/1‑275 (MQ=255)
                                                                                                                                                                                                                   tGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACGAGGTTTCAAACATCTCACGACCGTCAATACCTTTAATGAAAAAGCCTTTAAAGATCACTTTgg                                                >  2:193397/1‑275 (MQ=255)
                                                                                                                                                                                                                   tGACAGGAGTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGCATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAAAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCCGAATGAATCAGCTGCGCCAACTTGTATCAGGATTCAACCATCTTACGACCGTAAATACCTTTATAGAACAAACCATTAAAACACACTTTgg                                                >  4:667519/1‑275 (MQ=255)
                                                                                                                                                                                                                     aCAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTAAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACAGTAAATACCTTTAAAGCACAACCCTTTAAAGATAAATTTGGta                                              >  1:107211/1‑274 (MQ=255)
                                                                                                                                                                                                                      cAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAAc                                                                                                                                                                           <  7:77502/149‑1 (MQ=255)
                                                                                                                                                                                                                       aGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTc                                                                                                                                                                                                                                                                 <  3:48594/62‑1 (MQ=255)
                                                                                                                                                                                                                       aGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAg                                                                                                                                                                                                       <  2:369413/120‑1 (MQ=255)
                                                                                                                                                                                                                       aGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAAc                                                                                                                                                                           >  6:299826/1‑148 (MQ=255)
                                                                                                                                                                                                                         gATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGa                                                                                                                                                                         <  5:442628/148‑1 (MQ=255)
                                                                                                                                                                                                                            tttACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATggg                                                                                                                                                                     <  8:855716/149‑1 (MQ=255)
                                                                                                                                                                                                                            tttACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAAAGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGa                                       >  2:662704/1‑275 (MQ=255)
                                                                                                                                                                                                                             ttACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAATTCATCGATATAGACACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGACATCTATTACCAGGTTTCAAACATCTCCCGACCGTCAATACCTTTAATGAACAAGCCTTTAAAGATCAATTTGGTCCAATAGAt                                      >  4:432348/1‑275 (MQ=255)
                                                                                                                                                                                                                              tACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCaa                                                                                                                                                                                                         <  1:187159/111‑1 (MQ=255)
                                                                                                                                                                                                                               aCATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTaa                                                                                                                                                                  <  7:315841/149‑1 (MQ=255)
                                                                                                                                                                                                                                     ttAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATcc                                                                                                                                                                                                                                                       >  4:793902/1‑58 (MQ=255)
                                                                                                                                                                                                                                      tAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCg                                                                                                                                                            >  6:654913/1‑148 (MQ=255)
                                                                                                                                                                                                                                         cAGATTATCCTCAGCAAAAAGGGCTGGAATACCAGCCCTTTTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTTGCCCGAACGCATAGCGTCAAAGCCCTTCTGGACATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGAt                                      <  4:161013/263‑1 (MQ=255)
                                                                                                                                                                                                                                           gATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGaa                                                                                                                                                       <  5:891888/148‑1 (MQ=255)
                                                                                                                                                                                                                                           gATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGaa                                                                                                                                                       <  6:100341/148‑1 (MQ=255)
                                                                                                                                                                                                                                              tATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCg                                                                                                                                                                                                                                                                >  4:347936/1‑40 (MQ=255)
                                                                                                                                                                                                                                              tATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGc                                                                  >  4:220146/1‑230 (MQ=255)
                                                                                                                                                                                                                                               aTCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATc                                                                                                                                                  <  8:21607/149‑1 (MQ=255)
                                                                                                                                                                                                                                               aTCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGc                                                                                                                                             >  2:502764/1‑154 (MQ=255)
                                                                                                                                                                                                                                                 ccTCAGATAAAGNGGCGGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGAGCGAGGCCAGACTGAATCAGcgccg                                                                                                                           <  4:758101/170‑1 (MQ=255)
                                                                                                                                                                                                                                                 ccTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCg                                                                                                                                                 <  6:99655/148‑1 (MQ=255)
                                                                                                                                                                                                                                                  cTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGa                                                                                                                                                                                    <  2:259030/112‑1 (MQ=255)
                                                                                                                                                                                                                                                  cTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGa                                                                                                                                                <  5:67651/148‑1 (MQ=255)
                                                                                                                                                                                                                                                   tCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGa                                                                                                                                                                         <  7:670673/122‑1 (MQ=255)
                                                                                                                                                                                                                                                   tCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAgg                                                                                                                                              >  8:529027/1‑149 (MQ=255)
                                                                                                                                                                                                                                                   tCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAg                                                                                                                                               >  5:597859/1‑148 (MQ=255)
                                                                                                                                                                                                                                                    cAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCgg                                                                                                                                                           <  5:164329/135‑1 (MQ=255)
                                                                                                                                                                                                                                                    cAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGc                                                                                                                                             <  7:39526/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     aGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGcc                                                                                                                                            >  7:284702/1‑149 (MQ=255)
                                                                                                                                                                                                                                                      gATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGcc                                                                                                                                            >  6:601716/1‑148 (MQ=255)
                                                                                                                                                                                                                                                       aTAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGCCTGAATCAGCGCCGCCATCTTGTACCCGGTTTCAAACATCTCACGACCTTACATACCTTTAATGAACCAGCCTTTACAGATCACTTTGGTCCAGTCGATCGCCATATCAGCCTGCGGCATAccc             >  1:153771/1‑274 (MQ=255)
                                                                                                                                                                                                                                                          aaaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTaaa                                                                                    <  2:302072/200‑1 (MQ=255)
                                                                                                                                                                                                                                                           aaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACg                                                                                                                                                                          <  1:48873/113‑1 (MQ=255)
                                                                                                                                                                                                                                                            aGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTg                                                                                                                                      >  6:884560/1‑148 (MQ=255)
                                                                                                                                                                                                                                                             gggCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGaa                                                                                                                                    <  7:496185/149‑1 (MQ=255)
                                                                                                                                                                                                                                                              ggCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGaa                                                                                                                                    >  5:367275/1‑148 (MQ=255)
                                                                                                                                                                                                                                                               gCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATc                                                                                                                                  >  7:454273/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               gCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATc                                                                                                                                  >  7:80233/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                cTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATc                                                                                                                                  <  5:830206/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                   gAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAgcg                                                                                                                              >  8:261927/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    aaTACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCa                                                                                                                                 >  5:535790/1‑145 (MQ=255)
                                                                                                                                                                                                                                                                     aTACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACaa                                                                    <  4:481963/205‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aTACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGAt                                                        <  2:447597/217‑1 (MQ=255)
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                            
TAACGGAATAATAACACTAAGTTTATTGGTGCTGTTCATCATAAATTACCGACGACCTCTGAAGTTATCCTGACGTTTTACAGGCTTTGTATTAAGAAACTATTAAGCCTGTGAGATATTATTTTTACAGGTGAGAGAGTCACTTCAGGGTAGAGAGATAATGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAACGATGGGTAATGATCGGCGAAAGATCGAGGCCAGACTGAATCAGCGCCGCCATCTTGTACCAGGTTTCAAACATCTCACGACCGTAAATACCTTTAATGAACAAGCCTTTAAAGATCACTTTGGTCCAGTCGATAGACATATCAGACGGCGGAATACCCAGCATCGCAAT  >  NC_000913/3790040‑3790571

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