Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,288,032 C→A 6.6% intergenic (+22/‑58) yraH → / → yraI putative fimbrial‑like adhesin protein/putative periplasmic pilin chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,288,0320CA 6.6% 704.2 / 28.4 198intergenic (+22/‑58)yraH/yraIputative fimbrial‑like adhesin protein/putative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base C (94/91);  new base A (9/4);  total (103/95)
Fisher's exact test for biased strand distribution p-value = 2.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.64e-01

Only 100 of 201 total aligned reads displayed.

CCCGGCGGAATCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACA  >  NC_000913/3287884‑3288175
                                                                                                                                                    |                                                                                                                                               
    gcggAATCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTc                                                                                                                                             <  1:55851/149‑1 (MQ=255)
        aaTCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGAt                                                                                                                                         <  2:1554277/149‑1 (MQ=255)
        aaTCCAATCGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGAt                                                                                                                                         <  2:1372397/149‑1 (MQ=255)
          tCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATaa                                                                                                                                       >  1:78698/1‑149 (MQ=255)
          tCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATaa                                                                                                                                       <  4:1165192/149‑1 (MQ=255)
          tCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATaa                                                                                                                                       >  3:855641/1‑149 (MQ=255)
            cAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACAGGTCAGATAAc                                                                                                                                      >  4:869173/1‑148 (MQ=255)
            cAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACAGGTCAGATAAc                                                                                                                                      <  3:869173/148‑1 (MQ=255)
             aaTGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGc                                                                                                                                    >  3:1232343/1‑149 (MQ=255)
             aaTGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGc                                                                                                                                    >  3:361297/1‑149 (MQ=255)
                  tATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCttt                                                                                                                               >  3:1165152/1‑149 (MQ=255)
                      gTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGATGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCa                                                                                                                                            >  4:904500/1‑132 (MQ=255)
                       tCAATGTTGACCTCATCGCTTACTACGTATCAACCGTGGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTttcctt                                                                                                                          >  1:265891/1‑149 (MQ=255)
                        cAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTttccttc                                                                                                                         <  3:705490/149‑1 (MQ=255)
                           tGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCctc                                                                                                                      >  4:1433443/1‑149 (MQ=255)
                           tGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCctc                                                                                                                      <  4:884116/149‑1 (MQ=255)
                               gACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTtc                                                                                                                  >  3:97698/1‑149 (MQ=255)
                                 ccTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTtctc                                                                                                                <  4:1541649/149‑1 (MQ=255)
                                  cTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCg                                                                                                               >  2:1108495/1‑149 (MQ=255)
                                   tCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGt                                                                                                              >  1:1108126/1‑149 (MQ=255)
                                      tCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGt                                                                                                           <  4:278566/149‑1 (MQ=255)
                                      tCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGt                                                                                                           >  3:1493501/1‑149 (MQ=255)
                                       cGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTa                                                                                                          >  3:162858/1‑149 (MQ=255)
                                        gCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGtt                                                                                                             <  4:444809/145‑1 (MQ=255)
                                          ttACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATaa                                                                                                                                       >  4:621757/1‑117 (MQ=255)
                                          ttACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATaa                                                                                                                                       <  3:621757/117‑1 (MQ=255)
                                          ttACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCa                                                                                                       >  2:984473/1‑149 (MQ=255)
                                              tacGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGtt                                                                                                             <  3:896735/139‑1 (MQ=255)
                                              tacGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGtt                                                                                                   <  2:667643/149‑1 (MQ=255)
                                              tacGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGtt                                                                                                   >  1:402935/1‑149 (MQ=255)
                                                    tCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAAc                                                                                                                                      >  2:454523/1‑108 (MQ=255)
                                                      aaCCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATa                                                                                           >  3:926813/1‑149 (MQ=255)
                                                       aCCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATaa                                                                                          >  1:1275236/1‑149 (MQ=255)
                                                          gTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGtt                                                                                                   <  1:156967/137‑1 (MQ=255)
                                                          gTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTa                                                                                       >  2:1242986/1‑149 (MQ=255)
                                                              ttCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTc                                                                                   >  1:1415931/1‑149 (MQ=255)
                                                                ccATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCaa                                                                                 >  2:20121/1‑149 (MQ=255)
                                                                   tCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGaa                                                                                                >  3:1143773/1‑131 (MQ=255)
                                                                   tCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAAc                                                                              >  2:1455804/1‑149 (MQ=255)
                                                                   tCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAAc                                                                              >  2:1124096/1‑149 (MQ=255)
                                                                   tCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACAGGTCAGATAAc                                                                                                                                      >  2:1201700/1‑93 (MQ=255)
                                                                          tCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCaaaa                                                                               <  3:1389736/141‑1 (MQ=255)
                                                                          tCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGaa                                                                           >  1:1208776/1‑145 (MQ=255)
                                                                              cccAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTcgcg                                                                   >  2:222450/1‑149 (MQ=255)
                                                                               ccAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCgg                                                                  >  1:961123/1‑149 (MQ=255)
                                                                                cAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATGCGTTAGGGCGTTGTCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGt                                                                 >  4:1153732/1‑149 (MQ=255)
                                                                                cAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGt                                                                 <  4:1382446/149‑1 (MQ=255)
                                                                                        gCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTa                                                                                       <  4:1185708/119‑1 (MQ=255)
                                                                                        gCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTaa                                                         >  3:196659/1‑149 (MQ=255)
                                                                                         cAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAAc                                                        <  4:732459/149‑1 (MQ=255)
                                                                                           gACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCt                                                      >  1:1019042/1‑149 (MQ=255)
                                                                                           gACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCt                                                      >  2:517559/1‑149 (MQ=255)
                                                                                            aCGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCtt                                                     <  1:1495687/149‑1 (MQ=255)
                                                                                             cGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCttg                                                    >  4:1465525/1‑149 (MQ=255)
                                                                                              gCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCttgt                                                   <  3:609568/149‑1 (MQ=255)
                                                                                              gCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCttgt                                                   <  1:20121/149‑1 (MQ=255)
                                                                                               cTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCttgtt                                                  <  1:1455804/149‑1 (MQ=255)
                                                                                               cTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCttgtt                                                  <  1:517559/149‑1 (MQ=255)
                                                                                                tATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCttgttg                                                 >  3:459584/1‑149 (MQ=255)
                                                                                                   cGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTgtgt                                              <  3:1465525/149‑1 (MQ=255)
                                                                                                       aacgtaacgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCt                                          <  2:265891/149‑1 (MQ=255)
                                                                                                        acgtaacgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGtt                                                                                                             >  4:1133446/1‑81 (MQ=255)
                                                                                                          gtaacgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCt                                       <  1:1124096/149‑1 (MQ=255)
                                                                                                          gtaacgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCt                                       <  1:1187181/149‑1 (MQ=255)
                                                                                                            aacgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGt                                     >  2:1355966/1‑149 (MQ=255)
                                                                                                               gCTGGATTACCGTTAATACGTTAGGGCGTTATCAGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTAtt                                  <  2:484230/149‑1 (MQ=255)
                                                                                                                cTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTg                                 <  1:467190/149‑1 (MQ=255)
                                                                                                                  ggATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGc                               >  2:1589202/1‑149 (MQ=255)
                                                                                                                        ccGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGgcg                         >  3:1514596/1‑149 (MQ=255)
                                                                                                                         cGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGgcgc                        >  4:180043/1‑149 (MQ=255)
                                                                                                                         cGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGgcgc                        <  2:402935/149‑1 (MQ=255)
                                                                                                                          gTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCt                       >  4:158785/1‑149 (MQ=255)
                                                                                                                           ttAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGa                                                               >  3:801121/1‑108 (MQ=255)
                                                                                                                           ttAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCttgt                                                   >  2:1447372/1‑120 (MQ=255)
                                                                                                                              aTACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCt                                                      >  2:12759/1‑114 (MQ=255)
                                                                                                                              aTACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCa                   <  4:102482/149‑1 (MQ=255)
                                                                                                                               tACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTCTCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCt                       >  3:288920/1‑144 (MQ=255)
                                                                                                                               tACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCt                       <  4:288920/144‑1 (MQ=255)
                                                                                                                               tACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCag                  >  2:1211506/1‑149 (MQ=255)
                                                                                                                               tACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCag                  <  3:1003545/149‑1 (MQ=255)
                                                                                                                                 cGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCagga                <  4:1493501/149‑1 (MQ=255)
                                                                                                                                  gTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCaggag               <  2:584285/149‑1 (MQ=255)
                                                                                                                                  gTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCaggag               <  2:738113/149‑1 (MQ=255)
                                                                                                                                   ttAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCa                             >  4:848284/1‑134 (MQ=255)
                                                                                                                                     aGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAAGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGaa            >  3:1052947/1‑149 (MQ=255)
                                                                                                                                      gggCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAAt           >  1:162805/1‑149 (MQ=255)
                                                                                                                                        gCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCg         <  3:582082/149‑1 (MQ=255)
                                                                                                                                         cGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGa                                                               <  4:464889/94‑1 (MQ=255)
                                                                                                                                          gTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGtt       <  1:1355966/149‑1 (MQ=255)
                                                                                                                                           ttATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAg                            <  1:606572/127‑1 (MQ=255)
                                                                                                                                           ttATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAg                            >  2:606572/1‑127 (MQ=255)
                                                                                                                                           ttATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGttt      <  2:1382984/149‑1 (MQ=255)
                                                                                                                                           ttATCTGACAGGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAGAACGAACATTCGCGGTGATATTAACCTATTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCTtttt     >  4:1361544/1‑149 (MQ=255)
                                                                                                                                           ttATCTGACAGGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGttt      >  1:1314756/1‑149 (MQ=255)
                                                                                                                                           ttATCTGACAGGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGttt      >  2:823233/1‑149 (MQ=255)
                                                                                                                                              tCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTAtt                                  >  3:297879/1‑118 (MQ=255)
                                                                                                                                               cTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACa  >  2:1061966/1‑149 (MQ=255)
                                                                                                                                               cTGACAGGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACa  <  4:349343/149‑1 (MQ=255)
                                                                                                                                                    cTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGa                                                               <  2:122248/83‑1 (MQ=255)
                                                                                                                                                    cTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGa                                                               >  1:122248/1‑83 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CCCGGCGGAATCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACA  >  NC_000913/3287884‑3288175

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