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

GCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTC  >  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)
                        gCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCAtttt                                                                                                                    >  4:241497/1‑149 (MQ=255)
                        gCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCAtttt                                                                                                                    <  1:156275/149‑1 (MQ=255)
                        gCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCAtttt                                                                                                                    >  8:369465/1‑149 (MQ=255)
                         cAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCAttttt                                                                                                                   >  2:433196/1‑149 (MQ=255)
                          aaCATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCAtttttt                                                                                                                  >  4:361596/1‑149 (MQ=255)
                            cATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCAtttttttt                                                                                                                <  7:410684/149‑1 (MQ=255)
                              tttGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGactccactcc                                                                                                                           >  4:230262/1‑137 (MQ=255)
                              tttGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGactccactcc                                                                                                                           <  3:230262/137‑1 (MQ=255)
                                tGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCAGTTTTTTtgtgt                                                                                                            >  6:238821/1‑149 (MQ=255)
                                tGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCAGTTTTTTtgtgt                                                                                                            >  1:241046/1‑149 (MQ=255)
                                 gATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTtgtgtg                                                                                                           >  6:81507/1‑149 (MQ=255)
                                  aTCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCAtttttttt                                                                                                                >  6:72748/1‑143 (MQ=255)
                                  aTCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCAtttttttt                                                                                                                <  5:72748/143‑1 (MQ=255)
                                     gTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTtgt                                                                                                              >  1:167903/1‑142 (MQ=255)
                                     gTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTtgt                                                                                                              <  2:167903/142‑1 (MQ=255)
                                      tCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGc                                                                                                                                      >  7:7035/1‑118 (MQ=255)
                                      tCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGc                                                                                                                                      <  8:7035/118‑1 (MQ=255)
                                      tCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATc                                                                                                      >  7:102011/1‑149 (MQ=255)
                                      tCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATc                                                                                                      >  4:230301/1‑149 (MQ=255)
                                      tCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAAGAGCGACTCCACTCCCATTTTTTTtgt                                                                                                              >  1:60445/1‑141 (MQ=255)
                                      tCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAAAAGCGACTCCACTCCCATTTTTTTtgt                                                                                                              <  2:60445/141‑1 (MQ=255)
                                            cTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATcgccgc                                                                                                 >  3:141620/1‑148 (MQ=255)
                                               acacTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCaaa                                                                                             <  5:81507/149‑1 (MQ=255)
                                                 acTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCAttttttt                                                                                                                 <  8:2351/127‑1 (MQ=255)
                                                 acTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCAttttttt                                                                                                                 >  7:2351/1‑127 (MQ=255)
                                                     ccaccaACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCAt                                                                                        <  3:361596/148‑1 (MQ=255)
                                                     ccaccaACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATc                                                                                       <  3:241497/149‑1 (MQ=255)
                                                           aCATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGtcat                                                                                  >  7:93803/1‑148 (MQ=255)
                                                               aaTTAAAACCGTCGAAAATATCGGAATTGAGGTATTCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGtcatcatca                                                                             >  7:377359/1‑149 (MQ=255)
                                                               aaTTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGtcatcatca                                                                             >  1:74800/1‑149 (MQ=255)
                                                               aaTTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTATAATGACGGTAACGCCCGCCAGGCATAATAGCGACTACACTCCCATTATTTATGTGTGACATCGCCGCCAAATTCATCGtcataatca                                                                             >  6:327658/1‑149 (MQ=255)
                                                                 ttAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTcatcatcat                                                                            >  8:88076/1‑148 (MQ=255)
                                                                    aaaCCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACt                                                                        <  1:70755/149‑1 (MQ=255)
                                                                     aaCCGTCGAAATTATCGGACTTGAGGTAAGCTTCGCCCAGCGCCCGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATATTAGCGACTCCACTCCCATTTTTTTTGTGGGTCATCGCCGCCAAATTCATCGTCATCATCATTACTg                                                                       <  6:203820/149‑1 (MQ=255)
                                                                          tCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGcgcg                                                                   >  6:322149/1‑148 (MQ=255)
                                                                            gAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTAt                                                                >  8:323007/1‑149 (MQ=255)
                                                                              aaaTATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTAt                                                                <  5:327658/147‑1 (MQ=255)
                                                                               aaTATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCAttttt                                                                                                                   >  3:21006/1‑95 (MQ=255)
                                                                               aaTATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCAttttt                                                                                                                   <  4:21006/95‑1 (MQ=255)
                                                                                 tatCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGa                                                           <  3:131858/149‑1 (MQ=255)
                                                                                 tatCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCACACCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGa                                                           <  1:285773/149‑1 (MQ=255)
                                                                                     ggAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAg                                                        <  4:412920/148‑1 (MQ=255)
                                                                                     ggAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGt                                                       <  8:102011/149‑1 (MQ=255)
                                                                                         ttGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGt                                                   <  7:323007/149‑1 (MQ=255)
                                                                                            aGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCAt                                                                                                       >  8:58578/1‑94 (MQ=255)
                                                                                            aGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCAt                                                                                                       <  7:58578/94‑1 (MQ=255)
                                                                                            aGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCt                                                 <  2:74800/148‑1 (MQ=255)
                                                                                                aaTCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTAt                                                                >  7:212082/1‑130 (MQ=255)
                                                                                                aaTCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAg                                            >  8:192202/1‑149 (MQ=255)
                                                                                                aaTCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCACCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAg                                            >  7:251363/1‑149 (MQ=255)
                                                                                                aaTCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCACCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAg                                            >  7:258045/1‑149 (MQ=255)
                                                                                                aaTCTTCGCCCAGCGCCAGTTTGGCTTTTAATGAACGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTAt                                                                <  8:212082/130‑1 (MQ=255)
                                                                                                 aTCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCCTCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAg                                            <  5:238821/148‑1 (MQ=255)
                                                                                                 aTCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGt                                           <  8:93803/149‑1 (MQ=255)
                                                                                                      cGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAgc                                       <  8:388125/148‑1 (MQ=255)
                                                                                                      cGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAgc                                       <  3:230301/148‑1 (MQ=255)
                                                                                                      cGCCCAGCGCCAGTTTGGATTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAgc                                       <  8:377359/148‑1 (MQ=255)
                                                                                                       gCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAgcgc                                     <  1:433196/149‑1 (MQ=255)
                                                                                                        cccAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGt                                           >  1:11802/1‑142 (MQ=255)
                                                                                                        cccAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGt                                           <  2:11802/142‑1 (MQ=255)
                                                                                                         ccAGCGCCAGTTTGTCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGa                                                           <  2:76256/125‑1 (MQ=255)
                                                                                                         ccAGCGCCAGTTTGCCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGa                                                           >  1:76256/1‑125 (MQ=255)
                                                                                                              gcCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACg                              >  3:326633/1‑149 (MQ=255)
                                                                                                                 aGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTAt                                                                <  7:27351/112‑1 (MQ=255)
                                                                                                                 aGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTAt                                                                >  8:27351/1‑112 (MQ=255)
                                                                                                                  gTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCAt                          <  7:369465/149‑1 (MQ=255)
                                                                                                                  gTTTGGCTTTTAATGACGGTAACGCCCACCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCAt                          <  8:251363/149‑1 (MQ=255)
                                                                                                                  gTTTGGCTTTTAATGACGGTAACGCCCACCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCAt                          <  8:258045/149‑1 (MQ=255)
                                                                                                                     tGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCGATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATgg                        <  2:106027/148‑1 (MQ=255)
                                                                                                                      ggCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGc                       <  4:141620/148‑1 (MQ=255)
                                                                                                                        cTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTAGATAGTTGGTCTGCCAGTAAGCGCGCATCAGCAATGGccc                     >  2:162683/1‑148 (MQ=255)
                                                                                                                         ttttAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGccc                     >  7:42824/1‑147 (MQ=255)
                                                                                                                         ttttAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGccc                     >  1:66852/1‑147 (MQ=255)
                                                                                                                         ttttAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGACGCCAAATTCATCGTCATCATCATTACTGCGAGTATGTTGATAGTTGGTCTGCCAGTCCGCGCGCATACGCCATGGcac                     >  3:39102/1‑145 (MQ=255)
                                                                                                                         ttttAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGCGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGccc                     >  8:212849/1‑147 (MQ=255)
                                                                                                                         ttttAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGCGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGccc                     >  4:158406/1‑147 (MQ=255)
                                                                                                                             aaTGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAATTTCATCGACATCAGAATTCCTGCGAGTAGGTTAATAGTTGGTCTGACAGTCAGCGCGAATGCGCCATGGCCACACGt                >  6:67797/1‑148 (MQ=255)
                                                                                                                              aTGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGtt               <  7:192202/148‑1 (MQ=255)
                                                                                                                               tGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTaa             >  3:349314/1‑149 (MQ=255)
                                                                                                                               tGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTaa             >  1:19370/1‑149 (MQ=255)
                                                                                                                                  cGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAAcc           <  4:326633/148‑1 (MQ=255)
                                                                                                                                          cccgccAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGAATACGCCATGGCCACAGGTTAACCACGACCGTc  >  1:84352/1‑149 (MQ=255)
                                                                                                                                           ccgccAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGa                                                           <  5:249077/91‑1 (MQ=255)
                                                                                                                                           ccgccAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGa                                                           >  6:249077/1‑91 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
GCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTC  >  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: