Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,447,536 C→A 7.4% pseudogene (945/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
*CP0092732,447,5360CA 7.4% 976.1 / 35.8 265pseudogene (945/2646 nt)yfcUpseudogene, FimD fimbrial export usher family;putative membrane; Not classified; putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (106/137);  new base A (13/9);  total (119/146)
Fisher's exact test for biased strand distribution p-value = 1.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 269 total aligned reads displayed.

CGAGTGAGGTATGGCATCGAGGCGGTGCTGATGTCATATTCCTGCACCTGGCCGTTCTGTTCTTCAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGG  >  CP009273/2447387‑2447684
                                                                                                                                                     |                                                                                                                                                    
cGAGTGAGGTATGGCATCGAGGCGGTGCTGATGTCATATTCCTGCACCTGGCCGTTCTGTTCTTCAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATc                                                                                                                                                      <  1:4449076/150‑1 (MQ=255)
 gAGTGAGGTATGGCATCGAGGCGGTGCTGATGTCATATTCCTGCACCTGGCCGTTCTGTTCTTCAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATca                                                                                                                                                     <  1:1949611/150‑1 (MQ=255)
     gAGGTATGGCATCGAGGCGGTGCTGATGTCATATTCCTGCACCTGGCCGTTCTGTTCTTAAATGCGAATATTCAACGTTCCGGAGACGGAATCACCAAGATCCAGAAAACGAAACGGTCCAGCCGGAACCAGCGTTTCGTAAATAACACg                                                                                                                                                 >  1:4385277/1‑150 (MQ=255)
     gAGGTATGGCATCGAGGCGGTGCTGATGTCATATTCCTGCACCTGGCCGTGCTGTTGTTGAATGCGAATATTCAACGGACCGGAGACGGAAAAAACAACATCCAGAATACGAAACCGTCCCGCCCGAAACAGCGTTTCGTAAATCACACg                                                                                                                                                 >  1:1197910/1‑150 (MQ=255)
      aGGTATGGCATCGAGGCGGTGCTGATGTCATATTCCTGCACCTGGCCGTTCTGTTCTTCAATGCGAATATTCCACGTACCGGAGACGGAATCACCAAGATCCAGAAAACGAAAAGGTCCAGCCGGAACCAGCGTTTAGTAAAAAacaccc                                                                                                                                                >  1:1822494/1‑148 (MQ=255)
      aGGTATGGCATCGAGGCGGTGCTGATGTCATATTCCTGCACCTGGCCGTTCTGTTCTTCAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGc                                                                                                                                                >  2:4542408/1‑150 (MQ=255)
         tATGGCATCGAGGCGGTGCTGATGTCATATTCCTGCACCTGGCCGTTGTGTTCTTAAATTCGAATAAGCAAAGTACCCGAGACGGAATCACCAAGATCCAGAAAAAGAAAAGGTCCAGCCGGAACCAGCGTTTCGTAAATAAAAcccccc                                                                                                                                             >  1:3639123/1‑150 (MQ=255)
         tATGGCATCGAGGCGGTGCTGATGTCATATTCCTGCACCTGGCCGTTCTGTTCTTCAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGcccc                                                                                                                                             <  2:3018970/150‑1 (MQ=255)
           tGGCATCGAGGCGGTGCTGATGTCATATTCCTGCACCTGGCCGTTCTGTTCTTCAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCAt                                                                                                                                           <  2:1647978/150‑1 (MQ=255)
           tGGCATCGAGGCGGTGCTGATGTCATATTCCTGCACCTGGCCGTTCTGTTCTTCAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCAt                                                                                                                                           <  2:3223489/150‑1 (MQ=255)
           tGGCATCGAGGCGGTGCTGATGTCATATTCCTGCACCTGGCCGTTCTGTTCTTCAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCAt                                                                                                                                           <  2:4628078/150‑1 (MQ=255)
             gCATCGAGGCGGTGCTGATGTCATATTCCTGCACCTGGCCGTTCTGTTCTTCAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCt                                                                                                                                         <  2:1897947/150‑1 (MQ=255)
             gCATCGAGGCGGTGCTGATGTCATATTCCTGCACCTGGCCGTTCTGTTCTTCAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCt                                                                                                                                         <  2:3726864/150‑1 (MQ=255)
              cATCGAGGCGGTGCTGATGTCATATTCCTGCACCTGGCCGTTCTGTTCTTCAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAAACACACGCCCCATCTg                                                                                                                                        >  2:4436039/1‑150 (MQ=255)
                  gAGGCGGTGCTGATGTCATATTCCTGCACCTGGCCGTGCTGTTCTTCAATGCGCATATTCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCCGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTg                                                                                                                                    >  1:587501/1‑150 (MQ=255)
                      cGGTGCTGATGTCATATTCCTGCACCTGGCCGTTGTGTTGTTCAATTCGAATATTCAACGTACCGGAGACGGAATCACCAAGATCCTGAATAAGAAAAGGTCCAGCCGGAACCAGCGTTTCGTAAATCACACGCCCCCACTGGCTGACgg                                                                                                                                >  1:977241/1‑150 (MQ=255)
                       ggTGCTGATGTCATATTCCTGCACCTGGCCGTTCTGTTCTTCAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGt                                                                                                                               <  1:747942/150‑1 (MQ=255)
                       ggTGCTGATGTCATAGTCCTGCACCTGGCCGTTCTGTTCTTCAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGt                                                                                                                               <  1:3115830/150‑1 (MQ=255)
                        gTGCTGATGTCATATTCCTGCACCTGGCCGTTCTGTTCTTCAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTc                                                                                                                              <  1:2267919/150‑1 (MQ=255)
                            tGATGTCATATTCCTGCACCTGGCCGTTCTGTTATTCAAAGCGAATATTCAACGTAACGGAGACGGAATCACCAAGATACAGAATACGAAACGGTCCAGCCGGAAACAGCGTTTCGTAAATAACACGCCCCATCTGGCTGACGGTCAcat                                                                                                                          >  1:2292921/1‑148 (MQ=255)
                                 tCATATTCCTGCACCTGGCCGTTCTGTTCTTCAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCt                                                                                                                     <  1:469859/150‑1 (MQ=255)
                                 tCATATTCCTGCACCTGGCCGTTCTGTTCTTCAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCt                                                                                                                     <  1:1431578/150‑1 (MQ=255)
                                    tatTCCTGCACCTGGCCGTTCTGTTCTTCAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGACTACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATAACACGCCCCATCTGGCTGACGGTCACTTTTGCtgtg                                                                                                                  <  2:3396321/150‑1 (MQ=255)
                                      ttCCTGCACCTGGCCGTTCTGTTCTTCAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTgtggt                                                                                                                <  2:1456124/150‑1 (MQ=255)
                                       tCCTGCACCTGGCCGTTCTGTTCTTAAAAGCGAATATTCAACGTAACCGATACGGAATCACCAAGATACATAATAAGAAAAGGTCCAGCCCGAACCTGCGTTTAGTAAATAACACGCCCCCACTGGCTGAAGGTCACTTTTGCTgtggtg                                                                                                               >  1:1451768/1‑150 (MQ=255)
                                       tCCTGCACCTGGCCGTTCTGTTATTCAAAGCGAATATGCAACGTACCGGAGACGGAAACACCAAGATCCTGAAAACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTgtggtg                                                                                                               >  2:575808/1‑150 (MQ=255)
                                        ccTGCACCTGGCCGTTCTGTTCTTCAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGgtgt                                                                                                              <  2:1305865/150‑1 (MQ=255)
                                        ccTGCACCTGGCCGTTCTGTTCTTCAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGgtgt                                                                                                              <  2:2010364/150‑1 (MQ=255)
                                         cTGCACCTGGCCGTTCTGTTCTTCAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGgtgtg                                                                                                             <  2:3867372/150‑1 (MQ=255)
                                          tGCACCTGGCCGTTCTGTTATTCAATGCGAATATGCAACGTAACGGAGACGGAATCACCAAGATCCAGAATACGAAACCGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCCTCTGGCTGACGGTCACTTTTGCTGTGgtgtgt                                                                                                            >  2:2910758/1‑150 (MQ=255)
                                            cACCTGGCCGTTCTGTTCTTCAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGc                                                                                                          <  2:3585443/150‑1 (MQ=255)
                                             aCCTGGCCGTTCTGTTCTTCAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGcc                                                                                                         <  2:2994602/150‑1 (MQ=255)
                                              ccTGGCCGTTCTGTTCTTCAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCa                                                                                                        >  2:3533086/1‑150 (MQ=255)
                                               cTGGCCGTTCTGTTCTTCAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCAc                                                                                                       >  2:1996429/1‑150 (MQ=255)
                                                    cGTTCTGTTCTTCAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCgg                                                                                                  <  1:3009775/150‑1 (MQ=255)
                                                        ctgttctTCAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAAt                                                                                              >  1:241841/1‑150 (MQ=255)
                                                        ctgttctTCAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAAt                                                                                              <  2:71792/150‑1 (MQ=255)
                                                        ctgttctTCAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATAACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAAt                                                                                              >  2:1509124/1‑150 (MQ=255)
                                                             cttcAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTg                                                                                         <  2:2670588/150‑1 (MQ=255)
                                                             cttcAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTg                                                                                         <  2:1569269/150‑1 (MQ=255)
                                                               tcAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGgc                                                                                       <  1:830517/150‑1 (MQ=255)
                                                                cAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATAACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGgcg                                                                                      <  1:2075567/150‑1 (MQ=255)
                                                                    gCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTa                                                                                  <  2:3180941/150‑1 (MQ=255)
                                                                    gCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTa                                                                                  <  2:3340962/150‑1 (MQ=255)
                                                                        atatGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATAACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCcg                                                                              >  1:3920922/1‑150 (MQ=255)
                                                                             cAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAg                                                                         <  2:3930532/150‑1 (MQ=255)
                                                                             cAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAg                                                                         <  2:1092013/150‑1 (MQ=255)
                                                                              aaCGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGa                                                                        >  1:4577916/1‑150 (MQ=255)
                                                                              aaCGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCAACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGa                                                                        >  1:2196205/1‑150 (MQ=255)
                                                                               aCGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGAt                                                                       <  1:1423944/150‑1 (MQ=255)
                                                                                cGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGAtt                                                                      >  1:3521904/1‑150 (MQ=255)
                                                                                cGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGAtt                                                                      >  1:2156443/1‑150 (MQ=255)
                                                                                 gTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTg                                                                     <  1:785297/150‑1 (MQ=255)
                                                                                 gTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTg                                                                     >  2:193395/1‑150 (MQ=255)
                                                                                   aCCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTggg                                                                   >  2:1485883/1‑150 (MQ=255)
                                                                                    ccGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGa                                                                  <  1:1999885/150‑1 (MQ=255)
                                                                                        agaCGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCa                                                              >  1:1826032/1‑150 (MQ=255)
                                                                                         gaCGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCaa                                                             <  2:4357814/150‑1 (MQ=255)
                                                                                           cGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACa                                                           <  1:371604/150‑1 (MQ=255)
                                                                                            ggAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCCTCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACAt                                                          >  1:2515098/1‑150 (MQ=255)
                                                                                            ggAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACAt                                                          >  1:1028035/1‑150 (MQ=255)
                                                                                                 cACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCAATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGACGCGCAGATTGTGAGGCAACATTTGAt                                                     >  1:1784818/1‑150 (MQ=255)
                                                                                                  aCCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATc                                                    >  1:4212084/1‑150 (MQ=255)
                                                                                                  aCCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATc                                                    <  2:2643580/150‑1 (MQ=255)
                                                                                                   ccAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCg                                                   >  2:2358225/1‑150 (MQ=255)
                                                                                                    cAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGt                                                  <  2:1822504/150‑1 (MQ=255)
                                                                                                    cAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCCGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGt                                                  >  1:4614484/1‑150 (MQ=255)
                                                                                                        aTCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGt                                              >  2:1913513/1‑150 (MQ=255)
                                                                                                            tGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTg                                          >  2:2472284/1‑150 (MQ=255)
                                                                                                             gAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGa                                         >  2:4534612/1‑150 (MQ=255)
                                                                                                               aTACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGaca                                       >  1:2494745/1‑150 (MQ=255)
                                                                                                                  cGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTg                                    >  2:1753021/1‑150 (MQ=255)
                                                                                                                  cGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTg                                    <  2:4288227/150‑1 (MQ=255)
                                                                                                                  cGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTg                                    >  2:3440015/1‑150 (MQ=255)
                                                                                                                  cGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATAACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTg                                    <  2:4706080/150‑1 (MQ=255)
                                                                                                                    aaaCGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGcc                                  <  2:2627422/150‑1 (MQ=255)
                                                                                                                        ggTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGccacca                              >  1:1514115/1‑150 (MQ=255)
                                                                                                                        ggTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGccacca                              <  2:1818956/150‑1 (MQ=255)
                                                                                                                         gTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCaa                             <  1:4454394/150‑1 (MQ=255)
                                                                                                                           ccAGCCGGAACCTGCGTTTCGTAAATCACATGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACa                           <  1:4100600/150‑1 (MQ=255)
                                                                                                                           ccAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACa                           >  2:4172483/1‑150 (MQ=255)
                                                                                                                           ccAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCCGACACTGCCACCAACa                           <  1:2776832/150‑1 (MQ=255)
                                                                                                                                cGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAAtt                      <  1:263100/150‑1 (MQ=255)
                                                                                                                                cGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAAtt                      <  2:241841/150‑1 (MQ=255)
                                                                                                                                 ggAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTa                     <  1:3299607/150‑1 (MQ=255)
                                                                                                                                   aaCCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTATCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTaaa                   <  2:4134147/150‑1 (MQ=255)
                                                                                                                                   aaCCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTaaa                   >  2:4477233/1‑150 (MQ=255)
                                                                                                                                     ccTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAAc                 <  2:2747192/150‑1 (MQ=255)
                                                                                                                                     ccTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAAc                 <  1:3592723/150‑1 (MQ=255)
                                                                                                                                     ccTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAAc                 <  1:1061758/150‑1 (MQ=255)
                                                                                                                                          gTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTTATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCg            <  1:1380601/150‑1 (MQ=255)
                                                                                                                                          gTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCg            <  1:3204163/150‑1 (MQ=255)
                                                                                                                                           tttCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGa           <  1:2402907/150‑1 (MQ=255)
                                                                                                                                               gTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAat       <  1:2952367/150‑1 (MQ=255)
                                                                                                                                               gTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAat       >  1:3325886/1‑150 (MQ=255)
                                                                                                                                                tAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAata      <  2:1197910/150‑1 (MQ=255)
                                                                                                                                                 aaaTCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAACatat     >  2:3769352/1‑150 (MQ=255)
                                                                                                                                                 aaaTCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAatat     <  1:852404/150‑1 (MQ=255)
                                                                                                                                                   aTCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCCACATAATTAAAACCGTCGAAAATATCg   >  1:2467198/1‑150 (MQ=255)
                                                                                                                                                    tCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAATTGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCgg  <  2:579657/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGAGTGAGGTATGGCATCGAGGCGGTGCTGATGTCATATTCCTGCACCTGGCCGTTCTGTTCTTCAATGCGAATATGCAACGTACCGGAGACGGAATCACCAAGATCCTGAATACGAAACGGTCCAGCCGGAACCTGCGTTTCGTAAATCACACGCCCCATCTGGCTGACGGTCACTTTTGCTGTGGTGTGTGCCACGCCGGAAATGTCTGGCGCGTAGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGG  >  CP009273/2447387‑2447684

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: