Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913306,0230GA56.7% 2.9 / 84.1 61T265T (ACC→ACT) ecpDfimbrial adhesin EcpD
Reads supporting (aligned to +/- strand):  ref base G (16/10);  new base A (18/16);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.67e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GATATCGCTGATCCCGAGGTTCACCATCTGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGG  >  NC_000913/305780‑306273
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                          
gATATCGCTGATCCCGAGGTTCACCATCTGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCTGACGGCCGATGCTAACGACGAGTTGGCGATCTCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGAGTATAg                                                                                                                                                                                                                                                      <  1:173074/250‑1 (MQ=255)
gATATCGCTGATCCCGAGGTTCACCATCTGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGAGTATAg                                                                                                                                                                                                                                                      <  1:112933/250‑1 (MQ=255)
 atatCGCTGATCCCGAGGTTCACCATCTGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGAGTATAGtt                                                                                                                                                                                                                                                    <  1:83377/251‑1 (MQ=255)
   atCGCTGATCCCGAGGTTCACCATCTGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCTAGTCTGCGCTCTTCATCTCTTTGAAGGTGTATTTATAGGTGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAAGTGCATATCGTATGCCCCGACGGTTGAGGCTAACGACGAGGTGGTGATCGCGGTGAAGATATGGATTGAGGTGTTGCTGATTCCGTTTGTTTACAGTGTATAGTTaa                                                                                                                                                                                                                                                  <  2:275661/251‑1 (MQ=255)
        ggATCCCGAGGTTCACCATCTGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGAGTATAGTTAACCATc                                                                                                                                                                                                                                             <  2:119660/250‑1 (MQ=255)
         gATCCCGAGGTTCACCATCTGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGAGTATAGTTAACCATc                                                                                                                                                                                                                                             <  1:218976/250‑1 (MQ=255)
              cGAGGTTCACCATCTGCTTAAAGAAGTTGCTCGAGAAGAAAACTTAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGAGTATAGTTAACCATCTTACa                                                                                                                                                                                                                                        >  1:212142/1‑250 (MQ=255)
              cGAGGTTCACCATCTGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGCGGTGGTGCTGAGTCCGTTTGTTTGCAGAGTATAGTTAACCATCTTACa                                                                                                                                                                                                                                        >  2:115558/1‑250 (MQ=255)
               gAGGTTCACCATCTGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCTATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGAGTATAGTTAACCATCTTACAAc                                                                                                                                                                                                                                      <  1:274807/251‑1 (MQ=255)
                aGGTTCACCATCTGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACt                                                                                                                                                                                                                                     >  1:145280/1‑251 (MQ=255)
                aGGTTCACCATCTGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACt                                                                                                                                                                                                                                     <  1:18393/251‑1 (MQ=255)
                      aCCATCTGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATcc                                                                                                                                                                                                                               >  1:154485/1‑251 (MQ=255)
                       ccATCTGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGAGTATAGTTAACCATCTTACAACTTAATCCt                                                                                                                                                                                                                              >  1:790/1‑251 (MQ=255)
                            tGCTTAAAGACGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATTGATTTAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACg                                                                                                                                                                                                                         >  1:166333/1‑251 (MQ=255)
                            tGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGGATTGATCAACGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGGTGCTGAGTACGTTTGTTTGCAGGGTATAGTTAACCATCTCACAACTTAATCCTGaa                                                                                                                                                                                                                           >  1:166361/1‑249 (MQ=255)
                             gCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGAGTATAGTTAACCATCTTACAACTTAATCCTGAAc                                                                                                                                                                                                                          <  1:256964/249‑1 (MQ=255)
                                tAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGc                                                                                                                                                                                                                       <  1:72061/249‑1 (MQ=255)
                                   agaagTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCg                                                                                                                                                                                                                  <  1:68263/251‑1 (MQ=255)
                                          gCTCGAGAAGAAAACATAGAGCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGGTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGttt                                                                                                                                                                                                           <  2:229144/251‑1 (MQ=255)
                                              gagaAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGAGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGt                                                                                                                                                                                                       >  1:30312/1‑251 (MQ=255)
                                                          tAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGAGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGc                                                                                                                                                                                           >  2:25869/1‑251 (MQ=255)
                                                             aTCGAGTCTGCGCTCTTCATCTCGTTGAAGGGGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTTAGTCCGTTTGTTTGCAGGGGATAGTTAACCATCTTACAACTTACTCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGcatt                                                                                                                                                                                        >  2:62589/1‑248 (MQ=255)
                                                               cGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTcc                                                                                                                                                                                      >  2:162960/1‑251 (MQ=255)
                                                                gAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGAGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTccc                                                                                                                                                                                     >  2:298269/1‑251 (MQ=255)
                                                                gAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGAGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTAGGCGTTccc                                                                                                                                                                                     >  1:208569/1‑251 (MQ=255)
                                                                       cgcTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGccc                                                                                                                                                                              >  1:274839/1‑251 (MQ=255)
                                                                                 ctcGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGAGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGt                                                                                                                                                                     >  1:144245/1‑250 (MQ=255)
                                                                                   cGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGAGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTa                                                                                                                                                                    <  2:18602/249‑1 (MQ=255)
                                                                                    gTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGAGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTc                                                                                                                                                                 >  1:229953/1‑251 (MQ=255)
                                                                                        aaGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGCCGGAACTGCATATCGTACGCCCCGACGGCCGAGGCTCACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGGTGCTGAGTCCGTTTGTTTGCAGAGTATAGTTAACCATCTTACAAATTAATCCTGACCGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGGACTCAGTc                                                                                                                                                             >  2:12862/1‑251 (MQ=255)
                                                                                          ggTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGATAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCACAGCGTCGGTACTCCGTCGc                                                                                                                                                           >  2:284449/1‑251 (MQ=255)
                                                                                             gtagtaATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGAGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTg                                                                                                                                                         <  1:298269/250‑1 (MQ=255)
                                                                                               agtaATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGCGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGAGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCACAGAGTCGGTACTCCGTCGCTGTTg                                                                                                                                                      >  1:278727/1‑251 (MQ=255)
                                                                                                 taATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGAGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGAt                                                                                                                                                    >  1:73091/1‑251 (MQ=255)
                                                                                                               gCTCACCGCTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGAGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGGCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCa                                                                                                                                      <  2:94191/251‑1 (MQ=255)
                                                                                                                      ggTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTg                                                                                                                               >  2:274824/1‑251 (MQ=255)
                                                                                                                               gAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGAGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTAt                                                                                                                      >  2:135312/1‑251 (MQ=255)
                                                                                                                                  cTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGATTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGACAGAGCGTGGGGAATTATCa                                                                                                                    >  2:12760/1‑250 (MQ=255)
                                                                                                                                     ccATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCg                                                                                                                >  2:120359/1‑251 (MQ=255)
                                                                                                                                       aTTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGAGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGt                                                                                                               >  1:15908/1‑250 (MQ=255)
                                                                                                                                          cataCTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGt                                                                                                               >  1:215607/4‑247 (MQ=255)
                                                                                                                                          cAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGAGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGAt                                                                                                           <  1:45919/251‑1 (MQ=255)
                                                                                                                                                      cATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCg                                                                                               <  1:295713/251‑1 (MQ=255)
                                                                                                                                                               cGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACGCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGa                                                                                      <  2:154485/251‑1 (MQ=255)
                                                                                                                                                                      aCGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCg                                                                               >  1:150100/1‑251 (MQ=255)
                                                                                                                                                                        ggCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGAGTATAGTTAACCATCTTACAACGTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGaa                                                                             <  2:64415/251‑1 (MQ=255)
                                                                                                                                                                         gCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGAGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAAc                                                                            <  1:25869/251‑1 (MQ=255)
                                                                                                                                                                          ccGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGAGAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACa                                                                           <  1:293488/251‑1 (MQ=255)
                                                                                                                                                                              ggCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGAGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACAGTAGATTTGCTGCTTTCGTATGGGTGAAGTTGCGAACATAcc                                                                       >  1:38172/1‑251 (MQ=255)
                                                                                                                                                                                             ggCGATCGCCGGGAATATATGGATTGAGGTGTGGCTGAGTCCGTTTGTTTGCAGATTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGTGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGcc                                                         <  2:278727/250‑1 (MQ=255)
                                                                                                                                                                                                gATCGCCGGGAATATCTGGATTGAGGTGTTGCTGAGTCCGTTTGCCTGCAGAGCATAGTTAACCATCTCACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGAt                                                     >  1:304033/1‑251 (MQ=255)
                                                                                                                                                                                                                         gtgtTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTc                                                                                                                                           >  2:106553/1‑140 (MQ=255)
                                                                                                                                                                                                                         gtgtTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTc                                                                                                                                           <  1:106553/140‑1 (MQ=255)
                                                                                                                                                                                                                         gtgtTGCTGAGTCCGTTTGTTTGCAGAGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGc                            <  1:135312/251‑1 (MQ=255)
                                                                                                                                                                                                                                  aGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGtt                   >  1:174959/1‑251 (MQ=255)
                                                                                                                                                                                                                                    tCCGTTTGTTTGCAGAGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTc                  >  1:256991/1‑250 (MQ=255)
                                                                                                                                                                                                                                       gtttgtttGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACag              <  1:162960/251‑1 (MQ=255)
                                                                                                                                                                                                                                        tttgtttgCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTcac                >  1:177127/1‑248 (MQ=255)
                                                                                                                                                                                                                                              tgCAGGGTATAGGTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCTCCGCTGCTGGCGTCATAGTTCACAGAGGTCt        <  2:129181/250‑1 (MQ=255)
                                                                                                                                                                                                                                              tgCAGAGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTg       >  2:220775/1‑251 (MQ=255)
                                                                                                                                                                                                                                               gCAGAGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTg       <  2:48632/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   ggTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTggcagg  >  2:228262/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   ggTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTggcagg  >  1:33174/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                          
GATATCGCTGATCCCGAGGTTCACCATCTGCTTAAAGAAGTTGCTCGAGAAGAAAACATAGATCGAGTCTGCGCTCTTCATCTCGTTGAAGGTGTAGTAATAGGCGGTATTGCTCACCGGTTTCCATGAACTGCCATTCAGACTGAACTGCATATCGTACGCCCCGACGGCCGAGGCTAACGACGAGTTGGCGATCGCCGGGAATATATGGATTGAGGTGTTGCTGAGTCCGTTTGTTTGCAGGGTATAGTTAACCATCTTACAACTTAATCCTGAACGGCTGCCGATGGTTTGCGTCCGGCAGTCGGCGTTCCCTTCGCCCAGAGTCGGTACTCCGTCGCTGTTGATAAATACTTCCGCCAGCGAGTGGGTATTTATCAACCGTAGATTTGCTGCTTTCGTATGGGTGACGTTGCGAACATACCAGTTGCCGGAGGCCTGATCCTTACAGCGTGCGCCGCTGCTGGCGTCATAGTTCACAGAGGTCTGGCAGG  >  NC_000913/305780‑306273

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