Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,449,615 C→G 16.7% G145G (GGG→GGC) ‡ yfcO ← DUF2544 family putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,449,6150CG16.7% 503.2 / 85.7 192G145G (GGG→GGC) yfcODUF2544 family putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (79/80);  new base G (16/16);  total (95/97)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01

Only 100 of 202 total aligned reads displayed.

AGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACC  >  NC_000913/2449517‑2449714
                                                                                                   |                                                                                                   
aGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGcc                                                                                                    >  2:1848392/1‑100 (MQ=255)
  ggTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGcccc                                                                                                  <  2:1815633/100‑1 (MQ=255)
  ggTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGccc                                                                                                   <  1:4113040/99‑1 (MQ=255)
      aGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGAt                                                                                              <  1:1277050/100‑1 (MQ=255)
      aGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGAt                                                                                              <  1:1985417/100‑1 (MQ=255)
        cAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATAc                                                                                            <  1:3121240/100‑1 (MQ=255)
           ttttGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATagtc                                                                                          <  2:3104825/100‑4 (MQ=38)
            tttGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATagtca                                                                                         <  2:1731152/100‑5 (MQ=38)
            tttGACTTACTGACATAATTCTCATTCCATGAAACGTAATAAGGTTGAACAGAAACGGATGTACAAATCGTGTCGTGAC‑TATCAGGTCCCTGATACGtaa                                                                                        <  2:2591673/100‑3 (MQ=255)
             ttGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCAc                                                                                       <  1:1069676/100‑1 (MQ=38)
              tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCACg                                                                                      <  1:1698920/100‑1 (MQ=38)
              tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACg                                                                                      >  1:1551149/1‑100 (MQ=255)
               gACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGc                                                                                     >  1:2033564/1‑100 (MQ=255)
                aCTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACg                                                                                      <  1:480973/98‑1 (MQ=255)
                aCTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGcc                                                                                    <  1:2735294/100‑1 (MQ=255)
                 cTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGccc                                                                                   >  2:1169861/1‑100 (MQ=255)
                 cTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGcc                                                                                    <  1:3945402/99‑1 (MQ=255)
                  ttACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCACgaca                                                                                   <  1:1967592/99‑4 (MQ=37)
                   tACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCACgacacg                                                                                 <  1:2515884/100‑6 (MQ=25)
                   tACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCACgacacg                                                                                 <  1:2005300/100‑6 (MQ=25)
                    aCTGACATAATTCTCATTCCATGAAACGTCTTAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCACgacac                                                                                  <  2:1041865/98‑5 (MQ=25)
                       gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGga                                                                             >  1:802565/1‑100 (MQ=255)
                       gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCCTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGga                                                                             >  2:772365/1‑100 (MQ=255)
                         cATAATTCTCATTCCATGAAACGTCATAAGGTTGAACATAAACGGATGGACAAATCGTCTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGgaga                                                                           <  2:2285505/100‑1 (MQ=255)
                            aaTTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGc                                                                        >  2:2494146/1‑100 (MQ=255)
                                ctcATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCa                                                                    <  1:3070602/100‑1 (MQ=255)
                                ctcATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCCTGTCGTGAC‑TATCAGGCCGCTGATACGTCACGCGCCCTGGAGAGGCtt                                                                      >  2:1145093/1‑98 (MQ=255)
                                 tcATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCag                                                                   >  1:2095971/1‑100 (MQ=255)
                                 tcATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCag                                                                   <  1:3139250/100‑1 (MQ=255)
                                  cATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCaga                                                                  >  1:597311/1‑100 (MQ=255)
                                   aTTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCagag                                                                 <  1:3568510/100‑1 (MQ=255)
                                   aTTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCagag                                                                 <  1:2126370/100‑1 (MQ=255)
                                   aTTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCagag                                                                 >  1:1423251/1‑100 (MQ=255)
                                    ttCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGt                                                                <  2:730835/100‑1 (MQ=255)
                                      ccATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtca                                                              >  2:220751/1‑100 (MQ=255)
                                      ccATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtca                                                              >  2:2417707/1‑100 (MQ=255)
                                        aTGAACCGTTATCAGGTTGAACAGAAACGGAGGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATGCGTTACGCACCCTGGAGAGGCTTCAGAGtcat                                                             >  2:3990652/1‑99 (MQ=255)
                                        aTGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcatt                                                            >  1:2113393/1‑100 (MQ=255)
                                         tGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcattc                                                           >  1:2439243/1‑100 (MQ=255)
                                         tGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcattc                                                           >  2:1968363/1‑100 (MQ=255)
                                          gAAACGTCATAAGGTTGAACAGCAACGGATGTACAAATCGTGTCGTGTC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcattca                                                          <  2:2830242/100‑1 (MQ=255)
                                           aaaCGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcattcat                                                         <  1:2895584/100‑1 (MQ=255)
                                             aCGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCataa                                                       >  1:179957/1‑100 (MQ=255)
                                              cGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCataat                                                      >  2:427086/1‑100 (MQ=255)
                                               gTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCataata                                                     <  2:2699191/100‑1 (MQ=255)
                                               gTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCataata                                                     <  2:3538335/100‑1 (MQ=255)
                                                tCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCAtaataa                                                    >  1:2287381/1‑100 (MQ=255)
                                                  aTAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaacc                                                  <  1:219704/100‑1 (MQ=255)
                                                   tAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaacca                                                 >  2:3716213/1‑100 (MQ=255)
                                                   tAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaacca                                                 >  1:1588386/1‑100 (MQ=255)
                                                    aaGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaaccaa                                                >  2:2247833/1‑100 (MQ=255)
                                                      ggTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaaccaacc                                              >  1:133290/1‑100 (MQ=255)
                                                       gTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCataatagatcggaa                                             >  2:1567479/1‑92 (MQ=255)
                                                        ttGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCAt                                            <  2:1252560/100‑1 (MQ=255)
                                                        ttGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCAt                                            >  1:1176571/1‑100 (MQ=255)
                                                         tGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATg                                           <  1:1836507/100‑1 (MQ=255)
                                                          gAACAGAAACTGATGTACAAATCGTGTCGTGAC‑TATCAGGCCCCTTATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATgg                                          >  2:616441/1‑100 (MQ=255)
                                                          gAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATgg                                          <  1:275009/100‑1 (MQ=255)
                                                           aaCAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGa                                         <  1:1863922/100‑1 (MQ=255)
                                                           aaCAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGa                                         <  2:1527157/100‑1 (MQ=255)
                                                           aaCAGAAACGGATGCACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGa                                         <  2:2858189/100‑1 (MQ=255)
                                                            aCAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGAc                                        <  2:2114071/100‑1 (MQ=255)
                                                            aCAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGAc                                        <  2:1611304/100‑1 (MQ=255)
                                                             cAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACa                                       >  1:3030090/1‑100 (MQ=255)
                                                             cAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACa                                       <  2:1691687/100‑1 (MQ=255)
                                                                 aaCGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGaa                                   <  1:3459160/100‑1 (MQ=255)
                                                                   cGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGaaaa                                 <  1:2564960/100‑1 (MQ=255)
                                                                   cGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGaaaa                                 >  1:3991097/1‑100 (MQ=255)
                                                                    ggATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAg                                <  1:748317/100‑1 (MQ=255)
                                                                      aTGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTa                              >  1:3395007/1‑100 (MQ=255)
                                                                         gACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGa                           >  1:1006863/1‑100 (MQ=255)
                                                                         gACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCATGGAGAGGCTTCAGAGTCATTCAACATACCCCACCATGGACAAGAAAAGTAAg                            >  2:1489586/1‑99 (MQ=255)
                                                                          aCAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGAt                          <  1:3716234/100‑1 (MQ=255)
                                                                           cAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGAtt                         >  1:384715/1‑100 (MQ=255)
                                                                           cAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGAtt                         >  2:1988760/1‑100 (MQ=255)
                                                                           cAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGAtt                         <  2:668242/100‑1 (MQ=255)
                                                                           cAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGAt                          >  2:3021678/1‑99 (MQ=255)
                                                                               tcgtgtcgtgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAt                     <  1:377783/100‑1 (MQ=255)
                                                                                cgtgtcgtgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtc                    <  1:3351271/100‑1 (MQ=255)
                                                                                 gtgtcgtgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtct                   >  2:412255/1‑100 (MQ=255)
                                                                                   ggcgtgACGTATCAGGGGCCTGATA‑GTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATctct                 >  1:910866/3‑100 (MQ=37)
                                                                                    tcgtgAC‑TATGAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCCACTATGGACAAGAAAAGTAAGATTCAATCTCTg                <  2:2757795/100‑1 (MQ=255)
                                                                                     cgtgACGTATCAGGGGCCTGATA‑GTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc               >  1:727377/1‑100 (MQ=38)
                                                                                     cgtgACGTATCAGGGGCCTGATA‑GTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc               >  1:300978/1‑100 (MQ=38)
                                                                                     cgtgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc               >  1:207183/1‑100 (MQ=255)
                                                                                      gtgacgtATCAGGGGCCTGATA‑GTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCa              >  2:4015125/7‑100 (MQ=38)
                                                                                      gtgacgtATCAGGGGCCTGATA‑GTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCa              >  1:3126850/7‑100 (MQ=38)
                                                                                      gtgAC‑TATCAGGCCCGTGATACGTCACGCAGCCTGGAGAGGCTTCAGAGTCATTCATGATAACCAACCATGGACAAGAAAAGTAAGATTCAATctct                 <  2:3363505/97‑1 (MQ=255)
                                                                                       tgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCag             <  2:3943304/100‑1 (MQ=255)
                                                                                          c‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCagatc          >  2:3987059/1‑98 (MQ=255)
                                                                                          c‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAgaga          >  1:2800008/1‑100 (MQ=255)
                                                                                          c‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAgaga          >  1:3933605/1‑100 (MQ=255)
                                                                                            tATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGAc         >  2:2181664/1‑100 (MQ=255)
                                                                                                nggCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCACTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATa     >  2:78210/2‑100 (MQ=255)
                                                                                                aGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCagatcggaag     >  2:209811/1‑93 (MQ=255)
                                                                                                aGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATa     >  1:2935071/1‑100 (MQ=255)
                                                                                                 ggCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATaa    >  1:4052097/1‑100 (MQ=255)
                                                                                                 ggCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATaa    >  2:270986/1‑100 (MQ=255)
                                                                                                 ggCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATa     <  2:2641316/99‑1 (MQ=255)
                                                                                                   ccccTGATACGTCACGCCCCCTGGAGAAGCTTCACAGTCATTCATAATAACCAACCATGCACAAGAAAAGCAAGATTCCATCTCTGCAGAGACAATAAcc  <  2:2033969/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
AGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACC  >  NC_000913/2449517‑2449714

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: ‑