Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,548,8230AC41.3% 54.2 / 12.7 48T6P (ACC→CCC) fimHtype 1 fimbriae D‑mannose specific adhesin
Reads supporting (aligned to +/- strand):  ref base A (21/6);  new base C (1/18);  total (22/26)
Fisher's exact test for biased strand distribution p-value = 7.66e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTACAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTA  >  NC_000913/4548577‑4549061
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
ataAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTACAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTAccc                                                                                                                                                                                                                                               <  1:409333/250‑1 (MQ=255)
             gaacgcgCAACTCATCCAGTTAGCGCTACAGGATGACAGTGGCAACACATTGAATCCTGGCGCAACCAAACAAGTTCAGGCGGATGATTCCTCACAATCAGCGCACTTACCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGCGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTActg                                                                                                                                                                                                                                 >  2:256162/4‑251 (MQ=255)
                        aCATCCAGTTAGAGCTACAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTCCCCTGTTTGCTGTACTGCTGATGGGCTg                                                                                                                                                                                                                     <  1:201392/251‑1 (MQ=255)
                            ccAGTTGGAGATACAGGATGAAAGGGGCAACACATTGAATACTGACGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGAACTTCACGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTAACCTGTTTGCTGGACTGCTGATGGGCTGGt                                                                                                                                                                                                                   >  2:264256/1‑249 (MQ=255)
                                             aTGACAGTGGCTACACATTGAATTCTGGCGCAACCAAACCAGTTCAGGTGGATGATTCATCACCATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGACACCATTCAGGCCGTGATTAGCATCACCTATCCCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTCTTCACATGTTTGCTGTACTGCTGCTGGGCTGGCCGGTAAATGCCTGGTCAtt                                                                                                                                                                                               >  2:171853/1‑251 (MQ=255)
                                             aTGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTGCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCAGTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGGAAATGCCTGGTCAtt                                                                                                                                                                                               <  2:405943/251‑1 (MQ=255)
                                                       cAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTaa                                                                                                                                                                                      >  2:201910/1‑250 (MQ=255)
                                                         acacATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGACACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAAc                                                                                                                                                                                   >  1:294583/1‑251 (MQ=255)
                                                           acaTTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTAATTAGCATCACCTATACCTACAGCTGAACCCGCAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTCAATGCCTGGTCATTCGCCTGTAAAACCg                                                                                                                                                                                 >  1:248507/1‑251 (MQ=255)
                                                                       ggcgaaacaaaaccatttCAGGGGGATGATTCCCCAAAATCAGCGCACTCCCCGTTCCAGATAAGAGAATTGACAGTAAAGGGCGGGGCCCCTCAGGGACCCATTCAGGCTGTGATTAGCATCACCTATACCTCCAGCTGACCCCGAAGAGATGATGGTAATGAAACGAGTTATTCCCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGCCATCCGCCTGTAAAACCGCCAATGGTAc                                                                                                                                                                       <  2:122958/233‑1 (MQ=255)
                                                                           aaacaaaaaCAGTTCAGGTGGATGATCCCTCACAATCAGCGCACTCCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATCCCTACAGCTGACCCCGAAGAGATGATTGTAATGAAACGAGTTATTCCCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTAt                                                                                                                                                                 <  1:178870/246‑1 (MQ=255)
                                                                                 aaaCAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTAt                                                                                                                                                           >  1:130578/1‑251 (MQ=255)
                                                                                 aaaCAATTCAGGTGGATGATTCCCCAAAAACAGCGCACTCCCCTTTTCAGGTAAGAGCATGGACTGAAAATGGCGGAGCCACTCAGGGACCCATTCAGTCAGTGATTAGCATCCCCTATACCTACAGTTGAACCCGAAGAGAGGGTTGTAATGAAACGATTTATTCCCCTGTTTGCTGTACTGCTGAGGGGCTGGCCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTAt                                                                                                                                                           <  2:86121/251‑1 (MQ=255)
                                                                                     aGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGc                                                                                                                                                       >  1:179053/1‑251 (MQ=255)
                                                                                      gTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCg                                                                                                                                                      >  1:284849/1‑251 (MQ=255)
                                                                                      gTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCg                                                                                                                                                      >  1:101434/1‑251 (MQ=255)
                                                                                      gTTCAGGTGGATAATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGTGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGCGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCg                                                                                                                                                      >  1:318691/1‑251 (MQ=255)
                                                                                          ggggggATGATTCCTCACAATCAGCGCACTCCCCGTTACAGGTCAGAGCATTGACAGTAAAGGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCCCCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGATTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCTTGGTCTTTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTgg                                                                                                                                                  <  1:92215/247‑1 (MQ=255)
                                                                                          aGGTGGATGACTCCTCACAATCAGCGAACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGAAGCCACTCAGGGAACCCTTCAGGCAGTGATTAGCATCACCTATATCTACCGCTGAACCCGAAGAGATGATTGTAAGGAACCGAGTAATAACCAGATTTGCTGGACTGCTGATGGTCGGGTCGGTAAATGCTTGGTCATTCGCCTGTAACACCACCAATGACACCGCTATCCCTATTGACGGTgg                                                                                                                                                  >  2:224332/1‑251 (MQ=255)
                                                                                            gTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGc                                                                                                                                                 >  1:250125/1‑250 (MQ=255)
                                                                                             tGGATGATTCCTCACAATCAGCGCACTCCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTCCAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTggtcc                                                                                                                                               >  1:42293/1‑248 (MQ=255)
                                                                                                atgatTCATCACAATCAGCGCCCTTCCCGTTACAGGTCAGCGCATTGACAGTAAATGGCGGAGCCACTCAGTGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAATGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCACTGGTACCGCTATCCCTATTGGCGGTGGccgacc                                                                                                                                            >  1:322444/1‑246 (MQ=255)
                                                                                                  gatTCCTCACAATCAGCGCCCTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGAAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTACATGCCTGGTCATTCGCCTGTAAAACCGCCACTGGTACCGCTATCCCTCTTGGCGGttggcgctcc                                                                                                                                           >  2:95263/1‑241 (MQ=255)
                                                                                                  gatTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCGGCGccc                                                                                                                                          >  1:77574/1‑250 (MQ=255)
                                                                                                                  cgcATTTCCCGTTACAGGTAAGAGCATTGAAAGAAAATGGAGGAGCCACTCAGGGAACCATCCAGGAATTGATTAGAATCACTTATACCTACAGCTGACCCCGAAGAGATGATGGTAATGAACCGAGTTATTACCCTGTTTGCTGTACTTCTGATGGGCTGGTCGGTAATGGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCt                                                                                                                          <  1:284059/251‑1 (MQ=255)
                                                                                                                    cACTTCCCGTTACAGGTCAGAGCATTGAAAGTAAATGGCGGAGCCACTCAGGGACCCTTTCAGGCAGTGATTAGCATCCCCTATACCTACAGCTGACCCAGAGGAGATGATTGTAATAAAACGAGTTATTCCCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCTTTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTg                                                                                                                        <  2:109714/251‑1 (MQ=255)
                                                                                                                              tacagtccAGAGCATTGACAGTAAAGGGCGGAGCCACCAAGGGACCCATTCAGGCAGTGATTAGCATCACTTATACCTACAGCTGACCCCGAAGGGATGATTGTAATAAAACGAGTTATTCCCCTGTTGGCTGTACTGTTGATGGGCTGGTCGGTAAATGCTTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCcgtcg                                                                                                               <  1:117343/243‑1 (MQ=255)
                                                                                                                                  ggTCAGAGCATTGCCAGTAACTGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTCCAGCTGAACACGACGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCCGCGCCACTGTTTAAGTAAACCCTGCGCCCCTCGTgttt                                                                                                          >  2:73562/1‑248 (MQ=255)
                                                                                                                                     cAGATCATGGACATTAAAGGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATCCCTCCAGCTGAACCCGAAGTGATGATTGTAATGAACCGAGTTATTCCCCTGTTTGCTGTACTGCTGATGGGCTGGTCGTTAATTGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCTCCCGTCGTGAAtgt                                                                                                        <  1:118933/250‑1 (MQ=255)
                                                                                                                                           aTTGACAGTAAAGGGAGGCGCCTCCCAGGAAACAATGCAGGCGGTGATTACCATCCCCTATCCCTCCAGTGGACCCTAAATAGATGATTGTAAGGAAACGATTAATCTCCCTGTTCGCGTTACGGCTAAGGGTCTGGTCGGAAATTGCCTGGTCTTTCGCCAGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAGCCTTGCGCCCGTCGTGAATGTGGGGCaa                                                                                                 <  1:171853/251‑1 (MQ=255)
                                                                                                                                              gACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCCGCGCCAATGGTTTTGTAAACCTTTCTCCCGCCGGGGATGTggggggaaa                                                                                               >  2:346328/1‑245 (MQ=255)
                                                                                                                                                   aaaaTGGCGGAGCCACTCAGGGAACCATTCAGGCAGTAATAAGCATCACCTATCCCTACAGCGGAACCCGAAGAGATGATTGTAAGGAAACGAGTTATTCCCCTGTTTGTTGTACTGCTGATGGGCTGGTCGGTAAATGCTTGGTCATTCCCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGGGAATGTGGGGCAAAACCTGGt                                                                                         <  1:324419/250‑1 (MQ=255)
                                                                                                                                                    aaaTGGAGGAGCCACTCAGGGAACCATGCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGGTTGGAATGAAACGAGTTATTCCCTTGTTTGCTGTCCTGCGGATGGGCTGGTCGTTAAATGCGGGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTc                                                                                        <  1:61098/251‑1 (MQ=255)
                                                                                                                                                                cccTCAGGGAACCATCCAGGAAGTGATTAGCATCCCCTATCCCTACAGCTTAACCCGAAGAGATGATTGTAATGAAACGAGTTATTCCCCTGTTGGCTTTCCTGCTGATGGGCTGTTCGGTAAATGCCTGTTCATTCGCCTGAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATTTGGGGCAAAACCTGGTCGTGGATCTTTCg                                                                            <  1:346328/249‑1 (MQ=255)
                                                                                                                                                                    cagggatctaatCAGGCAGAGATGAGCATCGCCTATCCCTACATCTAACCCCGAAGAGAGGATGGTAATGAAACGAGTTATTCCCCTGTTTGGTGTACTGTGGAGGGTGGGGTCGGAAAATGCTGGCTCATTCGCCGGAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCGTGCTCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGAc                                                                          <  1:257740/238‑1 (MQ=255)
                                                                                                                                                                                   gcagtgttagacctCACCTATTCCTACAGCTAAACCCGAAGAGATTATTGAAATGAAAGGGTTTATTCCCCTGTTTGCTGTACTGCTGAGGGGCGGGTCGGTAAAGGCCTGTTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATa                                                         <  2:130578/238‑1 (MQ=255)
                                                                                                                                                                                             gaatccccTATACCTACAGCTGAACCCAAAGAGATGATTGTAATGAAACGAGTTATTCCCCTGTTTGCTTTACTGCGTAGGGGCGGGTCGGTAAATCCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATcc                                               <  1:131444/245‑1 (MQ=255)
                                                                                                                                                                                                  cctgttaatatcatgtGAACCCTGAAGGGATGATGTTAATGCACCGATTTATACCCCTGTTTTCTGTACTTCGGAGGGCCTGGTCTGAAAACTCCGGGTCTTCCGCCTTTAAAACCGCCAATGGTACCGCTATCCCTACTGGCGGTGGCAGCGCCAATGTTTATTTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACTCAAATCTGTTGCCATAACGATTAGCCGGaa                                           <  2:294583/236‑1 (MQ=255)
                                                                                                                                                                                                      atcccTACAGCTGACCCCGAAGAGAGGATTGAAATGAAACGTGTTATTCCCCTGTTTGCTGTTCGCCGTAGGGGCTGTTCGGTAAAGCCCTGGCCATCCGCCTGAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGGGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAAGCGTTTGCCATAACGATTATCCTGAAAcc                                        <  2:161573/246‑1 (MQ=255)
                                                                                                                                                                                                      ataCCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCCCAGTTTATTTAAAACTTGCGCCCGTCGTGGATGTGGGGCAAAAACTGGTCGTGGCTCTTTTTCCGCAAATATTTTTCCCGAACAATTATCCCGTAACCAt                                      >  2:248872/1‑251 (MQ=255)
                                                                                                                                                                                                        cccTACAGCTGACCCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCGGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTa                                    <  1:393754/250‑1 (MQ=255)
                                                                                                                                                                                                         ccTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCCATGTTTATGTTAACCTTTCGCCCGTCGTGAATGTGGGGCCAAACCTGGGCGTGGGTCTTTCCACGCAAATCTTTTTCCATAACGATTATCCGGAAACCATTAc                                   >  1:55933/1‑251 (MQ=255)
                                                                                                                                                                                                                   aaCCCGAAGAGATGATTTTAATAAAACGATTTTTTTCCCTGTTTGCGTTCCTGTTAAGGGGCTGTTCGGAAAATTCTTGGTCTTCCCCCTGAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAGCCATTACAGACTATGTc                         <  1:351177/251‑1 (MQ=255)
                                                                                                                                                                                                                   aaCCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCCGCGCCAATGTTTTTGTAAAACTTTCCCCCCTCGGGAATGTTGGGCAAAACCTGGTCGTGGATCTTTTCACGCCAAACTTTTTCCCTAACGATTTTCCCGAAAcccataaaaactatttc                         >  1:282577/1‑236 (MQ=255)
                                                                                                                                                                                                                     ccccgaatagagtatTGAAATAAAACGAGTTATTACCGTGTTTGCTGTACGGCGGATGGGCTGTTGGGTAAAGTCCTGGTCATCCGCCTGAAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTca                        <  2:211418/238‑1 (MQ=255)
                                                                                                                                                                                                                           agagATGATTGTAATGAAACGAGTTATTCCCCTGTTTGCTGTACTGTTGAGGGGCTGGTCGGTAAAGGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGc                  <  1:201910/251‑1 (MQ=255)
                                                                                                                                                                                                                                          gAAACGAGTTATTCCCTTGTTTGCTTTACTGCTGAGGGGCTGGGCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCt    <  1:131409/250‑1 (MQ=255)
                                                                                                                                                                                                                                           aaaCGAGTTATTCCCCTGTTTGCTGTACTGCTGATGGTCTGGGCGGTAAATCCCTGGTTATCCCCTTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTa  <  2:146927/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
ATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTACAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGAGCATTGACAGTAAATGGCGGAGCCACTCAGGGAACCATTCAGGCAGTGATTAGCATCACCTATACCTACAGCTGAACCCGAAGAGATGATTGTAATGAAACGAGTTATTACCCTGTTTGCTGTACTGCTGATGGGCTGGTCGGTAAATGCCTGGTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTA  >  NC_000913/4548577‑4549061

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