Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,452,339 A→G 24.0% intergenic (‑3/+17) yfcS ← / ← yfcU putative periplasmic pilin chaperone/pseudogene, FimD fimbrial export usher family;putative membrane; Not classified; putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,452,3390AG24.0% 220.4 / 76.5 129intergenic (‑3/+17)yfcS/yfcUputative periplasmic pilin chaperone/pseudogene, FimD fimbrial export usher family;putative membrane; Not classified; putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (42/55);  new base G (13/18);  total (56/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

Only 100 of 130 total aligned reads displayed.

GTTACAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACAGGCTTTATGTCAT  >  NC_000913/2452245‑2452438
                                                                                              |                                                                                                   
gTTACAGACGCCAGGGCAGACAAACGGGTAGCGCTAAGCAGTAACGCCCGCGTCATAGCCCCCAATTTGGCAGAGCAAAGTAAATCAGGCATGGCTCAtt                                                                                                >  2:3219708/1‑100 (MQ=255)
   aCAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCt                                                                                             <  2:1388874/100‑1 (MQ=255)
   aCAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCt                                                                                             <  2:3039483/100‑1 (MQ=255)
   aCAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCt                                                                                             >  1:426705/1‑100 (MQ=255)
     aGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTa                                                                                           >  1:1905173/1‑100 (MQ=255)
      gACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTaa                                                                                          >  1:254814/1‑100 (MQ=255)
      gACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTaa                                                                                          >  2:3065063/1‑100 (MQ=255)
        cGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATc                                                                                        >  1:1914142/1‑100 (MQ=255)
              ggCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTa                                                                                  >  2:2550970/1‑100 (MQ=255)
               gCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTAt                                                                                 >  1:366394/1‑100 (MQ=255)
                cAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAAGTTGGCAGAGCAAAGGAAATCAGACATGGATCATTCCTTAATCACGTTAtt                                                                                >  1:280883/1‑100 (MQ=255)
                 aGACAAACTGGTAGCGCTAAGCAGTAACGGCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAATTCAGACATGGATCATTCCTTAATCACGTTATTg                                                                               <  1:2037437/100‑1 (MQ=255)
                 aGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTg                                                                               <  1:2896429/100‑1 (MQ=255)
                 aGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTg                                                                               >  2:1352641/1‑100 (MQ=255)
                  gACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTgg                                                                              >  1:1970750/1‑100 (MQ=255)
                   aCAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGt                                                                             <  2:202728/100‑1 (MQ=255)
                    cAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGtt                                                                            <  2:669428/100‑1 (MQ=255)
                     aaaCTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTa                                                                           <  2:1451926/100‑1 (MQ=255)
                     aaaCTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTa                                                                           <  1:1148172/100‑1 (MQ=255)
                      aaCTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAg                                                                          <  1:405474/100‑1 (MQ=255)
                       aCTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGc                                                                         <  2:3119752/100‑1 (MQ=255)
                       aCTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGc                                                                         <  2:1035196/100‑1 (MQ=255)
                        cTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTCGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCt                                                                        <  2:1890798/100‑1 (MQ=255)
                         tGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTg                                                                       >  1:2200029/1‑100 (MQ=255)
                         tGGGAGCGCTAAGCAGTAACGCCAGCGCCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTg                                                                       <  1:270706/100‑1 (MQ=255)
                          ggTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTga                                                                      >  2:1376038/1‑100 (MQ=255)
                          ggTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTga                                                                      >  2:2327292/1‑100 (MQ=255)
                          ggTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTga                                                                      >  1:602015/1‑100 (MQ=255)
                          ggTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTga                                                                      <  2:2104621/100‑1 (MQ=255)
                           gTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTgaa                                                                     >  2:2609620/1‑100 (MQ=255)
                           gTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTgaa                                                                     <  2:280884/100‑1 (MQ=255)
                            tAGCGCTAAGCAGTAACGCCAACGTCATAGCACCCAATTTGGCAGACCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTgaaa                                                                    <  1:104515/100‑1 (MQ=255)
                              gcgcTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTgaaaca                                                                  >  1:2805428/1‑100 (MQ=255)
                              gcgcTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTgaaaca                                                                  >  1:2022935/1‑100 (MQ=255)
                               cgcTAAGCAGTAACGCGAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCACTCCTTAATCACGTTATTGGTTAGc                                                                         <  1:837516/92‑1 (MQ=255)
                               cgcTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTgaaacag                                                                 <  2:2772484/100‑1 (MQ=255)
                                gcTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTgaaacaga                                                                >  1:2193382/1‑100 (MQ=255)
                                 cTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTgaaacagaa                                                               >  1:2308539/1‑100 (MQ=255)
                                     gCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGaaacagaaacag                                                           >  1:2596/1‑100 (MQ=255)
                                      cAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAgg                                                          >  2:144856/1‑100 (MQ=255)
                                       ccgatctccAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAggg                                                         <  2:200673/93‑1 (MQ=255)
                                       aGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAggg                                                         <  2:1970763/100‑1 (MQ=255)
                                       aGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAggg                                                         >  1:2952878/1‑100 (MQ=255)
                                       aGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAggg                                                         >  1:2976430/1‑100 (MQ=255)
                                       aGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAggg                                                         >  2:2838701/1‑100 (MQ=255)
                                       aGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAggg                                                         <  2:2511212/100‑1 (MQ=255)
                                        gTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAgggg                                                        <  2:3124999/100‑1 (MQ=255)
                                        gTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAgggg                                                        <  2:3097814/100‑1 (MQ=255)
                                            cGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATc                                                    >  2:2561845/1‑100 (MQ=255)
                                               cAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAAGCAGACATGGGTCATTCCTTAATCACTTTATTGGTTAGCTGAAACAGAAACAGGGGGATCggt                                                 <  1:1802564/100‑1 (MQ=255)
                                                aGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCggtg                                                <  1:790732/100‑1 (MQ=255)
                                                  cGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCggtggg                                              <  2:602018/100‑1 (MQ=255)
                                                  cGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCggtggg                                              <  1:1275275/100‑1 (MQ=255)
                                                   gTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCggtgggt                                             >  1:1101323/1‑100 (MQ=255)
                                                    tCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCggtgggtg                                            >  1:970073/1‑100 (MQ=255)
                                                    tCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCggtgggtg                                            <  1:235211/100‑1 (MQ=255)
                                                    tCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCggtgggtg                                            <  2:2401599/100‑1 (MQ=255)
                                                      aTAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTc                                          <  2:2596/100‑1 (MQ=255)
                                                        aGCACCCAAGTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCCCGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCAc                                        <  2:1945926/100‑1 (MQ=255)
                                                         gCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCAcc                                       <  1:464286/100‑1 (MQ=255)
                                                         gCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCAcc                                       >  2:1847346/1‑100 (MQ=255)
                                                           aCCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTg                                     >  2:1010279/1‑100 (MQ=255)
                                                            cccAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGt                                    >  2:383157/1‑100 (MQ=255)
                                                            cccAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGt                                    <  2:1284997/100‑1 (MQ=255)
                                                              cAATTTGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTg                                  <  1:1376032/100‑1 (MQ=255)
                                                              cAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTg                                  <  1:864740/100‑1 (MQ=255)
                                                              cAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTg                                  <  2:1357019/100‑1 (MQ=255)
                                                                 tttGGCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGt                               <  1:1857173/100‑1 (MQ=255)
                                                                  ttGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTc                              >  1:1653245/1‑100 (MQ=255)
                                                                  ttGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTc                              <  2:1896540/100‑1 (MQ=255)
                                                                     gCAGAGCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGGCTGc                           >  1:1431308/1‑100 (MQ=255)
                                                                        gagCAAAGTAAATCAGACATGGATCATGCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTcc                        <  1:678272/100‑1 (MQ=255)
                                                                         agCAAAGTAAATCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCt                       <  2:1151166/100‑1 (MQ=255)
                                                                          gCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGATACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTg                      <  1:1847335/100‑1 (MQ=255)
                                                                          gCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTg                      >  1:2045428/1‑100 (MQ=255)
                                                                            aaaGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGaa                    >  1:2569016/1‑100 (MQ=255)
                                                                              aGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATg                  <  1:91231/100‑1 (MQ=255)
                                                                              aGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATg                  <  1:2561825/100‑1 (MQ=255)
                                                                              aGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATg                  >  2:1614367/1‑100 (MQ=255)
                                                                              aGTAAATCAGACATGGACCATTCCTTAATCACGTTATTTGTTAGCTGAAACAGAAACAGGGCGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATg                  <  1:1413948/100‑1 (MQ=255)
                                                                               gTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGa                 >  1:2168710/1‑100 (MQ=255)
                                                                                tAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGAc                <  1:267453/100‑1 (MQ=255)
                                                                                tAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGAc                <  1:2630752/100‑1 (MQ=255)
                                                                                tAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGAc                <  1:213536/100‑1 (MQ=255)
                                                                                 aaaTCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACa               <  1:2925862/100‑1 (MQ=255)
                                                                                 aaaTCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACa               <  1:1709935/100‑1 (MQ=255)
                                                                                  aaTCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACAg              >  2:1841504/1‑100 (MQ=255)
                                                                                  aaTCAGACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACAg              <  2:592139/100‑1 (MQ=255)
                                                                                  aaTCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGGCAg              <  2:366396/100‑1 (MQ=255)
                                                                                  aaTCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGCATCGGTGGGTGTCACAGGTTGCGTCTGCTCCTGAATGACAg              >  2:2875611/1‑100 (MQ=255)
                                                                                    tCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACAGGc            >  1:3138123/1‑100 (MQ=255)
                                                                                       gACATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACAGGCttt         <  2:1715311/100‑1 (MQ=255)
                                                                                        aCATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACAGGCTTTa        <  2:2901501/100‑1 (MQ=255)
                                                                                         cATGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACAGGCTTTAt       >  1:2534435/1‑100 (MQ=255)
                                                                                           tGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACAGGCTTTATGt     >  1:1510444/1‑100 (MQ=255)
                                                                                           tGGGTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACAGGCTTTATGt     <  1:998009/100‑1 (MQ=255)
                                                                                           tGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACCGGCTTTAtat     >  2:142626/1‑98 (MQ=255)
                                                                                           tGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACAGGCTTTATGt     >  2:1564166/1‑100 (MQ=255)
                                                                                             gATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAAAAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACAGGCTTTATGtca   >  1:1851024/1‑100 (MQ=255)
                                                                                              aTCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACAGGCTTTATGtcat  <  2:1851036/100‑1 (MQ=255)
                                                                                              |                                                                                                   
GTTACAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACAGGCTTTATGTCAT  >  NC_000913/2452245‑2452438

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: