Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,790,303 A→T 16.1% 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,3030AT16.1% 509.9 / 73.2 198intergenic (‑222/+17)waaH/tdhLPS(HepIII)‑glucuronic acid glycosyltransferase/L‑threonine 3‑dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (86/79);  new base T (13/19);  total (99/99)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

Only 100 of 201 total aligned reads displayed.

GAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAAC  >  NC_000913/3790204‑3790402
                                                                                                   |                                                                                                   
  gcagcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATAcc                                                                                                   <  1:688541/100‑1 (MQ=255)
    agcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATAcc                                                                                                   <  2:2705986/98‑1 (MQ=255)
    agcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAg                                                                                                 >  1:3833736/1‑100 (MQ=255)
     gcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGc                                                                                                <  2:1580214/100‑1 (MQ=255)
       tatataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGccc                                                                                              >  2:3332768/1‑100 (MQ=255)
       tatataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGCCTGGAATACCAGccc                                                                                              >  1:3867027/1‑100 (MQ=255)
        atataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCt                                                                                             >  1:3888564/1‑100 (MQ=255)
         tataTGAGATTAGTATAATTAATATTTTTTTATTACATTGACAGTATTTTACATATTTAACAGATTATCTTTAGATAAAAGGTCTGGTATTCCAGcctt                                                                                             <  2:2371199/99‑3 (MQ=255)
         tataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCtt                                                                                            <  2:1103234/100‑1 (MQ=255)
          ataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATTCTCAGATAAAAGGGCTGGTATCCCAGCCCttt                                                                                           <  2:1912984/100‑2 (MQ=255)
          ataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTg                                                                                           <  2:3146715/100‑1 (MQ=255)
           taAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCtttt                                                                                          <  2:2847221/100‑3 (MQ=255)
           taAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGt                                                                                          >  2:351994/1‑100 (MQ=255)
           taAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTAAGATAACAAGGCTGGAATATCAGCACTTGt                                                                                          >  2:1167873/1‑100 (MQ=255)
            aaGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGAATATCCTCAGATAAAAGGGCTGGAATAACAGCCCTTGtt                                                                                         >  1:1630243/1‑100 (MQ=255)
             agagATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttat                                                                                        <  1:620031/100‑5 (MQ=255)
             agagATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTc                                                                                        >  1:1132802/1‑100 (MQ=255)
             agagATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTc                                                                                        >  1:242778/1‑100 (MQ=255)
              gagaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatc                                                                                       <  2:153252/100‑6 (MQ=39)
              gagaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCg                                                                                       <  1:3300177/100‑1 (MQ=255)
                naTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtg                                                                                     >  2:676614/2‑100 (MQ=255)
                gaTCAGCACAATTACTAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtg                                                                                     <  1:3876312/100‑1 (MQ=255)
                 aTCAGCATAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtgt                                                                                    <  1:3738852/100‑1 (MQ=255)
                 aTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTATATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtgt                                                                                    <  2:2501990/100‑1 (MQ=255)
                 aTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatctga                                                                                    <  1:552074/100‑9 (MQ=38)
                 aTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatctga                                                                                    <  1:1986416/100‑9 (MQ=38)
                 aTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtgt                                                                                    <  1:231781/100‑1 (MQ=255)
                 aTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtg                                                                                     <  2:3804802/99‑1 (MQ=255)
                  tCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGtt                                                                                   <  2:2789147/100‑1 (MQ=255)
                   cAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTa                                                                                  >  2:2141957/1‑100 (MQ=255)
                    aGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTaa                                                                                 <  2:3914823/100‑1 (MQ=255)
                    aGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTaa                                                                                 >  1:3768729/1‑100 (MQ=255)
                     gCACAATTACGAACTCTTTATTACAATGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCCGGAATACCAACCCTTGTTCGTGTTaa                                                                                 >  2:3235310/1‑99 (MQ=255)
                      cacaATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATc                                                                               <  1:2928492/100‑1 (MQ=255)
                       acaATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATcc                                                                              <  1:2956799/100‑1 (MQ=255)
                        caATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTAAAATCAAATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATccc                                                                             >  2:1318422/1‑100 (MQ=255)
                          aTTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAg                                                                           >  2:630045/1‑100 (MQ=255)
                            tACGAACTATCTAATACACTGAAAAAATTCTACATATTTAACTTAATATCCTCAGATAAAAGGGCTGCAATACCAGCCCTTGTTCGTGTTAATCACAGCt                                                                         <  2:2747545/100‑1 (MQ=255)
                                aaCTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGa                                                                     <  1:3069494/100‑1 (MQ=255)
                                  ttccgatctaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAAt                                                                   <  1:2478330/92‑1 (MQ=255)
                                  ctctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAAt                                                                   >  1:526454/1‑100 (MQ=255)
                                  ctctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAAt                                                                   <  2:2522868/100‑1 (MQ=255)
                                  ctctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAAt                                                                   >  1:2204128/1‑100 (MQ=255)
                                   tctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATa                                                                  >  1:2792626/1‑100 (MQ=255)
                                   tctGTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAAt                                                                   <  1:2275109/99‑1 (MQ=255)
                                    ctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATaa                                                                 >  1:2451721/1‑100 (MQ=255)
                                    cgttattaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATGAAAGAGCTGGAATACTAGCCCTTGTTCGTTTTAATCCCAGCTgagact                                                                   >  2:3725508/3‑96 (MQ=255)
                                       tattaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACtt                                                              >  2:28817/1‑100 (MQ=255)
                                         ttaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACtttn                                                            <  1:676622/100‑2 (MQ=255)
                                           acacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTccc                                                          >  1:1516188/1‑100 (MQ=255)
                                           acacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTccc                                                          >  1:2565821/1‑100 (MQ=255)
                                           acacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTAGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTccc                                                          >  1:7407/1‑100 (MQ=255)
                                            cacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCg                                                         <  2:228430/100‑1 (MQ=255)
                                            cacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCg                                                         <  1:1363048/100‑1 (MQ=255)
                                            cacTGAAAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGTCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCg                                                         <  2:1260479/100‑1 (MQ=255)
                                                   aGGATTTTACATATTTAACAGATTATCCTCAGATAAAGGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGc                                                  >  1:3083010/1‑100 (MQ=255)
                                                   aGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGc                                                  <  2:1354710/100‑1 (MQ=255)
                                                     gATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGccc                                                <  1:3306101/100‑1 (MQ=255)
                                                      aTTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCg                                               <  2:3078061/100‑1 (MQ=255)
                                                            cATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACGAGCCCTTGTTCGTGTTAATGCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCa                                         >  2:2127892/1‑100 (MQ=255)
                                                            cATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCa                                         >  1:3839277/1‑100 (MQ=255)
                                                            cATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATACCAGAACAGAATAACTTTCCCGGACTGGCCCGAACGCa                                         >  1:1404418/1‑100 (MQ=255)
                                                            cATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCATGATAATCCAAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCa                                         >  1:1016434/1‑100 (MQ=255)
                                                              tatTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCa                                         <  2:250917/98‑1 (MQ=255)
                                                              tatTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATa                                       <  2:2842059/100‑1 (MQ=255)
                                                               atTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAg                                      <  2:867330/100‑1 (MQ=255)
                                                               atTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAg                                      <  1:1942576/100‑1 (MQ=255)
                                                                tttAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGc                                     <  1:1842012/100‑1 (MQ=255)
                                                                    aCAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCa                                 <  1:3236027/100‑1 (MQ=255)
                                                                       gATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAg                              <  2:3357824/100‑1 (MQ=255)
                                                                       gATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAg                              >  2:2874711/1‑100 (MQ=255)
                                                                         ttATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTATTCCCAGCTCAGACTAACTTTCCCGGACTGGCCCGAACGCATAGCTTCAAAGc                             >  1:1265710/1‑99 (MQ=255)
                                                                              cTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCtt                         >  1:3627870/1‑98 (MQ=255)
                                                                              cTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCt                       >  2:2299163/1‑100 (MQ=255)
                                                                               tCAGATAAAAGGGCTGGAATACGAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCACGGACTGGCCCGAACGCATAGAGTCAAAGCCCGTCTg                      >  2:447340/1‑100 (MQ=255)
                                                                                 acgaacaaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTAGCCCGAACACATAGCGTCAAAGCCCTTCTGGa                    >  1:1232958/7‑100 (MQ=38)
                                                                                  cgaacaaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaa                   >  2:2453342/6‑100 (MQ=39)
                                                                                  cgaacaaGGGCTGGTATTCCAGCCCATGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGa                    >  2:1416271/6‑99 (MQ=38)
                                                                                  cgaacaaGGGCTGGTATTCCAGCCATTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTATGGaa                   >  2:138057/6‑100 (MQ=255)
                                                                                   aTAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaaa                  <  1:1642906/100‑1 (MQ=255)
                                                                                   aTAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaaa                  <  2:3589203/100‑1 (MQ=255)
                                                                                   aTAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGa                    <  1:3554020/98‑1 (MQ=255)
                                                                                    tAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAat                 <  1:2775048/100‑1 (MQ=255)
                                                                                     aaaaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatc                >  2:1759612/1‑100 (MQ=255)
                                                                                     aaaaGGGCTGGAATACCAGACCTTGTTCGTGTTAATCAAAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGGCATTCTGGAAat                 >  2:786141/1‑99 (MQ=255)
                                                                                      caaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTTGCCCGAACCCCTAGCGTCAAAGCCCTTCTGGCAatca               >  2:3016929/2‑100 (MQ=255)
                                                                                      aaaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatca               >  1:3686149/1‑100 (MQ=255)
                                                                                      aaaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatca               <  2:1850053/100‑1 (MQ=255)
                                                                                       aaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAACGCCCTTCCGGaa                   >  2:3443935/1‑95 (MQ=255)
                                                                                       aaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGGCAAAGCCCTTCTGGAAatcat              >  2:3384575/1‑100 (MQ=255)
                                                                                        aGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcatc             >  1:836088/1‑100 (MQ=255)
                                                                                        aGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcatc             <  1:267552/100‑1 (MQ=255)
                                                                                        aGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcatc             <  1:2127911/100‑1 (MQ=255)
                                                                                        aGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcatc             <  1:3521671/100‑1 (MQ=255)
                                                                                         gggCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCg            <  2:576294/100‑1 (MQ=255)
                                                                                              ggAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATaga       >  1:1465434/1‑100 (MQ=255)
                                                                                              ggAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACACATAGCGTCAAAGCCCTTCTGGAAATCATCGATaga       >  1:2253997/1‑100 (MQ=255)
                                                                                               gAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATagag      <  1:878963/100‑1 (MQ=255)
                                                                                                 aTACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGaa    >  1:2354142/1‑100 (MQ=255)
                                                                                                  tACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGaaa   <  1:988698/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAGAAAC  >  NC_000913/3790204‑3790402

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