Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,288,032 C→A 16.2% intergenic (+22/‑58) yraH → / → yraI putative fimbrial‑like adhesin protein/putative periplasmic pilin chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,288,0320CA16.2% 487.0 / 71.5 179intergenic (+22/‑58)yraH/yraIputative fimbrial‑like adhesin protein/putative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base C (81/68);  new base A (15/14);  total (96/83)
Fisher's exact test for biased strand distribution p-value = 8.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 191 total aligned reads displayed.

GTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGT  >  NC_000913/3287933‑3288131
                                                                                                   |                                                                                                   
gTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGAcc                                                                                                     <  1:1088141/100‑1 (MQ=255)
gTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGAcc                                                                                                     <  1:4169225/100‑1 (MQ=255)
 tATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCt                                                                                                    >  1:3695196/1‑100 (MQ=255)
 tATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCt                                                                                                    >  1:3695441/1‑100 (MQ=255)
  aTCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGAcc                                                                                                     <  2:2842048/98‑1 (MQ=255)
   tCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGt                                                                                                  >  1:2206570/1‑100 (MQ=255)
    cAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTc                                                                                                 >  1:1011205/1‑100 (MQ=255)
       ccGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGAccagatcgg                                                                                              >  1:2227988/1‑93 (MQ=255)
       ccGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGa                                                                                              <  1:450697/100‑1 (MQ=255)
       ccGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGa                                                                                              >  2:3871058/1‑100 (MQ=255)
       ccGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGa                                                                                              <  1:3747491/100‑1 (MQ=255)
          tAGTTCCCTCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTATCGTAACGCTGGATTACCGTTAATACGTTAGCGCGTTATCTGACAGGTCAGATaa                                                                                           <  2:2650475/100‑1 (MQ=255)
            gTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGATTACCGTTAATACGTTAGGGCGTTATCTGACAGGTCAGATAACGc                                                                                        <  2:3909807/100‑1 (MQ=38)
             ttCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACAGGTCAGATAACg                                                                                         <  2:1944067/99‑1 (MQ=255)
             ttCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACAGGTCAGATAACGc                                                                                        <  1:3275602/100‑1 (MQ=255)
              tCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGcc                                                                                       <  2:217048/100‑1 (MQ=255)
              tCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACAGGTCAGATAACGcc                                                                                       <  2:1379914/100‑1 (MQ=255)
               ccATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGccc                                                                                      <  1:3589861/100‑1 (MQ=255)
                cATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCt                                                                                     <  1:807532/100‑1 (MQ=255)
                cATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCt                                                                                     >  1:3797272/1‑100 (MQ=255)
                cATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCt                                                                                     >  1:2610375/1‑100 (MQ=255)
                cATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCt                                                                                     >  1:536365/1‑100 (MQ=255)
                 aTCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCtt                                                                                    >  1:3015769/1‑100 (MQ=255)
                 aTCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCtt                                                                                    >  2:2709320/1‑100 (MQ=255)
                     gAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTttcc                                                                                >  1:2177014/1‑100 (MQ=255)
                      aaaTCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTttcct                                                                               <  2:3814904/100‑1 (MQ=255)
                       aaTCACCCCAGGGGAAGCAGAGGCTATCGTGAACGTAACGCTGGATGACCGTTAATACGTGAGGGCGTTATCTGACAGGTCAGATAACGCCCtaacg                                                                                 <  2:4181386/97‑5 (MQ=25)
                         tCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACAGGTCAGATAACGCCCtaacgtatta                                                                            <  1:2028413/100‑10 (MQ=25)
                         tCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCATTATCTGACCTGTCAGATAACGCCCTTttccttcc                                                                            >  2:1538134/1‑100 (MQ=255)
                           aCCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCctc                                                                          <  1:1648843/100‑1 (MQ=255)
                            ccccAGGGGAAGCAGACGCTATCGTTAACGTCACGCTAGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCctct                                                                         <  2:674117/100‑1 (MQ=255)
                            ccccAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCctct                                                                         <  2:1342326/100‑1 (MQ=255)
                            ccccAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCctct                                                                         >  2:234228/1‑100 (MQ=255)
                            ccccAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGACCTTTTCCTTCctct                                                                         >  1:1860572/1‑100 (MQ=255)
                               cAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTtc                                                                      >  1:281813/1‑100 (MQ=255)
                               cAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTtc                                                                      >  1:3053208/1‑100 (MQ=255)
                                 ggggAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTtctc                                                                    >  2:2403395/1‑100 (MQ=255)
                                  gggAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCg                                                                   >  1:1854405/1‑100 (MQ=255)
                                   ggAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTTACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCg                                                                   <  2:2242676/99‑1 (MQ=255)
                                    gAAGCAGACGATATCGTTAACGTAACGATGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATCACGCCCTTTTCCTTCCTCTTTCTCGtt                                                                 <  2:1659018/100‑1 (MQ=255)
                                     aaGCAGGCGCTATCGTTAACGTAACGGTGGATTACAGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTg                                                                <  2:4011196/100‑1 (MQ=255)
                                       gCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTa                                                              <  1:1058245/100‑1 (MQ=255)
                                       gCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTa                                                              >  1:294473/1‑100 (MQ=255)
                                       gCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTa                                                              >  1:6767/1‑100 (MQ=255)
                                           aCGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAg                                                          <  2:3023021/100‑1 (MQ=255)
                                             gCTATCGTTAACGTAACGCTTGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGt                                                        <  2:1831845/100‑1 (MQ=255)
                                             gCTATCGTTAACGTAACGCTGGATTACCGTTAATTCGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGt                                                        <  2:3949528/100‑1 (MQ=255)
                                             gCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGt                                                        <  2:1303614/100‑1 (MQ=255)
                                             gCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGt                                                        <  2:440427/100‑1 (MQ=255)
                                                aTCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGa                                                     >  1:1271739/1‑100 (MQ=255)
                                                   gTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGaaaa                                                  >  2:696910/1‑100 (MQ=255)
                                                   gTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGaaaa                                                  >  1:1684016/1‑100 (MQ=255)
                                                        cgtaacgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAAc                                             <  1:2865003/100‑1 (MQ=255)
                                                         gtaacgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATagat                                            >  1:1004334/1‑97 (MQ=255)
                                                         gtaacgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACt                                            >  1:3460882/1‑100 (MQ=255)
                                                           aacgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTAt                                          <  2:3695175/100‑1 (MQ=255)
                                                            acgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATg                                         <  1:78929/100‑1 (MQ=255)
                                                             cgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGt                                        >  1:703546/1‑100 (MQ=255)
                                                                 ggATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCaaa                                    >  1:994628/1‑100 (MQ=255)
                                                                   aTTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAAc                                  >  1:1511850/1‑100 (MQ=255)
                                                                    ttACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACg                                 <  1:892063/100‑1 (MQ=255)
                                                                     tACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGa                                >  1:1458267/1‑100 (MQ=255)
                                                                       ccGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAAc                              >  1:3864145/1‑100 (MQ=255)
                                                                       ccGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAAc                              >  2:1430436/1‑100 (MQ=255)
                                                                        cGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACa                             >  1:2674226/1‑100 (MQ=255)
                                                                        cGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACa                             <  1:2428798/100‑1 (MQ=255)
                                                                         gTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAAAGCCCTTTTCCTTCCTCTTTCTCGTTGTAAAAAATTGAAAAATAACTATGTCAAAACGAACa                             <  2:118828/99‑1 (MQ=255)
                                                                          ttAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACAtt                           >  1:1601524/1‑100 (MQ=255)
                                                                           tAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTc                          >  2:3508381/1‑100 (MQ=255)
                                                                             gtttcatAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTcgc                        <  1:4163708/94‑1 (MQ=255)
                                                                                cGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGt                     >  1:2262857/1‑100 (MQ=255)
                                                                                cGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGt                     >  2:2548805/1‑100 (MQ=255)
                                                                                 gTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTg                    <  2:3243105/100‑1 (MQ=255)
                                                                                   tAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGat                  <  1:3833183/100‑1 (MQ=255)
                                                                                      ggCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATAtt               >  2:3359267/1‑100 (MQ=255)
                                                                                      ggCGTTATCTGACAGGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATAtt               >  2:1035251/1‑100 (MQ=255)
                                                                                      ggCGTTATCTGACAGGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATAtt               >  1:4027473/1‑100 (MQ=255)
                                                                                       gCGTTATCTGACCTGTCAGATAGCGCCCTTTTCCTTCCTCTTTCTCGTTGTATTAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGttatta              >  2:2433348/1‑100 (MQ=255)
                                                                                       gCGTTATCTGACAGGTCAGATAACGCCCTTTTCCTTTATCTTTCTCGTTGGATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGGGGTGATATTa              >  2:751555/1‑100 (MQ=38)
                                                                                       gCGTTATCTGACAGGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTa              >  1:450584/1‑100 (MQ=255)
                                                                                       gCGTTATCTGACAGGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTa              >  1:596553/1‑100 (MQ=255)
                                                                                       gCGTTATCTGACAGGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTa              >  2:2092979/1‑100 (MQ=255)
                                                                                       gCGTTATCTGACAGGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTa              >  2:3706178/1‑100 (MQ=255)
                                                                                        cGTTATCTGACAGGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTaa             >  1:2981012/1‑100 (MQ=255)
                                                                                         gTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAAc            >  1:4094231/1‑100 (MQ=255)
                                                                                         gTTATCTGACAGGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTTTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCCCGGTGATATTAAc            >  2:3863543/1‑100 (MQ=38)
                                                                                          ttATCTGACCTGTCATATAACGCCCTTTTACTACATCTTTCTCGTTGTATTAGATTGAAAAATAACTATGTCAAAACGCACATTCGCGGTGATATTAAcc           <  2:345416/100‑1 (MQ=255)
                                                                                          ttATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAAcc           <  1:2772091/100‑1 (MQ=255)
                                                                                          ttATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAAcc           >  1:4132407/1‑100 (MQ=255)
                                                                                          ttATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAAcc           <  2:1596755/100‑1 (MQ=255)
                                                                                            aTCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCtt         >  2:1759500/1‑100 (MQ=255)
                                                                                            aTCTGACCTGTCAGATAACGCACTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCtt         <  2:947605/100‑1 (MQ=255)
                                                                                             tCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCttg        >  1:1967656/1‑100 (MQ=255)
                                                                                               tGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCttgtt      >  1:1533341/1‑100 (MQ=255)
                                                                                               tGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCttgtt      >  2:1514328/1‑100 (MQ=255)
                                                                                                gACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCttgttg     <  1:1754590/100‑1 (MQ=255)
                                                                                                gACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCttgttg     <  1:877828/100‑1 (MQ=255)
                                                                                                gACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATAAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCttgttg     <  2:3107184/100‑1 (MQ=255)
                                                                                                  ccTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAACCGAACATTCGCGGTGATATTAACCTTGTtgtg   >  2:611516/1‑100 (MQ=255)
                                                                                                   cTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTgtgt  >  2:1828034/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGT  >  NC_000913/3287933‑3288131

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