Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,454,402 C→T 5.5% pseudogene (600/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,4020CT 5.5% 226.3 / 2.3 76pseudogene (600/2646 nt)yfcUpseudogene, FimD fimbrial export usher family;putative membrane; Not classified; putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (37/32);  new base T (2/2);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.17e-01

ACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAG  >  NC_000913/2454254‑2454548
                                                                                                                                                    |                                                                                                                                                  
aCGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGAcc                                                                                                                                                    >  1:143624/1‑149 (MQ=255)
    ccgccAGGCATAATAGCGACTCCACTCCCATTTATTTTGTGTGTCATCGCCGCCAAATTCATCGTCTTCATCATTACTGCTCGTATGTTTATAGTTGGTCTGCCAGACAGCGCGCTGACGCCATGGCCCCAGGTTAACCCAGACCGTcc                                                                                                                                                >  5:107414/1‑149 (MQ=255)
                  aGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGAtc                                                                                                                                  >  2:119701/1‑149 (MQ=255)
                    cGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGCCCGTCCCGTTGCCGCTGAtctc                                                                                                                                >  7:146601/1‑149 (MQ=255)
                          cactccCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCAACATCATTACTCCGCGTATGGTAATAGGTAGTCTGCCAGTCAGCGGGCATACGCCAAGGCCACAGTTTAACACAGACCGTCACGTTGCGGCGGATCTCGTTa                                                                                                                            >  6:226072/1‑147 (MQ=255)
                          cacccacATTTTTATTGTGTGTCATAGCCGCCAAATGCATCGTCATCATCATTACTGAGCGTATGTTGATAGGTGGTCTGCCAGTCAGCTCGCATACTCCAGGGCCCCAGGTTCACCCCGACCTTCCCGTTGCCGCTGATCTCGTTa                                                                                                                            <  6:387313/141‑1 (MQ=255)
                            ctccCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGt                                                                                                                        <  2:107631/149‑1 (MQ=255)
                             tccCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCTTCATTAGAGGTCGTATGTTTATAGTTGGTCTACCAGTCAGCTGGCATACGCCATGGACAAAGGTTAACCCCGCCCGTCACGTAGCCGCTGATCTCGTTACTGt                                                                                                                        >  8:246772/1‑148 (MQ=255)
                             tccCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGt                                                                                                                        >  3:87806/1‑148 (MQ=255)
                               ccATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACTGTCCCGTTGCCGCTGATCTCGTTACTGTCa                                                                                                                      >  1:279370/1‑148 (MQ=255)
                               ccATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACTGTCCCGTTGCCGCTGATCTCGTTACTGTCAt                                                                                                                     >  7:69377/1‑149 (MQ=255)
                               ccATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTg                                                                                                                                           <  1:291134/127‑1 (MQ=255)
                               ccATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGATAACCCCGACCGTCCCGTTg                                                                                                                                           >  2:291134/1‑127 (MQ=255)
                                cATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATc                                                                                                                    <  4:140747/149‑1 (MQ=255)
                                                    gccgccAAATTCATCGTCATCATAATTACTGCTCGTATGTTGATAGTTGGTCTGCCAGTCAGCTCGCATACGCCATGGCCACAGGTTAACCCCGACCGTCCCTTTTCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTtcgtgtc                                                                                                <  7:100704/149‑1 (MQ=255)
                                                      cgccAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTtcgtgtcgt                                                                                              <  2:143624/149‑1 (MQ=255)
                                                        ccAAATTCATCGTCATCATCATTACTGCCAGTATGTTGATAGTTTGACTGCCAGTCAGCTCGCATACGCCATGGCCACAGGTTACACCCGACCGTACCGATGCCGCTGATATCGTTACTGTCATCGCCGCCATTTTCTTcgtgtcgtgt                                                                                            >  1:373993/1‑149 (MQ=255)
                                                        ccAAAGACATCGTCATCATCATTACTGCTCGTATGTTGAAAGTTGGGCTGCCAGGCTGCGCGCATACGCCATGGCCCCAGGTTACCCCCGACCGTCCCGTAGCCGCTGATCTAGTTAGTGTCATCGCCGCCATTTTCTTcgtgtcgtgt                                                                                            <  4:65121/149‑1 (MQ=255)
                                                          aaaTTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGACAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTTATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGttt                                                                                          >  3:96617/1‑149 (MQ=255)
                                                            aTTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTgc                                                                                        <  3:44426/149‑1 (MQ=255)
                                                             ttCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTgcg                                                                                       >  4:282530/1‑149 (MQ=255)
                                                                aTCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACTGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCa                                                                                     <  8:69377/148‑1 (MQ=255)
                                                                aTCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACTGTCCCGTTGCCGCTGATCTCGTTACTGTCATAGCCGCCATTTTCTTCGTGTCGTGTTTGCGCa                                                                                     <  2:279370/148‑1 (MQ=255)
                                                                   gTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCTCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTTATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGa                                                                                 <  2:219054/149‑1 (MQ=255)
                                                                    tcatcatcaTTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGAt                                                                                <  7:46021/149‑1 (MQ=255)
                                                                       tcatcatTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGACCCAGGTTAACACCGACCGTCACGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCt                                                                             >  1:279500/1‑149 (MQ=255)
                                                                          tcatTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTTATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTa                                                                          <  7:246772/149‑1 (MQ=255)
                                                                           catTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAg                                                                         >  7:375378/1‑149 (MQ=255)
                                                                                   cgcgTATGTTGCTAGTAGGGCTGCCAGTCAGCGGGCATACGCCATGGCCCCAGGGTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCa                                                                   <  5:226072/147‑1 (MQ=255)
                                                                                   cgcgTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGt                                                                                   <  5:270894/131‑1 (MQ=255)
                                                                                   cgcgTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGt                                                                                   >  6:270894/1‑131 (MQ=255)
                                                                                   cgcgTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCa                                                                   <  8:375378/147‑1 (MQ=255)
                                                                                   cgcgTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAAt                                                                 <  1:119701/149‑1 (MQ=255)
                                                                                      gTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGAt                                                              >  8:84761/1‑149 (MQ=255)
                                                                                          gTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCg                                                          >  1:389176/1‑149 (MQ=255)
                                                                                           ttGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTg                                                                                  <  1:327430/124‑1 (MQ=255)
                                                                                           ttGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTg                                                                                  >  2:327430/1‑124 (MQ=255)
                                                                                            tGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCg                                                          >  2:88510/1‑147 (MQ=255)
                                                                                              aTAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGCTTAACCCCGACCGTCCCGTTGCCGCTGATCTCCTTACTGTCATAGCCGCCATATTCTTCGTGTCGTGGTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCg                                                          >  6:268695/1‑145 (MQ=255)
                                                                                                 gTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCaa                                                                  >  4:168292/1‑134 (MQ=255)
                                                                                                 gTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCaa                                                                  <  3:168292/134‑1 (MQ=255)
                                                                                                  ttGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGCTCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGc                                                   >  5:182613/1‑148 (MQ=255)
                                                                                                     gTCTTCCAGTCCGCGAGCACACTCCATGGCACCCGGATACCCCAGACCGTACCGATGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCcatc                                               <  3:206101/149‑1 (MQ=255)
                                                                                                     gTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCGCGTCACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCcat                                                >  1:73741/1‑148 (MQ=255)
                                                                                                      tCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCcatc                                               >  7:16743/1‑148 (MQ=255)
                                                                                                        tGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCCatcatc                                            >  4:80648/1‑149 (MQ=255)
                                                                                                        tGCCAGTCAGCGCGCATACGCCATGGCCCCAGCTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCCatcatc                                            >  5:336189/1‑149 (MQ=255)
                                                                                                         gCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGCCCGTCCCGTTGCCGCTGATCTCGATACTGTCATCGCCCCCATTTTCTTCGTCTCGTGTATGCGCAGTGATGCTGTAGTCCCCAATGATCCCGGACATGCCATCATcc                                           >  4:212247/1‑149 (MQ=255)
                                                                                                         gCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCCATCATcc                                           >  8:288691/1‑149 (MQ=255)
                                                                                                          ccAGTCAGCGCGCATACGCCATGGCCCCAGGTTCACCCCGACCGGCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCCatcatc                                            <  8:289732/147‑1 (MQ=255)
                                                                                                          ccAGTAAGCGCGCATACGACATGGCCCAATGTTAACCCAGACCGTCCCGTTTCCTCTGATCTCGTTACTGTCATCGCCGCAATTTTCTTCGTGTCGTGTTTTCGAAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCCATCATcc                                           >  3:260003/1‑148 (MQ=255)
                                                                                                           cAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCCATCATccc                                          >  3:190849/1‑148 (MQ=255)
                                                                                                            aGTCAGCGCGCTTACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGGTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCCATCATCCCaa                                        <  8:146601/149‑1 (MQ=255)
                                                                                                             gTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCATTGATCCCGGAGATGCCATCATCCCAAc                                       <  4:96617/149‑1 (MQ=255)
                                                                                                             gTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCcatca                                              >  1:357180/1‑142 (MQ=255)
                                                                                                             gTCAGCGAGCATACGCTATGGCCCCAGGTTAACCCCGACCGTCGCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCcatca                                              <  2:357180/142‑1 (MQ=255)
                                                                                                                 gcgcgcATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCCATCATcc                                           >  4:53414/1‑141 (MQ=255)
                                                                                                                 gcgcgcATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCCATCATcc                                           <  3:53414/141‑1 (MQ=255)
                                                                                                                 gcgcgcATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCCATCATCCCAACgaga                                   <  7:84761/149‑1 (MQ=255)
                                                                                                                   gcgcATACGCCTTTGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCCATCATCCCAACGAGAAg                                 <  6:336189/149‑1 (MQ=255)
                                                                                                                     gcATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCCATCATCCCCACGAGAAGGt                               >  7:257770/1‑149 (MQ=255)
                                                                                                                     gcATACGCCATGGCCCCAGGGTAACCCCTACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCCATCATCCCAACGAGAAGGt                               <  4:276624/149‑1 (MQ=255)
                                                                                                                          cGCCATTTCCCCAGGTTAACCCCGACAGTCCCGTATCATCGGATCTCGTTACTGTCATCGCCGCAATTTGATTCTGGTCGTGTTTGCGATTTGATGCTGTATTCCGCAATGCTCCCGGAGATGCCATCATCCCAACGAGCAGGGGGag                           <  7:274653/148‑2 (MQ=255)
                                                                                                                          cGCCATGGCCCCAGGGTAACCCAGACCGTCCCGTTGCCGCTTATCTCGTTAGTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCCATCATCCCAACGAGAAGGTGGATc                          <  8:213990/149‑1 (MQ=255)
                                                                                                                           gCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCCATCATCCCAACGAGAAGGTGGATcc                         >  5:71477/1‑149 (MQ=255)
                                                                                                                             cATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTAATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCGGAGCTGCCATCATCCCAACGAGAAGGTGGATcc                         >  2:393940/1‑147 (MQ=255)
                                                                                                                                  ccccCGGTTATCCCCCACGGTCCCGTTGCCGCTGATGTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCCATCATCCCAACGAGAAGGTGGATCCCAGTCgg                  <  6:258921/149‑1 (MQ=255)
                                                                                                                                   cccAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCGATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCCATCATCCCAACGAGAAGGTGGATCCCAGTCggg                 >  7:385450/1‑149 (MQ=255)
                                                                                                                                   cccAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCCATCATCCCAACgaga                                   >  6:44400/1‑131 (MQ=255)
                                                                                                                                   cccAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCCATCATCCCAACgaga                                   <  5:44400/131‑1 (MQ=255)
                                                                                                                                   cccAGCTTAGACCCGACCGTCCCCTTGCCGCTGATCTCTTTACTGTCATCGCCGACATGTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTTGTCCGCAATGATCCCGGAGATGCCATCATCGCAACGAGAAGGTGGATCCCTGTCggg                 <  6:107414/149‑1 (MQ=255)
                                                                                                                                       ggTTAACTCCGACGGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCATCGTGGCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCCATCATCCCAACGATAAGGTGGATCCCAGTCGGGCCAg             <  2:373993/149‑1 (MQ=255)
                                                                                                                                               ccGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCCATCATCCCAACGAGAAGGTGGAt                           <  4:83176/127‑1 (MQ=255)
                                                                                                                                               ccGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCCATCATCCCAACGAGAAGGTGGAt                           >  3:83176/1‑127 (MQ=255)
                                                                                                                                                cGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCCATCATCCCAACGAGAAGGTGGATACAAGTCGGGCCAGGGAGATTc     >  2:247055/1‑148 (MQ=255)
                                                                                                                                                cGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGAAATGATCCCGGAGATGCCATCATCCCTACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCg    <  2:279500/149‑1 (MQ=255)
                                                                                                                                                  aCCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAg  <  3:282530/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGCTGATCTCGTTACTGTCATCGCCGCCATTTTCTTCGTGTCGTGTTTGCGCAGTGATGCTGTAGTCCGCAATGATCCCGGAGATGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAG  >  NC_000913/2454254‑2454548

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