Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 966,368 A→G 100% A683A (GCA→GCG ycaI → ComEC family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913966,3680AG99.1% 357.9 / ‑4.9 102A683A (GCA→GCGycaIComEC family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (49/52);  total (50/52)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 102 total aligned reads displayed.

GGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCC  >  NC_000913/966271‑966463
                                                                                                 |                                                                                               
ggCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCt                                                                                               <  1:198134/100‑1 (MQ=255)
ggCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACt                                                                                               >  2:2152169/1‑100 (MQ=255)
 gcagcaAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTc                                                                                              <  1:2388979/100‑1 (MQ=255)
  cagcaAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCg                                                                                             >  2:118388/1‑100 (MQ=255)
   agcaAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATAGAAAAGTGGCGCTCgc                                                                                            >  2:1497711/1‑100 (MQ=255)
     caAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgcat                                                                                          >  1:1647452/1‑100 (MQ=255)
      aaGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgcatc                                                                                         <  1:298928/100‑1 (MQ=255)
         tGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgcatcggc                                                                                      <  2:2514104/100‑1 (MQ=255)
         tGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgcatcggc                                                                                      <  1:2362285/100‑1 (MQ=255)
         tGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgcatcggc                                                                                      <  1:1888705/100‑1 (MQ=255)
         tGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgcatcgg                                                                                       <  2:616919/99‑1 (MQ=255)
           cAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgcatcggcat                                                                                    >  2:1413685/1‑100 (MQ=255)
            aGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgcatcggca                                                                                     <  2:454853/98‑1 (MQ=255)
             ggCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgcatcggcatcg                                                                                  >  2:1478857/1‑100 (MQ=255)
              gCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATcgc                                                                                 >  2:1706786/1‑100 (MQ=255)
              gCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCACATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATcgc                                                                                 <  2:2391375/100‑1 (MQ=255)
               cAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATcgcg                                                                                >  1:1213204/1‑100 (MQ=255)
                  cATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCta                                                                             <  2:1755710/100‑1 (MQ=255)
                  cATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCta                                                                             >  1:1526958/1‑100 (MQ=255)
                  cATTGCGTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCAGTAATTCAGAGAGTGAATGGAAAAGGGGCGCTCGCATCGGCATCGCGCta                                                                             <  2:2256473/100‑1 (MQ=255)
                     tGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATaa                                                                          >  1:925862/1‑100 (MQ=255)
                      gCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAAc                                                                         >  1:1048016/1‑100 (MQ=255)
                       cTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACg                                                                        <  1:1191311/100‑1 (MQ=255)
                       cTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACg                                                                        <  1:67374/100‑1 (MQ=255)
                       cTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACg                                                                        >  1:568609/1‑100 (MQ=255)
                             gTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGc                                                                  >  2:15155/1‑100 (MQ=255)
                               aCCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGa                                                                >  1:1956492/1‑100 (MQ=255)
                                ccTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGAc                                                               >  2:1122646/1‑100 (MQ=255)
                                ccTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGAc                                                               <  2:1647440/100‑1 (MQ=255)
                                 cTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACt                                                              >  1:2008345/1‑100 (MQ=255)
                                 cTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACt                                                              <  2:1506457/100‑1 (MQ=255)
                                  tCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTg                                                             >  1:712508/1‑100 (MQ=255)
                                   cACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGc                                                            >  2:1150808/1‑100 (MQ=255)
                                    aCCATGGGAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGc                                                            <  2:740688/99‑1 (MQ=255)
                                    aCCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGTATCGGCATCGCGCTATAACGCATGGCGACTGcc                                                           >  2:233528/1‑100 (MQ=255)
                                    aCCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGcc                                                           <  2:1589315/100‑1 (MQ=255)
                                     ccATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGccc                                                          <  2:1438790/100‑1 (MQ=255)
                                     ccATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGccc                                                          >  1:1648250/1‑100 (MQ=255)
                                     ccATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGccc                                                          >  2:580063/1‑100 (MQ=255)
                                       aTGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCctc                                                        >  2:1950665/1‑100 (MQ=255)
                                       aTGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCctc                                                        <  2:2495241/100‑1 (MQ=255)
                                       aTGGCAGTAATACCACATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCctc                                                        <  2:1649872/100‑1 (MQ=255)
                                         ggCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTa                                                      <  1:2484248/100‑1 (MQ=255)
                                         ggCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTa                                                      <  1:580070/100‑1 (MQ=255)
                                         ggCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTa                                                      >  1:1388880/1‑100 (MQ=255)
                                          gCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTaa                                                     <  2:1435636/100‑1 (MQ=255)
                                           cAGTAATACCACATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAAc                                                    >  2:1368367/1‑100 (MQ=255)
                                             gTAATACCTCATCAGCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACaa                                                  >  2:242332/1‑100 (MQ=255)
                                               aaTACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAg                                                <  1:401897/100‑1 (MQ=255)
                                               aaTACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAg                                                <  1:1356161/100‑1 (MQ=255)
                                                aTACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGt                                               >  1:1603590/1‑100 (MQ=255)
                                                aTACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGt                                               >  2:1515862/1‑100 (MQ=255)
                                                  aCCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTa                                             <  1:1478868/100‑1 (MQ=255)
                                                  aCCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTa                                             <  2:1460467/100‑1 (MQ=255)
                                                   ccTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGCAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTaa                                            >  2:1974723/1‑100 (MQ=255)
                                                   ccTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTaa                                            <  2:116958/100‑1 (MQ=255)
                                                   ccTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTaa                                            <  2:116959/100‑1 (MQ=255)
                                                   ccTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTaa                                            >  2:1176976/1‑100 (MQ=255)
                                                   ccTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTaa                                            >  1:730843/1‑100 (MQ=255)
                                                    cTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAg                                           >  1:2004157/1‑100 (MQ=255)
                                                    cTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAg                                           <  2:455458/100‑1 (MQ=255)
                                                     tcatcaTCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGc                                          <  1:291599/100‑1 (MQ=255)
                                                      catcatCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCa                                         <  2:185649/100‑1 (MQ=255)
                                                        tcatcGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATc                                       >  2:648820/1‑100 (MQ=255)
                                                         catcGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCg                                      >  1:2200551/1‑100 (MQ=255)
                                                            cGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTa                                   <  2:2284194/100‑1 (MQ=255)
                                                            cGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTa                                   >  1:949738/1‑100 (MQ=255)
                                                               tGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCa                                <  1:2228987/100‑1 (MQ=255)
                                                               tGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCa                                <  1:1928098/100‑1 (MQ=255)
                                                                 ccATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAAc                              >  2:1978101/1‑100 (MQ=255)
                                                                 ccATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAAc                              <  2:2060211/100‑1 (MQ=255)
                                                                  ccttAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACt                             >  1:1740989/3‑100 (MQ=255)
                                                                   aTTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTg                            >  2:191130/1‑100 (MQ=255)
                                                                    ttAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGc                           <  2:316928/100‑1 (MQ=255)
                                                                    ttAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGc                           <  2:820991/100‑1 (MQ=255)
                                                                      aaTTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCaa                         >  1:1035279/1‑100 (MQ=255)
                                                                       aTTCAGCGAGTGAATGGGAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAg                        <  1:516112/100‑1 (MQ=255)
                                                                       aTTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAg                        <  2:1123005/100‑1 (MQ=255)
                                                                        ttCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCCCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAgg                       >  1:1255042/1‑100 (MQ=255)
                                                                         tCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGa                      <  2:502741/100‑1 (MQ=255)
                                                                           aGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGcta                    >  2:864526/1‑100 (MQ=255)
                                                                           aGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGata                    >  1:1454908/1‑100 (MQ=255)
                                                                           aGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGata                    <  1:532989/100‑1 (MQ=255)
                                                                               aGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCaa                >  1:142544/1‑100 (MQ=255)
                                                                                 tGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATg              <  1:450869/100‑1 (MQ=255)
                                                                                    aTGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGAt           >  1:436812/1‑100 (MQ=255)
                                                                                    aTGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGAt           >  2:2096198/1‑100 (MQ=255)
                                                                                    aTGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGAt           <  1:1413696/100‑1 (MQ=255)
                                                                                    aTGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGAt           >  1:752597/1‑100 (MQ=255)
                                                                                     tGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACCAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGAtt          <  2:1401150/100‑1 (MQ=255)
                                                                                     tGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGAtt          >  1:2032043/1‑100 (MQ=255)
                                                                                     tGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGAtt          <  1:607782/100‑1 (MQ=255)
                                                                                     tGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGAtt          <  2:45811/100‑1 (MQ=255)
                                                                                       gAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGa        <  1:965066/100‑1 (MQ=255)
                                                                                       gAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGa        >  1:2381396/1‑100 (MQ=255)
                                                                                       gAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGa        <  2:2387492/100‑1 (MQ=255)
                                                                                        aaaaGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGAt       <  1:191132/100‑1 (MQ=255)
                                                                                         aaaGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATa      >  1:1177175/1‑100 (MQ=255)
                                                                                         aaaGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATa      <  2:1255034/100‑1 (MQ=255)
                                                                                             tGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTcc  >  2:2153781/1‑100 (MQ=255)
                                                                                                 |                                                                                               
GGCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCC  >  NC_000913/966271‑966463

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