Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,580,094 C→T 100% D233N (GAT→AAT)  ydeM ← putative YdeN‑specific sulfatase‑maturating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,580,0940CT98.4% 404.2 / ‑2.8 122D233N (GAT→AAT) ydeMputative YdeN‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (60/60);  total (60/62)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

Only 100 of 122 total aligned reads displayed.

CCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATG  >  NC_000913/1579998‑1580190
                                                                                                |                                                                                                
ccTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGtt                                                                                               >  2:1367262/1‑100 (MQ=255)
 cTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGttt                                                                                              >  2:129769/1‑100 (MQ=255)
  ttGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGtttt                                                                                             <  2:2175009/100‑1 (MQ=255)
   tGCAGGACTCCTGGAAAATACATCTGGTATGCCCATTCCCAAAAAAACAGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGttttt                                                                                            >  1:1916220/1‑100 (MQ=255)
    gCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTa                                                                                           <  2:2528758/100‑1 (MQ=255)
    gCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTa                                                                                           >  2:3445400/1‑100 (MQ=255)
    gCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTa                                                                                           >  2:1996269/1‑100 (MQ=255)
    gCAGGACTCCTGGAAAATACAACTGGTATACCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGttttt                                                                                            <  1:672241/99‑1 (MQ=255)
     cAGGACTCCTGGCAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTa                                                                                           >  1:2547735/1‑99 (MQ=255)
     cAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTa                                                                                           <  1:481366/99‑1 (MQ=255)
      aGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAAc                                                                                         <  2:837656/100‑1 (MQ=255)
      aGGACTCCTGGAAAATACAACTGGTATACAAAAAACAAAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAAc                                                                                         <  1:1519730/100‑1 (MQ=255)
           tCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCAtt                                                                                    <  2:3531453/100‑1 (MQ=255)
            ccTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTg                                                                                   <  2:1613212/100‑1 (MQ=255)
             cTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGAAAAAGCTTTCAAACTGACGGATGAAAATGTCACCCACATTGTTTTTAACACATTGc                                                                                  >  1:3286260/1‑100 (MQ=255)
              tGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCa                                                                                 <  2:1986238/100‑1 (MQ=255)
              tGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCa                                                                                 >  1:914451/1‑100 (MQ=255)
              tGGAAAATACAACTGGTATGAAAAAAAAAAAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCa                                                                                 <  1:880254/100‑1 (MQ=255)
               ggAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCAt                                                                                >  2:129083/1‑100 (MQ=255)
               ggAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCAt                                                                                >  2:1643513/1‑100 (MQ=255)
               ggAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCAt                                                                                >  2:1989379/1‑100 (MQ=255)
               ggAAAATACAACTGGTATGCCCATTACCCAAAAGACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGc                                                                                  >  1:3918886/1‑98 (MQ=255)
                gAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATa                                                                               >  2:1869469/1‑100 (MQ=255)
                  aaaTACAACTGGTATGCCCATTCAAAAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATaaa                                                                             <  1:4092265/100‑1 (MQ=255)
                    aTACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATaaaaa                                                                           >  2:904796/1‑100 (MQ=255)
                       cAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATgg                                                                        >  1:3011430/1‑100 (MQ=255)
                        aaCTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGt                                                                       <  2:934166/100‑1 (MQ=255)
                        aaCTGATATTCCCATTCCCCAAAAAACGTCTTACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTTTTTTTAACCCATTGCATAAAAATGGt                                                                       <  2:1848290/100‑1 (MQ=255)
                         aCTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGtt                                                                      <  1:3753340/100‑1 (MQ=255)
                         aCTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGtt                                                                      >  2:2074979/1‑100 (MQ=255)
                          cTGGTATGCCCATTCCCCAAAAATCGGCTGACAAAGCTTTCAAACTGTCGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTg                                                                     <  2:597582/100‑1 (MQ=255)
                          cTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTg                                                                     <  2:1167625/100‑1 (MQ=255)
                          cTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTg                                                                     <  1:3066506/100‑1 (MQ=255)
                          cTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTg                                                                     >  2:3461713/1‑100 (MQ=255)
                             gTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACa                                                                  <  2:2342239/100‑1 (MQ=255)
                              tATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACAt                                                                 >  1:2084164/1‑100 (MQ=255)
                               aTGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACAtg                                                                <  2:437640/100‑2 (MQ=255)
                                tGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATaa                                                               <  2:3032535/100‑1 (MQ=255)
                                 gCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATaaa                                                              >  2:2187648/1‑100 (MQ=255)
                                  cccATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAAc                                                             >  2:1396535/1‑100 (MQ=255)
                                     aTTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTg                                                          <  1:1375552/100‑1 (MQ=255)
                                      ttCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGc                                                         <  1:1496104/100‑1 (MQ=255)
                                       tCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTggc                                                        <  2:3961402/100‑3 (MQ=255)
                                        ccccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCa                                                       >  2:1342489/1‑100 (MQ=255)
                                         cccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCAt                                                      >  1:1084997/1‑100 (MQ=255)
                                          ccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATa                                                     >  1:253203/1‑100 (MQ=255)
                                             aaaaaCGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGc                                                  <  2:3791031/100‑1 (MQ=255)
                                              aaaaCGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGcc                                                 >  2:376150/1‑100 (MQ=255)
                                              aaaaCGGCTGACAAAACTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTCGACATAAACTTGCCATAAGcc                                                 <  2:3651560/100‑1 (MQ=255)
                                               aaaCGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCg                                                >  1:2465698/1‑100 (MQ=255)
                                                  cGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTg                                              >  1:2941225/1‑99 (MQ=255)
                                                   ggCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTggg                                            >  2:3293017/1‑100 (MQ=255)
                                                   ggCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTggg                                            >  1:1694498/1‑100 (MQ=255)
                                                    gCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGa                                           >  1:3962486/1‑100 (MQ=255)
                                                    gCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGa                                           <  2:3963128/100‑1 (MQ=255)
                                                     cTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAg                                          >  2:2167710/1‑100 (MQ=255)
                                                       gACAAAGCTTTCAAACTGACGGATGAAAATTTGACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGAGATAAACTTGCCATAAGCCGTGGGAGGc                                        >  2:4058287/1‑100 (MQ=255)
                                                       gACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGc                                        >  1:1241205/1‑100 (MQ=255)
                                                         cAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGcac                                      <  2:2756559/100‑1 (MQ=255)
                                                          aaaGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGcaca                                     <  2:3210953/100‑1 (MQ=255)
                                                           aaGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAg                                    >  1:1915450/1‑100 (MQ=255)
                                                           aaGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAg                                    >  2:3831818/1‑100 (MQ=255)
                                                            aGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGa                                   >  2:1865563/1‑100 (MQ=255)
                                                              cTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGaaa                                 >  2:128390/1‑100 (MQ=255)
                                                               tttCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGaaaa                                >  2:3121048/1‑100 (MQ=255)
                                                                 tCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAAt                              <  1:2462174/100‑1 (MQ=255)
                                                                  cAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATc                             <  2:966641/100‑1 (MQ=255)
                                                                  cAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATc                             <  2:4175750/100‑1 (MQ=255)
                                                                   aaaCTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCa                            <  1:18168/100‑1 (MQ=255)
                                                                   aaaCTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCa                            <  1:2818576/100‑1 (MQ=255)
                                                                   aaaCTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCa                            <  1:1989375/100‑1 (MQ=255)
                                                                   aaaCTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCa                            <  2:1451078/100‑1 (MQ=255)
                                                                   aaaCTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCa                            <  2:571047/100‑1 (MQ=255)
                                                                   aaaCTGACGGATGAAAATTTCACCCACATTGCTTTTAAACCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCa                            >  1:1005186/1‑100 (MQ=255)
                                                                     aCTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAAt                          <  1:622790/100‑1 (MQ=255)
                                                                      cTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATg                         <  2:2941237/100‑1 (MQ=255)
                                                                       tGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGa                        >  1:165960/1‑100 (MQ=255)
                                                                        gACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGa                        >  1:3692959/1‑99 (MQ=255)
                                                                        gACAGATGAAAATTTCAGCCACATTGTTTTAAACCCATGGTATAAGATGGGTTGACATGACCTTGCCGTAAGGCAGTGAAGGCACAGAAAAATCAAt                          <  1:395604/97‑1 (MQ=255)
                                                                         aCGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATg                         <  1:1528269/97‑1 (MQ=255)
                                                                            gATGAAAATTTCACCCACATTGTTTTGAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACgg                   <  1:2991416/100‑1 (MQ=255)
                                                                             aTGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGa                  >  2:3654535/1‑100 (MQ=255)
                                                                              tGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGaa                 >  1:2865233/1‑100 (MQ=255)
                                                                              tGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGaa                 >  1:3299221/1‑100 (MQ=255)
                                                                               gAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAAt                >  2:2883900/1‑100 (MQ=255)
                                                                               gAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAAt                <  1:3111239/100‑1 (MQ=255)
                                                                               gAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAAt                >  2:2796904/1‑100 (MQ=255)
                                                                                 aaaTTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAAtgt              >  2:3216512/1‑100 (MQ=255)
                                                                                 aaaTTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAAtgt              >  2:3236127/1‑100 (MQ=255)
                                                                                  aaTTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAAtgtg             >  2:583302/1‑100 (MQ=255)
                                                                                    tttCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGtt           >  2:1386429/1‑100 (MQ=255)
                                                                                     ttCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTtc          >  1:446300/1‑100 (MQ=255)
                                                                                       cACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTtctc        <  1:129083/100‑1 (MQ=255)
                                                                                        aCCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCa       >  1:995383/1‑100 (MQ=255)
                                                                                         cccACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGGGGGAGGCACAGAAAAATCAATGTTACGGAATGTGTtatc        >  1:3873607/1‑98 (MQ=255)
                                                                                          ccACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACt     <  2:3930693/100‑1 (MQ=255)
                                                                                          ccACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGAAGAAACAAAAAAATCAATGATACGGAATGTGTTCTCACt     >  1:2948534/1‑100 (MQ=255)
                                                                                            acaTTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTAt   <  1:750231/100‑1 (MQ=255)
                                                                                            acaTTGTTTTTAACCCATTGCATAAAAAAGGTTGACATGAACTTGCCATAAGCCGTGAGAGAAAAAAAAAAATCAATGATACGGAATGTGTTCTCACTa    <  1:756921/99‑1 (MQ=255)
                                                                                             caTTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATg  <  1:3029157/100‑1 (MQ=255)
                                                                                                |                                                                                                
CCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATG  >  NC_000913/1579998‑1580190

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