Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,790,300 A→T 16.9% intergenic (‑219/+20) 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,3000AT16.9% 485.3 / 81.2 195intergenic (‑219/+20)waaH/tdhLPS(HepIII)‑glucuronic acid glycosyltransferase/L‑threonine 3‑dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (85/77);  new base T (13/20);  total (98/97)
Fisher's exact test for biased strand distribution p-value = 1.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.44e-01

Only 100 of 197 total aligned reads displayed.

TGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAG  >  NC_000913/3790201‑3790398
                                                                                                   |                                                                                                  
gatctgcagcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGa                                                                                                    <  1:626925/95‑1 (MQ=255)
   gagCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATa                                                                                                 <  1:3591820/100‑1 (MQ=255)
     gcagcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGaagatc                                                                                               >  2:626918/1‑96 (MQ=255)
     gcagcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATAcc                                                                                               >  1:2348693/1‑100 (MQ=255)
     gcagcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATAcc                                                                                               <  1:688541/100‑1 (MQ=255)
       agcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATAcc                                                                                               <  2:2705986/98‑1 (MQ=255)
       agcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAg                                                                                             >  1:156075/1‑100 (MQ=255)
       agcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAg                                                                                             >  1:1850068/1‑100 (MQ=255)
        gcTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGc                                                                                            >  1:3642516/1‑100 (MQ=255)
         cTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGcc                                                                                           >  2:1363033/1‑100 (MQ=255)
          tatataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGccc                                                                                          >  2:3332768/1‑100 (MQ=255)
           atataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCt                                                                                         <  1:3775702/100‑1 (MQ=255)
            tataTGAGATTAGTATAATTAATATTTTTTTATTACATTGACAGTATTTTACATATTTAACAGATTATCTTTAGATAAAAGGTCTGGTATTCCAGcctt                                                                                         <  2:2371199/99‑3 (MQ=255)
            tataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCtt                                                                                        <  2:1103234/100‑1 (MQ=255)
            tataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCt                                                                                         >  1:3763768/1‑99 (MQ=255)
             ataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTg                                                                                       <  2:3146715/100‑1 (MQ=255)
             ataAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATGTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTg                                                                                       >  2:2867359/1‑100 (MQ=255)
              taAGAGATCAGCACAATTACTAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCtttt                                                                                      <  1:781984/100‑3 (MQ=255)
              taAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGt                                                                                      >  2:351994/1‑100 (MQ=255)
              taAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTAAGATAACAAGGCTGGAATATCAGCACTTGt                                                                                      >  2:1167873/1‑100 (MQ=255)
                agagATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttat                                                                                    <  2:1562203/100‑5 (MQ=39)
                agagATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTc                                                                                    >  1:242778/1‑100 (MQ=255)
                 gagaTCAGCACAATTACGATCTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAATGGCTGTAATACCAGCGCTTGTTCg                                                                                   <  2:1407745/100‑1 (MQ=255)
                 gagaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatc                                                                                   <  2:949264/100‑6 (MQ=39)
                 gagaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCg                                                                                   <  1:3300177/100‑1 (MQ=255)
                  agaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatc                                                                                   <  1:3862389/99‑6 (MQ=39)
                  agaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgt                                                                                  >  2:762402/1‑100 (MQ=255)
                   gaTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtg                                                                                 <  1:1647895/100‑1 (MQ=255)
                    aTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTATATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtgt                                                                                <  2:2501990/100‑1 (MQ=255)
                    aTCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtg                                                                                 <  2:3804802/99‑1 (MQ=255)
                    aTCAGCACAATTACCAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCgtg                                                                                 <  1:3876395/99‑1 (MQ=255)
                     tCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGTATTCCAGCCCttttatctgag                                                                               <  1:1484760/100‑10 (MQ=38)
                     tCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGtt                                                                               <  2:2789147/100‑1 (MQ=255)
                       aGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTaa                                                                             <  2:3914823/100‑1 (MQ=255)
                        gCACAATTACGAACTCTTTATTACAATGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCCGGAATACCAACCCTTGTTCGTGTTaa                                                                             >  2:3235310/1‑99 (MQ=255)
                         cacaATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATc                                                                           <  1:2928492/100‑1 (MQ=255)
                          acaATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATcc                                                                          <  1:2956799/100‑1 (MQ=255)
                          acaATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATc                                                                           >  2:806201/1‑99 (MQ=255)
                            aaTTACGACCTCTTGATTACACTGACAGGATGTTACATGTTTAACAGATTATTCTCAGATTAAAGGGCTGGTATTCCAGCCCttttatctgag                                                                               <  2:1766543/93‑10 (MQ=21)
                               tACGAACTATCTAATACACTGAAAAAATTCTACATATTTAACTTAATATCCTCAGATAAAAGGGCTGCAATACCAGCCCTTGTTCGTGTTAATCACAGCt                                                                     <  2:2747545/100‑1 (MQ=255)
                                aCGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTc                                                                    >  1:2217928/1‑100 (MQ=255)
                                 cGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTc                                                                    <  2:1443281/99‑1 (MQ=255)
                                   aaCTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGa                                                                 <  1:3069494/100‑1 (MQ=255)
                                    aCTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGaa                                                                >  1:1945927/1‑100 (MQ=255)
                                     ttccgatctaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAAt                                                               <  1:2478330/92‑1 (MQ=255)
                                     ctctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAAt                                                               <  2:2522868/100‑1 (MQ=255)
                                      tctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATa                                                              >  1:2792626/1‑100 (MQ=255)
                                       ctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATaa                                                             >  1:2451721/1‑100 (MQ=255)
                                       ctTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTAGTGTTAATCCCAGCTCAGAATaa                                                             >  1:1642171/1‑100 (MQ=255)
                                       cgttattaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATGAAAGAGCTGGAATACTAGCCCTTGTTCGTTTTAATCCCAGCTgagact                                                               >  2:3725508/3‑96 (MQ=255)
                                            ttaCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACtttn                                                        <  1:676622/100‑2 (MQ=255)
                                             taCACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCATCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTcc                                                       <  2:1691642/100‑1 (MQ=255)
                                              acacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTccc                                                      >  1:2565821/1‑100 (MQ=255)
                                              acacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCACTTGTTCGTATTAATACCAGCTCAGAATAACTTTcc                                                       >  1:954083/1‑99 (MQ=255)
                                               cacTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCg                                                     <  2:228430/100‑1 (MQ=255)
                                               cacTGAAAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGTCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCg                                                     <  2:1260479/100‑1 (MQ=255)
                                                 cTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCATCCCTTGTTCGTGTTAATCCTAACTCAGAATaaaataaaaaaa                                                   >  1:1457069/1‑90 (MQ=255)
                                                   gACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACt                                                 >  2:2030698/1‑100 (MQ=255)
                                                     cAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTgg                                               <  1:541886/100‑1 (MQ=255)
                                                      aGGATTTTACATATTTAACAGATTATCCTCAGATAAAGGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGc                                              >  1:3083010/1‑100 (MQ=255)
                                                       ggATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGcc                                             <  1:2113766/100‑1 (MQ=255)
                                                        gATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGccc                                            <  1:3306101/100‑1 (MQ=255)
                                                         aTTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCg                                           <  2:3078061/100‑1 (MQ=255)
                                                            ttACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAAc                                        <  1:515033/100‑1 (MQ=255)
                                                              aCATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGc                                      >  1:1390206/1‑100 (MQ=255)
                                                               cATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGACTAACTTTCCCGCACTGGCCCGAACGCa                                     >  2:214112/1‑100 (MQ=255)
                                                               cATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCATGATAATCCAAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCa                                     >  1:1016434/1‑100 (MQ=255)
                                                                 tatTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCa                                     <  2:250917/98‑1 (MQ=255)
                                                                 tatTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATa                                   >  1:226360/1‑100 (MQ=255)
                                                                 tatTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATa                                   <  2:2842059/100‑1 (MQ=255)
                                                                   tttAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATa                                   >  1:860878/1‑98 (MQ=255)
                                                                    ttAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCg                                <  1:2185943/100‑1 (MQ=255)
                                                                    ttAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCg                                <  2:2876828/100‑1 (MQ=255)
                                                                       aCAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCa                             <  1:3236027/100‑1 (MQ=255)
                                                                          gATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAg                          <  2:3357824/100‑1 (MQ=255)
                                                                           aTTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGGTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGc                         <  2:13357/100‑1 (MQ=255)
                                                                              aTCCTCAGATAAAAGGGCTAGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGccc                       >  2:635165/1‑99 (MQ=255)
                                                                                 cTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCt                   >  2:2299163/1‑100 (MQ=255)
                                                                                 cTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCt                   >  1:576301/1‑100 (MQ=255)
                                                                                 cTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCt                   <  1:2299208/100‑1 (MQ=255)
                                                                                  tCAGATAAAAGGGCTGGAATACGAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCACGGACTGGCCCGAACGCATAGAGTCAAAGCCCGTCTg                  >  2:447340/1‑100 (MQ=255)
                                                                                   cAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTgg                 >  1:3542574/1‑100 (MQ=255)
                                                                                     gATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaa               >  1:1312237/1‑100 (MQ=255)
                                                                                     cgaacaaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaa               >  2:2453342/6‑100 (MQ=39)
                                                                                     cgaacaaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaa               >  1:3737291/6‑100 (MQ=39)
                                                                                      aTAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaaa              <  2:3589203/100‑1 (MQ=255)
                                                                                      aTAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGaaa              <  1:1230008/100‑1 (MQ=255)
                                                                                       tAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAat             >  1:3835146/1‑100 (MQ=255)
                                                                                       tAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAat             <  1:2775048/100‑1 (MQ=255)
                                                                                       tAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAat             <  2:2181754/100‑1 (MQ=255)
                                                                                        aaaaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCAGGACTAGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatc            >  1:1065678/1‑100 (MQ=255)
                                                                                         caaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTTGCCCGAACCCCTAGCGTCAAAGCCCTTCTGGCAatca           >  2:3016929/2‑100 (MQ=255)
                                                                                          aaGGGCTGGTATTCCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCCTCTGGAAatcat          >  2:1908229/1‑100 (MQ=255)
                                                                                          aaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAACGCCCTTCCGGaa               >  2:3443935/1‑95 (MQ=255)
                                                                                          aaGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGGCAAAGCCCTTCTGGAAatcat          >  2:3384575/1‑100 (MQ=255)
                                                                                           aGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcatc         <  1:3521671/100‑1 (MQ=255)
                                                                                           aGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAatcatc         <  1:267552/100‑1 (MQ=255)
                                                                                            gggCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCg        <  2:576294/100‑1 (MQ=255)
                                                                                               cTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATa     >  1:3871127/1‑100 (MQ=255)
                                                                                                 ggAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATaga   <  1:1258765/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
TGAGAGCAGCTATATAAGAGATCAGCACAATTACGAACTCTTTATTACACTGACAGGATTTTACATATTTAACAGATTATCCTCAGATAAAAGGGCTGGAATACCAGCCCTTGTTCGTGTTAATCCCAGCTCAGAATAACTTTCCCGGACTGGCCCGAACGCATAGCGTCAAAGCCCTTCTGGAAATCATCGATAGAG  >  NC_000913/3790201‑3790398

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: