Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,454,260 G→A 5.5% pseudogene (742/2646 nt) yfcU ← 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,454,2600GA 5.5% 283.0 / 4.4 90pseudogene (742/2646 nt)yfcUpseudogene, FimD fimbrial export usher family;putative membrane; Not classified; putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (45/40);  new base A (2/3);  total (47/43)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.02e-01

GCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTC  >  NC_000913/2454119‑2454405
                                                                                                                                             |                                                                                                                                                  
gcgcGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGc                                                                                                                                              <  8:53366/148‑1 (MQ=255)
          gcgcAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATATTTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGa                                                                                                                                    >  1:324117/1‑148 (MQ=255)
          gcgcAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGa                                                                                                                                    <  2:324117/148‑1 (MQ=255)
          gcgcAGATTGCGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACAAAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGac                                                                                                                                   >  1:106027/1‑149 (MQ=255)
            gcAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGATTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGactc                                                                                                                                 >  4:374281/1‑149 (MQ=255)
                aTTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGc                                                                                                                                              <  7:410727/132‑1 (MQ=255)
                aTTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGc                                                                                                                                              >  8:410727/1‑132 (MQ=255)
                 ttGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGactccactc                                                                                                                            <  6:209693/149‑1 (MQ=255)
                  tGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGactccactcc                                                                                                                           >  3:412920/1‑149 (MQ=255)
                        gCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑tttt                                                                                                                    >  4:241497/1‑149 (MQ=255)
                        gCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑tttt                                                                                                                    <  1:156275/149‑1 (MQ=255)
                        gCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑tttt                                                                                                                    >  8:369465/1‑149 (MQ=255)
                         cAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑ttttt                                                                                                                   >  2:433196/1‑149 (MQ=255)
                          aaCATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑tttttt                                                                                                                  >  4:361596/1‑149 (MQ=255)
                            cATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑tttttttt                                                                                                                <  7:410684/149‑1 (MQ=255)
                              tttGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGactccactcc                                                                                                                           >  4:230262/1‑137 (MQ=255)
                              tttGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGactccactcc                                                                                                                           <  3:230262/137‑1 (MQ=255)
                                tGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑GTTTTTTtgtgt                                                                                                            >  6:238821/1‑149 (MQ=255)
                                tGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑GTTTTTTtgtgt                                                                                                            >  1:241046/1‑149 (MQ=255)
                                 gATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTtgtgtg                                                                                                           >  6:81507/1‑149 (MQ=255)
                                  aTCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑tttttttt                                                                                                                >  6:72748/1‑143 (MQ=255)
                                  aTCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑tttttttt                                                                                                                <  5:72748/143‑1 (MQ=255)
                                     gTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTtgt                                                                                                              >  1:167903/1‑142 (MQ=255)
                                     gTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTtgt                                                                                                              <  2:167903/142‑1 (MQ=255)
                                      tCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGc                                                                                                                                      >  7:7035/1‑118 (MQ=255)
                                      tCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGc                                                                                                                                      <  8:7035/118‑1 (MQ=255)
                                      tCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATc                                                                                                      >  7:102011/1‑149 (MQ=255)
                                      tCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATc                                                                                                      >  4:230301/1‑149 (MQ=255)
                                      tCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAAGAGCGACTCCACTCCCA‑TTTTTTTtgt                                                                                                              >  1:60445/1‑141 (MQ=255)
                                      tCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAAAAGCGACTCCACTCCCA‑TTTTTTTtgt                                                                                                              <  2:60445/141‑1 (MQ=255)
                                            cTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATcgccgc                                                                                                 >  3:141620/1‑148 (MQ=255)
                                               acacTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCaaa                                                                                             <  5:81507/149‑1 (MQ=255)
                                                 acTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑ttttttt                                                                                                                 <  8:2351/127‑1 (MQ=255)
                                                 acTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑ttttttt                                                                                                                 >  7:2351/1‑127 (MQ=255)
                                                     ccaccaACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCAt                                                                                        <  3:361596/148‑1 (MQ=255)
                                                     ccaccaACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATc                                                                                       <  3:241497/149‑1 (MQ=255)
                                                           aCATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGtcat                                                                                  >  7:93803/1‑148 (MQ=255)
                                                               aaTTAAAACCGTCGAAAATATCGGAATTGAGGTATTCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGtcatcatca                                                                             >  7:377359/1‑149 (MQ=255)
                                                               aaTTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGtcatcatca                                                                             >  1:74800/1‑149 (MQ=255)
                                                               aaTTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTATAATGACGGTAACGCCCGCCAGGCATAATAGCGACTACACTCCCA‑TTATTTATGTGTGACATCGCCGCCAAATTCATCGtcataatca                                                                             >  6:327658/1‑149 (MQ=255)
                                                                 ttAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTcatcatcat                                                                            >  8:88076/1‑148 (MQ=255)
                                                                    aaaCCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACt                                                                        <  1:70755/149‑1 (MQ=255)
                                                                     aaCCGTCGAAATTATCGGACTTGAGGTAAGCTTCGCCCAGCGCCCGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATATTAGCGACTCCACTCCCA‑TTTTTTTTGTGGGTCATCGCCGCCAAATTCATCGTCATCATCATTACTg                                                                       <  6:203820/149‑1 (MQ=255)
                                                                          tCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGcgcg                                                                   >  6:322149/1‑148 (MQ=255)
                                                                            gAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTAt                                                                >  8:323007/1‑149 (MQ=255)
                                                                              aaaTATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTAt                                                                <  5:327658/147‑1 (MQ=255)
                                                                               aaTATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑ttttt                                                                                                                   >  3:21006/1‑95 (MQ=255)
                                                                               aaTATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑ttttt                                                                                                                   <  4:21006/95‑1 (MQ=255)
                                                                                 tatCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGa                                                           <  3:131858/149‑1 (MQ=255)
                                                                                 tatCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCACACCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGa                                                           <  1:285773/149‑1 (MQ=255)
                                                                                     ggAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAg                                                        <  4:412920/148‑1 (MQ=255)
                                                                                     ggAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGt                                                       <  8:102011/149‑1 (MQ=255)
                                                                                         ttGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGt                                                   <  7:323007/149‑1 (MQ=255)
                                                                                            aGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCAt                                                                                                       >  8:58578/1‑94 (MQ=255)
                                                                                            aGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCAt                                                                                                       <  7:58578/94‑1 (MQ=255)
                                                                                            aGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCt                                                 <  2:74800/148‑1 (MQ=255)
                                                                                                aaTCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTAt                                                                >  7:212082/1‑130 (MQ=255)
                                                                                                aaTCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAg                                            >  8:192202/1‑149 (MQ=255)
                                                                                                aaTCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCACCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAg                                            >  7:251363/1‑149 (MQ=255)
                                                                                                aaTCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCACCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAg                                            >  7:258045/1‑149 (MQ=255)
                                                                                                aaTCTTCGCCCAGCGCCAGTTTGGCTTTTAATGAACGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTAt                                                                <  8:212082/130‑1 (MQ=255)
                                                                                                 aTCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCCTCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAg                                            <  5:238821/148‑1 (MQ=255)
                                                                                                 aTCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGt                                           <  8:93803/149‑1 (MQ=255)
                                                                                                      cGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAgc                                       <  8:388125/148‑1 (MQ=255)
                                                                                                      cGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAgc                                       <  3:230301/148‑1 (MQ=255)
                                                                                                      cGCCCAGCGCCAGTTTGGATTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAgc                                       <  8:377359/148‑1 (MQ=255)
                                                                                                       gCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAgcgc                                     <  1:433196/149‑1 (MQ=255)
                                                                                                        cccAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGt                                           >  1:11802/1‑142 (MQ=255)
                                                                                                        cccAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGt                                           <  2:11802/142‑1 (MQ=255)
                                                                                                         ccAGCGCCAGTTTGTCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGa                                                           <  2:76256/125‑1 (MQ=255)
                                                                                                         ccAGCGCCAGTTTGCCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGa                                                           >  1:76256/1‑125 (MQ=255)
                                                                                                              gcCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACg                              >  3:326633/1‑149 (MQ=255)
                                                                                                                 aGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTAt                                                                <  7:27351/112‑1 (MQ=255)
                                                                                                                 aGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTAt                                                                >  8:27351/1‑112 (MQ=255)
                                                                                                                  gTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCAt                          <  7:369465/149‑1 (MQ=255)
                                                                                                                  gTTTGGCTTTTAATGACGGTAACGCCCACCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCAt                          <  8:251363/149‑1 (MQ=255)
                                                                                                                  gTTTGGCTTTTAATGACGGTAACGCCCACCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCAt                          <  8:258045/149‑1 (MQ=255)
                                                                                                                     tGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCGA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATgg                        <  2:106027/148‑1 (MQ=255)
                                                                                                                      ggCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGc                       <  4:141620/148‑1 (MQ=255)
                                                                                                                        cTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTAGATAGTTGGTCTGCCAGTAAGCGCGCATCAGCAATGGccc                     >  2:162683/1‑148 (MQ=255)
                                                                                                                         ttttAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGccc                     >  7:42824/1‑147 (MQ=255)
                                                                                                                         ttttAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGccc                     >  1:66852/1‑147 (MQ=255)
                                                                                                                         ttttAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGACGCCAAATTCATCGTCATCATCATTACTGCGAGTATGTTGATAGTTGGTCTGCCAGTCCGCGCGCATACGCCATGGcac                     >  3:39102/1‑145 (MQ=255)
                                                                                                                         ttttAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGCGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGccc                     >  8:212849/1‑147 (MQ=255)
                                                                                                                         ttttAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGCGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGccc                     >  4:158406/1‑147 (MQ=255)
                                                                                                                             aaTGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAATTTCATCGACATCAGAATTCCTGCGAGTAGGTTAATAGTTGGTCTGACAGTCAGCGCGAATGCGCCATGGCCACACGt                >  6:67797/1‑148 (MQ=255)
                                                                                                                              aTGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGtt               <  7:192202/148‑1 (MQ=255)
                                                                                                                               tGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTaa             >  3:349314/1‑149 (MQ=255)
                                                                                                                               tGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTaa             >  1:19370/1‑149 (MQ=255)
                                                                                                                                  cGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAAcc           <  4:326633/148‑1 (MQ=255)
                                                                                                                                          cccgccAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGAATACGCCATGGCCACAGGTTAACCACGACCGTc  >  1:84352/1‑149 (MQ=255)
                                                                                                                                           ccgccAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGa                                                           <  5:249077/91‑1 (MQ=255)
                                                                                                                                           ccgccAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGa                                                           >  6:249077/1‑91 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
GCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCA‑TTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTC  >  NC_000913/2454119‑2454405

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