Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,549,932 T→A 17.5% intergenic (+222/+21) fimH → / ← gntP minor component of type 1 fimbriae/fructuronate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,549,9320TA17.5% 327.8 / 55.3 136intergenic (+222/+21)fimH/gntPminor component of type 1 fimbriae/fructuronate transporter
Reads supporting (aligned to +/- strand):  ref base T (60/52);  new base A (15/9);  total (75/61)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 137 total aligned reads displayed.

CTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTT  >  NC_000913/4549833‑4550028
                                                                                                   |                                                                                                
cTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCAt                                                                                                  <  1:1805850/100‑1 (MQ=255)
    gCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGc                                                                                              >  2:1864777/1‑100 (MQ=255)
     cTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCa                                                                                             >  1:818089/1‑100 (MQ=255)
     cTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCa                                                                                             >  1:2310923/1‑100 (MQ=255)
     cTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCa                                                                                             <  1:1154414/100‑1 (MQ=255)
     cTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCa                                                                                             >  1:875764/1‑100 (MQ=255)
     cTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCa                                                                                             <  1:549305/100‑1 (MQ=255)
      tGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAg                                                                                            >  1:1160964/1‑100 (MQ=255)
      tGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAg                                                                                            >  1:2577672/1‑100 (MQ=255)
       ggTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAgg                                                                                           <  1:1869535/100‑1 (MQ=255)
       ggTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACAGTTTTAGACTTCGCCATCAGTCTGAGCCCGCCATTGGCAgg                                                                                           >  1:812685/1‑100 (MQ=25)
        gTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGCGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGGGCCTGCCATTGGCAGGc                                                                                          <  1:325118/100‑1 (MQ=255)
          cTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGTGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCAATGGCAGGctc                                                                                        <  1:2612161/100‑1 (MQ=255)
          cTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCAATGGCAGGctc                                                                                        <  2:2164421/100‑1 (MQ=255)
          cTGTTCCCCTCACCCTAACCCGCTCCCCTGAGCGGCGAGAGGACTTTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGTCTGCCATTGGCAgg                                                                                           <  1:378562/97‑1 (MQ=37)
           tgtgCCCCTCACCCTAACCCTCTCCCCGGAGTGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCAATGGCAGGctca                                                                                       <  1:2864006/100‑2 (MQ=37)
            gTTCCCCTCACCCTAACCCTCTCCCCGGGGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCAGAGCCTGCCAATGGCAGGctcag                                                                                      <  2:2314905/100‑3 (MQ=25)
              tCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGgt                                                                                    <  1:26222/100‑1 (MQ=255)
              tCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGgt                                                                                    <  2:599737/100‑1 (MQ=255)
              tCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGgt                                                                                    >  2:1405986/1‑100 (MQ=255)
               ccccTCCCCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCAATGGCAGGctcagact                                                                                   <  1:2014582/100‑6 (MQ=37)
               ccccTCACCCTACCCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGgtg                                                                                   <  1:218560/100‑1 (MQ=255)
                cccTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGgtgt                                                                                  >  1:1599449/1‑100 (MQ=255)
                cccTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGgt                                                                                    >  1:2773680/1‑98 (MQ=255)
                   tCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGCCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCCGGCTCTGGCGTcc                                                                                >  1:1691650/1‑99 (MQ=255)
                      cccTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCtttt                                                                            >  2:2466382/1‑100 (MQ=255)
                      cccTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCtttt                                                                            >  1:1220597/1‑100 (MQ=255)
                         tAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGTGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACg                                                                         <  1:1196630/100‑1 (MQ=255)
                           aCCCTCGCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCt                                                                       <  1:555825/100‑1 (MQ=255)
                            cccTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTa                                                                      >  1:307681/1‑100 (MQ=255)
                            cccTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTa                                                                      >  1:1276755/1‑100 (MQ=255)
                               tctcCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGTCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCa                                                                   <  2:499031/100‑1 (MQ=255)
                                ctcCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTCAGCCTGCCATCGGCAGGCTCTGGTGTCCCTTTACGCTACCAt                                                                  <  2:1969651/100‑1 (MQ=255)
                                 tcCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCAt                                                                  >  2:1073690/1‑99 (MQ=255)
                                  ccccGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGc                                                                >  2:2359252/1‑100 (MQ=255)
                                  ccccGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGc                                                                >  1:257754/1‑100 (MQ=255)
                                   cccGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCt                                                               <  1:665073/100‑1 (MQ=255)
                                   cccGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCt                                                               <  2:2197725/100‑1 (MQ=255)
                                      ggAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCGTCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCtaat                                                            <  1:1864768/100‑1 (MQ=255)
                                      ggAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCtaat                                                            >  2:2391986/1‑100 (MQ=255)
                                      ggAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCtaat                                                            <  2:14359/100‑1 (MQ=255)
                                       gAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCtaata                                                           <  1:738796/100‑1 (MQ=255)
                                       gAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCtaata                                                           <  1:1496707/100‑1 (MQ=255)
                                           ggCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCa                                                       <  2:257755/100‑1 (MQ=255)
                                           ggCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCa                                                       >  2:1471479/1‑100 (MQ=255)
                                            gCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAg                                                      >  1:2601477/1‑100 (MQ=255)
                                             cGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGc                                                     <  2:1801232/100‑1 (MQ=255)
                                             cGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGc                                                     <  2:1622439/100‑1 (MQ=255)
                                               aGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGcac                                                   <  2:929716/100‑1 (MQ=255)
                                                ggggACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGcaca                                                  <  2:2111677/100‑1 (MQ=255)
                                                  ggACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACaat                                                >  1:437030/1‑100 (MQ=255)
                                                      tGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATc                                            <  2:2196810/100‑1 (MQ=255)
                                                       gtacGGGCACATTTTTAGACTTTGCGATTGGTCTAAGCCTCCCATTGGAAGACGTTGTTTTCCTTTTACGAAAACATTCTAACAATCAGCCCaataacc                                            >  1:1982945/4‑97 (MQ=255)
                                                       gTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCa                                           >  2:1068595/1‑100 (MQ=255)
                                                       gTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCa                                           >  2:1752156/1‑100 (MQ=255)
                                                         ccGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGc                                         <  2:2573599/100‑1 (MQ=255)
                                                          cgtgCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGcc                                        <  2:682610/97‑1 (MQ=255)
                                                          cGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGcc                                        >  2:1095666/1‑100 (MQ=255)
                                                          cGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGcc                                        >  1:1155521/1‑100 (MQ=255)
                                                           gggCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGccc                                       >  1:1832295/1‑100 (MQ=255)
                                                             gCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCaa                                     >  2:117256/1‑100 (MQ=255)
                                                              cacaTTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCatc                                    <  2:1966716/100‑3 (MQ=255)
                                                                caTTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCa                                  <  2:2844782/100‑1 (MQ=255)
                                                                caTTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCa                                  <  2:960569/100‑1 (MQ=255)
                                                                caTTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCa                                  >  1:892239/1‑100 (MQ=255)
                                                                 aTTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTAACATGCTAATAATCAGCACAATAATCAGCCCAACCAc                                 >  2:2777152/1‑100 (MQ=255)
                                                                  tttttAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACg                                >  2:2103030/1‑100 (MQ=255)
                                                                   ttttAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCCGCACAATAATCAGCCCAACCACg                                >  1:982749/1‑99 (MQ=255)
                                                                    tttAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACgg                               <  1:1816790/99‑1 (MQ=255)
                                                                    tttAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGa                              >  1:2477247/1‑100 (MQ=255)
                                                                        gACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTg                          >  1:817538/1‑100 (MQ=255)
                                                                         aCTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGa                         <  2:2587756/100‑1 (MQ=255)
                                                                          cTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCCCGGAGTTGAc                        >  1:2032430/1‑100 (MQ=255)
                                                                          cTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGAc                        >  1:1812619/1‑100 (MQ=255)
                                                                           tttGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGAcc                       <  2:1277263/100‑1 (MQ=255)
                                                                           tttGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGAcc                       <  2:806196/100‑1 (MQ=255)
                                                                           tttGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGAcc                       <  2:2120264/100‑1 (MQ=255)
                                                                             tGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCCTGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAg                     <  1:1611472/100‑1 (MQ=255)
                                                                               tcatcaGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCt                   >  2:2455130/1‑100 (MQ=255)
                                                                                catcaGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTc                  >  1:1564301/1‑100 (MQ=255)
                                                                                 atcaGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTc                  >  1:2038990/1‑99 (MQ=255)
                                                                                  tcaGTCTGAGCCTGCCATTGGCAGGCTCTGGTGCCCTTTTACACTACCATGCTAATAATCAGAACAATAATTAGCCCAACCACGGAGTGGACCAGcagca                >  1:1390815/1‑96 (MQ=255)
                                                                                  gacaccagAGCCTGCCAATGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCca                >  2:1989734/8‑100 (MQ=39)
                                                                                   acaccagAGCCTGCCAATGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCcag               >  2:1547110/7‑100 (MQ=255)
                                                                                     gTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCcagca             <  1:1471473/100‑1 (MQ=255)
                                                                                     gTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCcagca             >  1:1701903/1‑100 (MQ=255)
                                                                                     accagAGCCTGCCAATGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCcagca             >  1:178059/5‑100 (MQ=255)
                                                                                     accagAGCCTGCCAATGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCcagca             >  2:202680/5‑100 (MQ=255)
                                                                                      ccagAGCCTGCCAATGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCcagcag            >  2:577114/4‑100 (MQ=255)
                                                                                      ccagAGCCTGCCAATGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCcagcag            >  1:2654829/4‑100 (MQ=255)
                                                                                       nagAGCCTGGCAATGGCAGGCAATGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATAAGACCAACCACGGAGTTGACCAGCTCCAGCAGa           >  1:594835/3‑100 (MQ=255)
                                                                                       cagAGCCTGCCAATGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGa           >  1:1988985/3‑100 (MQ=255)
                                                                                        agAGCCTGCCAATGGCAGGCTCTGGTGTCCTTCTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGAc          >  2:2767353/2‑100 (MQ=255)
                                                                                           gCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGAcccc       >  1:2256796/1‑100 (MQ=255)
                                                                                             cTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAg     >  2:206152/1‑100 (MQ=255)
                                                                                              tGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTTACCAGCTCCAGCAGACCCCAgg    >  2:737698/1‑100 (MQ=255)
                                                                                              tGCCAATGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATGATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAgg    >  1:1862137/1‑100 (MQ=255)
                                                                                              tGCCAATGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAgg    >  1:364179/1‑100 (MQ=255)
                                                                                                ccATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGtt  >  2:366602/1‑100 (MQ=255)
                                                                                                ccATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGtt  <  2:583536/100‑1 (MQ=255)
                                                                                                   |                                                                                                
CTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTT  >  NC_000913/4549833‑4550028

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: