Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,546,601 C→T I502I (ATC→ATT)  fimD → type I fimbriae usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,546,6010CT96.6% 277.9 / 0.8 91I502I (ATC→ATT) fimDtype I fimbriae usher protein
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (42/44);  total (44/47)
Fisher's exact test for biased strand distribution p-value = 6.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.89e-01

GTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAA  >  NC_000913/4546357‑4546846
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
gTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAacc                                                                                                                                                                                                                                                 >  1:816499/1‑250 (MQ=255)
  gCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTCCAACATTG‑AAacacc                                                                                                                                                                                                                                               >  2:273411/1‑250 (MQ=255)
  gCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAacac                                                                                                                                                                                                                                                <  1:17863/250‑1 (MQ=255)
                          aaCATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGAACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGcc                                                                                                                                                                                                                       >  1:472989/1‑251 (MQ=255)
                              tGGGGACACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCTAATGAATGCCCACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGaa                                                                                                                                                                                                                    >  2:836753/1‑250 (MQ=255)
                                  ggCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTc                                                                                                                                                                                                                >  1:299163/1‑250 (MQ=255)
                                    cACTTGGCGCTCTGTCTGTGGATATGATGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTTCGTTTTCTCTATAACAAATCGCTCAATGAATCAGTCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACAGTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCAcc                                                                                                                                                                                                             <  2:578588/251‑1 (MQ=255)
                                    cACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTGATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCTCGAATATTCAGTGAGTGGGTTACCGTTAGTCGACGAGCTGATATTTTAATTTCGCTGATACAACATACAGTCGGATGAATGGCTACAACAGGG‑AAACACAGGACGGAGTTATTCAGGTTAATCCGAAATTCAcc                                                                                                                                                                                                             <  2:331867/251‑1 (MQ=255)
                                     ccTGGGTGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCg                                                                                                                                                                                                            <  2:590615/250‑1 (MQ=255)
                                     aCTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGACACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCg                                                                                                                                                                                                            >  2:478632/1‑251 (MQ=255)
                                      cTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGACAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATCTTCAGTTAGTGGGTTACCGTTATTCGACCAACGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAAAATTA‑AAACACAGGCCGGAGTTATTCAGGTAAAGCCGAAATTCACCg                                                                                                                                                                                                            >  2:307548/1‑250 (MQ=255)
                                        gggCGCTCTGTCTGGGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGTACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGa                                                                                                                                                                                                           <  2:438978/249‑1 (MQ=255)
                                         ggCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATACCAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACt                                                                                                                                                                                                         >  1:854741/1‑250 (MQ=255)
                                         ggCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACt                                                                                                                                                                                                         >  1:66651/1‑250 (MQ=255)
                                           cgcTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGCCGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTAtt                                                                                                                                                                                                      >  2:370876/1‑251 (MQ=255)
                                           cgcTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTAtt                                                                                                                                                                                                      <  1:38567/251‑1 (MQ=255)
                                            gcTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTATATAACAACTCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAac                                                                                                                                                                                                                                                  >  1:388623/1‑206 (MQ=255)
                                            gcTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGGGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGGTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGCCGAATGAATGGCTACAACATTG‑AAac                                                                                                                                                                                                                                                  <  2:388623/206‑1 (MQ=255)
                                                 gtctgtGGATATGACGCAGGCTAATTCCACACTTACCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACCTACAGTCGAATGAATGGCTACAACATTG‑AAACACATGACGGAGTTATTCAGGTTAAGCCGAACTTCACCGACTATAACAAcc                                                                                                                                                                                                >  1:194076/1‑251 (MQ=255)
                                                     gtgGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTc                                                                                                                                                                                              <  1:360283/249‑1 (MQ=255)
                                                     gtgGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGc                                                                                                                                                                                            >  1:232535/1‑251 (MQ=255)
                                                             gACGCAGGCTAATTCCACAATTCCCGATGACAGTCAGCATTACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATTGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCCCCGACTATTACAACCTCGCTTATAACa                                                                                                                                                                                    <  2:66108/251‑1 (MQ=255)
                                                                  aGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTAGCTTATAACAAAcgc                                                                                                                                                                               >  2:376781/1‑251 (MQ=255)
                                                                           ccACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTa                                                                                                                                                                      <  2:154023/251‑1 (MQ=255)
                                                                             acacTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTAc                                                                                                                                                                     >  1:95130/1‑250 (MQ=255)
                                                                              cacTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACaa                                                                                                                                                                   <  2:237630/251‑1 (MQ=255)
                                                                                   cccGATGACAGTCAGCATGACTGACAATCGGTGCGTTTTCTCTATAATAAATCGCTCAATGAATCAGGCATGAATATTCAGTTAGTGGGTTACCGTTAGTCGACCAGCGGATATTTTAAGTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACATCCTCGCTTATAACAACCGCGGGAAATTACATCTCAc                                                                                                                                                              <  2:498465/251‑1 (MQ=255)
                                                                                   cccGATGACAGTCAGCATGACGGTCAAGCGGTGCGTTGTCTCCATAACAAATCGCTCAATGAATCAGGCCCGAATATTCTGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACTGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCAc                                                                                                                                                              <  1:413014/251‑1 (MQ=255)
                                                                                    ccGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATAGTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACAATG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCAca                                                                                                                                                             <  1:70533/251‑2 (MQ=255)
                                                                                    ccGATGACAGTCAGCATAACGGACAATCGGCGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTAAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACTTTT‑AAACACAGGACGGAGTTATTGATGTTAAGCCGAAATTCACCGACTATTAACACATCGATTATAACCAAcgtggaacattaccaaccacc                                                                                                                                                             >  2:238304/1‑232 (MQ=255)
                                                                                     cGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACCTTG‑AACCACAGGACGGAGTTATTCAGGTTAAGCCGCAATTCACCGACACTTACAACCTCGCTTATAACAAACGCGGGAAATTACACCTCACCg                                                                                                                                                            >  2:458494/1‑251 (MQ=255)
                                                                                      gtagACAGTTAGCATGACGGACAAACGGTGCGTTTTTTCCATAACAAATCGCTTAATGAATCAGGCACGAATTTTCAGTTAGTGGGTTACCGGTATTCGACCAGCGGAGATTTTAATTTCGCTGATACAACATACAGTCGAATGAATTGCTACAACATTG‑AAACACAGGACGTAGTTATTCAGGTTAAGCCGAAATGCACCGACTATTACAACCTCGCTTATAACAACCGCGGGAAAttcctacgcaccg                                                                                                                                                            <  2:812393/247‑12 (MQ=255)
                                                                                      gATGACAGTCAGCAAGAAGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCACTGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATAGCTACAACATTGCAAAC‑CAGGACGGATTTCTTCAGGTTAAGCCGAACTTCCCCGAATATTACCACCTCGCTTCTCACAAACGCGGGACATTACCACtaacatt                                                                                                                                                           >  2:810170/1‑245 (MQ=255)
                                                                                                     gACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCg                                                                                                                                             <  1:478632/250‑1 (MQ=255)
                                                                                                       cGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGgc                                                                                                                                          <  1:862928/251‑1 (MQ=255)
                                                                                                         gACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGAACGGAGTTATTAAGGTTAAGCCGAAATTCAACGACTATTACAACCTCGCTTATAAAAAACGCGGGAAAATAAAACTCACCGTTACTAAGCaaataggacgc                                                                                                                                        >  2:495050/1‑242 (MQ=255)
                                                                                                                ggTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCa                                                                                                                                  >  1:708344/1‑250 (MQ=255)
                                                                                                                   gCGTTTTCTCTATAACAAATCGCCCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATTCAGTCGAATTAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGTAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAacac                                                                                                                              <  2:91342/251‑1 (MQ=255)
                                                                                                                    cGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCCGCAACTCGGGCGCACATCAACACt                                                                                                                             >  2:119954/1‑251 (MQ=255)
                                                                                                                         tctcTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTAtt                                                                                                                        >  1:589880/1‑251 (MQ=255)
                                                                                                                          ctctATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTAttt                                                                                                                       >  1:461609/1‑251 (MQ=255)
                                                                                                                              ataACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGt                                                                                                                   >  2:492223/1‑251 (MQ=255)
                                                                                                                              ataAAAAATCGCTCAATGAATCAGGCCCGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTCCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCCCCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCTCCGTTGCTCAGCACCTCGGGCGCACATCAACACTGTATTTGAg                                                                                                                    <  2:445098/250‑1 (MQ=255)
                                                                                                                                    aaTCGCTCAATGAATCAGGCACGAATATTCCGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGACATTACAACTCCCCGTTCCTCAGCAACTCTGTCGCACATCAACACTGTATTTGAGTGGTacc                                                                                                             >  2:833813/1‑249 (MQ=255)
                                                                                                                                       cGCTCAATGAATCAGGCCCGAATATTCAGTTAGTGGGTTACCGTTATTCGCCCAGCGGATATTTTAATTTCGATGATACAACATACAGTCAAATGAATGGCTACAACATTG‑AACCACATAACGGAGTTATTCAGATTAAACCGAAATTCACCGACTCTTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCTTTACTCAGCAACTCGGGCGCACATCAACCCTGTATTTGCGTGGTACCCAt                                                                                                          >  1:382767/1‑251 (MQ=255)
                                                                                                                                           cAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCaaa                                                                                                      <  1:815866/251‑1 (MQ=255)
                                                                                                                                           cAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCaaa                                                                                                      <  1:432165/251‑1 (MQ=255)
                                                                                                                                                  tCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTg                                                                                               >  1:561944/1‑251 (MQ=255)
                                                                                                                                                          taaTATTCAGATAGTGGGTTACAGTTATTCGACCCACGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGGATGGCTAAAACATTG‑CAACACAGGAAGGCATTAATCAGGTTAAGCCGAAATTCACAGACTATTACAACCTCGCTTATAACCAAAGCGGGAAATTACAACCCAACGTACCACCGCCACTCGGGCGCACATCAACCCTGTATTTGCGGGGTAGCCATCACACTTAtgtgggaccc                                                                                        >  2:92556/2‑241 (MQ=255)
                                                                                                                                                          gAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAAcg                                                                                                                                                                                <  1:566272/162‑1 (MQ=255)
                                                                                                                                                          gAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAAcg                                                                                                                                                                                >  2:566272/1‑162 (MQ=255)
                                                                                                                                                           aaTATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAg                                                                                      <  1:502508/251‑1 (MQ=255)
                                                                                                                                                           aaTATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTACGCTGATGCAACAAACAGTCTAATGAAGTGCGACAACATTG‑AAACACAGTACGGAGTTATACAGGTTAAGCCGAAATTCACCGATTATTACAACCTCGCGTATAACAATCGCGGGAAATTACAACTGACCGTTACTCAGCATCTCGGTCGCACGTCAACACTGTATTTGAGTGGTAGTCATCAAACTTATTGGGGAACGAg                                                                                      <  1:49782/251‑1 (MQ=255)
                                                                                                                                                            atatTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGt                                                                                     >  1:680585/1‑251 (MQ=255)
                                                                                                                                                               ttCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGAGCGCACATCAACACTGGATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGGAAt                                                                                  >  1:35161/1‑251 (MQ=255)
                                                                                                                                                                tCAGTTTGTGGTTTACCGTTATTTGCTCAGCGGATATTTGAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAAt                                                                                  <  2:464903/250‑1 (MQ=255)
                                                                                                                                                                    ttAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGAGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGTGGAACGAGTAATGTCg                                                                              >  1:189817/1‑250 (MQ=255)
                                                                                                                                                                      aGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATg                                                                           <  2:273909/251‑1 (MQ=255)
                                                                                                                                                                       gTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCTACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGTGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGGAATGTCGATGa                                                                          <  2:561420/251‑1 (MQ=255)
                                                                                                                                                                           gTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCaa                                                                      <  2:593316/251‑1 (MQ=255)
                                                                                                                                                                           gTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCaa                                                                      <  1:629731/251‑1 (MQ=255)
                                                                                                                                                                                      tCGACCAGCGGATATGTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGTCTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTg                                                            <  1:827489/250‑1 (MQ=255)
                                                                                                                                                                                       cGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGa                                                          <  2:499119/251‑1 (MQ=255)
                                                                                                                                                                                       cGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCGTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTAGTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGa                                                          <  2:223581/251‑1 (MQ=255)
                                                                                                                                                                                        gACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGAt                                                         <  1:54826/251‑1 (MQ=255)
                                                                                                                                                                                            aGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTaaa                                                     <  2:442944/251‑1 (MQ=255)
                                                                                                                                                                                                gATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTGACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACt                                                 <  2:299163/251‑1 (MQ=255)
                                                                                                                                                                                                 atatTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTg                                                >  2:303453/1‑251 (MQ=255)
                                                                                                                                                                                                     tttAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGtt                                            >  1:324977/1‑251 (MQ=255)
                                                                                                                                                                                                       tAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTc                                           <  1:441439/250‑1 (MQ=255)
                                                                                                                                                                                                         aTTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGa                                         >  1:770979/1‑250 (MQ=255)
                                                                                                                                                                                                                 gagaaaaaaTACAGTTGAAGGAATGGCTACAACATTG‑AAACACAGGACGGAGTTAATCAGGTTAAGCCGAAATTCACCGTCTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACGGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGatagca                                 <  2:206626/242‑4 (MQ=255)
                                                                                                                                                                                                                      aacaTACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTgg                            >  1:509120/1‑250 (MQ=255)
                                                                                                                                                                                                                       acaTACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGAc                          <  2:472989/251‑1 (MQ=255)
                                                                                                                                                                                                                       acaTACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGAc                          <  1:320369/251‑1 (MQ=255)
                                                                                                                                                                                                                        cattcAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACg                         <  1:416940/247‑1 (MQ=255)
                                                                                                                                                                                                                         aTACAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGc                        >  1:613640/1‑251 (MQ=255)
                                                                                                                                                                                                                           aCAGTCGAATGAATGGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCt                       <  1:119954/250‑1 (MQ=255)
                                                                                                                                                                                                                            cAGTCGAATGAATGGCTACAACATTG‑ACACACAGGCCGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTCTTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTc                      >  2:443862/1‑250 (MQ=255)
                                                                                                                                                                                                                             aGTCGAATGAATGGCTACAACATCG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCTGGAAATTACAACTCACCGTTACTCAGCAACt                                                                                                                                               <  2:283961/128‑1 (MQ=255)
                                                                                                                                                                                                                             aGTCGAATGAATGGCTACAACATCG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCTGGAAATTACAACTCACCGTTACTCAGCAACt                                                                                                                                               >  1:283961/1‑128 (MQ=255)
                                                                                                                                                                                                                                    tgaatgGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTCTTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGCGTAATGTCGGTGAGACATTCCAGGCTGGATTAAATACTGCGTACGAAGATATCAACTGTACCCTCCGCTATAGc             >  2:333833/1‑251 (MQ=255)
                                                                                                                                                                                                                                     gaatgTCTCCAACATTG‑AAACACAGGTCAGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCAGGAAATTACAACCCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCtatatc             <  1:419449/250‑3 (MQ=255)
                                                                                                                                                                                                                                       atgGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTg          <  1:239794/251‑1 (MQ=255)
                                                                                                                                                                                                                                        tgGCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAACTAATGCGTTCGAAGATATCAACTGGCCGCTCAGCTATAGCCTg          >  1:312000/1‑250 (MQ=255)
                                                                                                                                                                                                                                         ggCTACTACATGG‑AAACACAGGTCGGAGTTATTCAGGGTAAGCCGAAATTCTCCGACTATTACAACCTCTCTTATAGCAAACTCGTGAAATTACTACGCTCCGTTACTCAGCAACTCGGGCGCACATCAACACGATATTTGAGTGGTAGCCATCAAACTTATGGGTGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCTAAGATATCAACTGGACGCCCAGCTATAGCCtta         <  2:346046/250‑3 (MQ=255)
                                                                                                                                                                                                                                         ggCTACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATACCAACAGGACGCTCAGCTATAgaatgac        >  2:580613/1‑245 (MQ=255)
                                                                                                                                                                                                                                            tACAACATTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGa      >  1:497816/1‑250 (MQ=255)
                                                                                                                                                                                                                                             acaacaTTG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGaa     <  2:497816/250‑1 (MQ=255)
                                                                                                                                                                                                                                              caacaTTG‑AAACACAGGCCGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGaaa    >  2:17879/1‑250 (MQ=255)
                                                                                                                                                                                                                                               aacaTTG‑AAACACAGGACGGAGTTATTCAGGTTAATCCGAAATTCACCGACGATAACAGTCCCGCATAGAACAAACGCTGGAAATGACAACTCACGGCTTCTCAGCAACTCTGGCGCTCATCAGCACTGTATTTGAGCGGTATCGAGCAAACTGATTGGGGAACGAGTAAGGTCGATGAGCAAATCCAGGCTGGAGTAAATACTGCGTTCGAAGATAGCAACTGGACGCTCAGCTATGGCCTGACGaaaaa  <  2:867791/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
GTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCG‑AAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAA  >  NC_000913/4546357‑4546846

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