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→ATTfimDtype 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

GTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAA  >  NC_000913/4546357‑4546846
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
gTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAacc                                                                                                                                                                                                                                                 >  1:816499/1‑250 (MQ=255)
  gCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTCCAACATTGAAacacc                                                                                                                                                                                                                                               >  2:273411/1‑250 (MQ=255)
  gCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAacac                                                                                                                                                                                                                                                <  1:17863/250‑1 (MQ=255)
                          aaCATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGAACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGcc                                                                                                                                                                                                                       >  1:472989/1‑251 (MQ=255)
                              tGGGGACACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCTAATGAATGCCCACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGaa                                                                                                                                                                                                                    >  2:836753/1‑250 (MQ=255)
                                  ggCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTc                                                                                                                                                                                                                >  1:299163/1‑250 (MQ=255)
                                    cACTTGGCGCTCTGTCTGTGGATATGATGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTTCGTTTTCTCTATAACAAATCGCTCAATGAATCAGTCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACAGTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCAcc                                                                                                                                                                                                             <  2:578588/251‑1 (MQ=255)
                                    cACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTGATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCTCGAATATTCAGTGAGTGGGTTACCGTTAGTCGACGAGCTGATATTTTAATTTCGCTGATACAACATACAGTCGGATGAATGGCTACAACAGGGAAACACAGGACGGAGTTATTCAGGTTAATCCGAAATTCAcc                                                                                                                                                                                                             <  2:331867/251‑1 (MQ=255)
                                     ccTGGGTGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCg                                                                                                                                                                                                            <  2:590615/250‑1 (MQ=255)
                                     aCTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGACACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCg                                                                                                                                                                                                            >  2:478632/1‑251 (MQ=255)
                                      cTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGACAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATCTTCAGTTAGTGGGTTACCGTTATTCGACCAACGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAAAATTAAAACACAGGCCGGAGTTATTCAGGTAAAGCCGAAATTCACCg                                                                                                                                                                                                            >  2:307548/1‑250 (MQ=255)
                                        gggCGCTCTGTCTGGGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGTACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGa                                                                                                                                                                                                           <  2:438978/249‑1 (MQ=255)
                                         ggCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATACCAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACt                                                                                                                                                                                                         >  1:854741/1‑250 (MQ=255)
                                         ggCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACt                                                                                                                                                                                                         >  1:66651/1‑250 (MQ=255)
                                           cgcTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGCCGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTAtt                                                                                                                                                                                                      >  2:370876/1‑251 (MQ=255)
                                           cgcTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTAtt                                                                                                                                                                                                      <  1:38567/251‑1 (MQ=255)
                                            gcTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTATATAACAACTCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAac                                                                                                                                                                                                                                                  >  1:388623/1‑206 (MQ=255)
                                            gcTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGGGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGGTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGCCGAATGAATGGCTACAACATTGAAac                                                                                                                                                                                                                                                  <  2:388623/206‑1 (MQ=255)
                                                 gtctgtGGATATGACGCAGGCTAATTCCACACTTACCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACCTACAGTCGAATGAATGGCTACAACATTGAAACACATGACGGAGTTATTCAGGTTAAGCCGAACTTCACCGACTATAACAAcc                                                                                                                                                                                                >  1:194076/1‑251 (MQ=255)
                                                     gtgGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTc                                                                                                                                                                                              <  1:360283/249‑1 (MQ=255)
                                                     gtgGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGc                                                                                                                                                                                            >  1:232535/1‑251 (MQ=255)
                                                             gACGCAGGCTAATTCCACAATTCCCGATGACAGTCAGCATTACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATTGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCCCCGACTATTACAACCTCGCTTATAACa                                                                                                                                                                                    <  2:66108/251‑1 (MQ=255)
                                                                  aGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTAGCTTATAACAAAcgc                                                                                                                                                                               >  2:376781/1‑251 (MQ=255)
                                                                           ccACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTa                                                                                                                                                                      <  2:154023/251‑1 (MQ=255)
                                                                             acacTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTAc                                                                                                                                                                     >  1:95130/1‑250 (MQ=255)
                                                                              cacTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACaa                                                                                                                                                                   <  2:237630/251‑1 (MQ=255)
                                                                                   cccGATGACAGTCAGCATGACTGACAATCGGTGCGTTTTCTCTATAATAAATCGCTCAATGAATCAGGCATGAATATTCAGTTAGTGGGTTACCGTTAGTCGACCAGCGGATATTTTAAGTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACATCCTCGCTTATAACAACCGCGGGAAATTACATCTCAc                                                                                                                                                              <  2:498465/251‑1 (MQ=255)
                                                                                   cccGATGACAGTCAGCATGACGGTCAAGCGGTGCGTTGTCTCCATAACAAATCGCTCAATGAATCAGGCCCGAATATTCTGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACTGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCAc                                                                                                                                                              <  1:413014/251‑1 (MQ=255)
                                                                                    ccGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATAGTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACAATGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCAca                                                                                                                                                             <  1:70533/251‑2 (MQ=255)
                                                                                    ccGATGACAGTCAGCATAACGGACAATCGGCGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTAAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACTTTTAAACACAGGACGGAGTTATTGATGTTAAGCCGAAATTCACCGACTATTAACACATCGATTATAACCAAcgtggaacattaccaaccacc                                                                                                                                                             >  2:238304/1‑232 (MQ=255)
                                                                                     cGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACCTTGAACCACAGGACGGAGTTATTCAGGTTAAGCCGCAATTCACCGACACTTACAACCTCGCTTATAACAAACGCGGGAAATTACACCTCACCg                                                                                                                                                            >  2:458494/1‑251 (MQ=255)
                                                                                      gtagACAGTTAGCATGACGGACAAACGGTGCGTTTTTTCCATAACAAATCGCTTAATGAATCAGGCACGAATTTTCAGTTAGTGGGTTACCGGTATTCGACCAGCGGAGATTTTAATTTCGCTGATACAACATACAGTCGAATGAATTGCTACAACATTGAAACACAGGACGTAGTTATTCAGGTTAAGCCGAAATGCACCGACTATTACAACCTCGCTTATAACAACCGCGGGAAAttcctacgcaccg                                                                                                                                                            <  2:812393/247‑12 (MQ=255)
                                                                                      gATGACAGTCAGCAAGAAGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCACTGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATAGCTACAACATTGCAAACCAGGACGGATTTCTTCAGGTTAAGCCGAACTTCCCCGAATATTACCACCTCGCTTCTCACAAACGCGGGACATTACCACtaacatt                                                                                                                                                           >  2:810170/1‑245 (MQ=255)
                                                                                                     gACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCg                                                                                                                                             <  1:478632/250‑1 (MQ=255)
                                                                                                       cGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGgc                                                                                                                                          <  1:862928/251‑1 (MQ=255)
                                                                                                         gACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGAACGGAGTTATTAAGGTTAAGCCGAAATTCAACGACTATTACAACCTCGCTTATAAAAAACGCGGGAAAATAAAACTCACCGTTACTAAGCaaataggacgc                                                                                                                                        >  2:495050/1‑242 (MQ=255)
                                                                                                                ggTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCa                                                                                                                                  >  1:708344/1‑250 (MQ=255)
                                                                                                                   gCGTTTTCTCTATAACAAATCGCCCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATTCAGTCGAATTAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGTAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAacac                                                                                                                              <  2:91342/251‑1 (MQ=255)
                                                                                                                    cGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCCGCAACTCGGGCGCACATCAACACt                                                                                                                             >  2:119954/1‑251 (MQ=255)
                                                                                                                         tctcTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTAtt                                                                                                                        >  1:589880/1‑251 (MQ=255)
                                                                                                                          ctctATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTAttt                                                                                                                       >  1:461609/1‑251 (MQ=255)
                                                                                                                              ataACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGt                                                                                                                   >  2:492223/1‑251 (MQ=255)
                                                                                                                              ataAAAAATCGCTCAATGAATCAGGCCCGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTCCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCCCCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCTCCGTTGCTCAGCACCTCGGGCGCACATCAACACTGTATTTGAg                                                                                                                    <  2:445098/250‑1 (MQ=255)
                                                                                                                                    aaTCGCTCAATGAATCAGGCACGAATATTCCGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGACATTACAACTCCCCGTTCCTCAGCAACTCTGTCGCACATCAACACTGTATTTGAGTGGTacc                                                                                                             >  2:833813/1‑249 (MQ=255)
                                                                                                                                       cGCTCAATGAATCAGGCCCGAATATTCAGTTAGTGGGTTACCGTTATTCGCCCAGCGGATATTTTAATTTCGATGATACAACATACAGTCAAATGAATGGCTACAACATTGAACCACATAACGGAGTTATTCAGATTAAACCGAAATTCACCGACTCTTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCTTTACTCAGCAACTCGGGCGCACATCAACCCTGTATTTGCGTGGTACCCAt                                                                                                          >  1:382767/1‑251 (MQ=255)
                                                                                                                                           cAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCaaa                                                                                                      <  1:815866/251‑1 (MQ=255)
                                                                                                                                           cAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCaaa                                                                                                      <  1:432165/251‑1 (MQ=255)
                                                                                                                                                  tCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTg                                                                                               >  1:561944/1‑251 (MQ=255)
                                                                                                                                                          taaTATTCAGATAGTGGGTTACAGTTATTCGACCCACGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGGATGGCTAAAACATTGCAACACAGGAAGGCATTAATCAGGTTAAGCCGAAATTCACAGACTATTACAACCTCGCTTATAACCAAAGCGGGAAATTACAACCCAACGTACCACCGCCACTCGGGCGCACATCAACCCTGTATTTGCGGGGTAGCCATCACACTTAtgtgggaccc                                                                                        >  2:92556/2‑241 (MQ=255)
                                                                                                                                                          gAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAAcg                                                                                                                                                                                <  1:566272/162‑1 (MQ=255)
                                                                                                                                                          gAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAAcg                                                                                                                                                                                >  2:566272/1‑162 (MQ=255)
                                                                                                                                                           aaTATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAg                                                                                      <  1:502508/251‑1 (MQ=255)
                                                                                                                                                           aaTATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTACGCTGATGCAACAAACAGTCTAATGAAGTGCGACAACATTGAAACACAGTACGGAGTTATACAGGTTAAGCCGAAATTCACCGATTATTACAACCTCGCGTATAACAATCGCGGGAAATTACAACTGACCGTTACTCAGCATCTCGGTCGCACGTCAACACTGTATTTGAGTGGTAGTCATCAAACTTATTGGGGAACGAg                                                                                      <  1:49782/251‑1 (MQ=255)
                                                                                                                                                            atatTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGt                                                                                     >  1:680585/1‑251 (MQ=255)
                                                                                                                                                               ttCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGAGCGCACATCAACACTGGATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGGAAt                                                                                  >  1:35161/1‑251 (MQ=255)
                                                                                                                                                                tCAGTTTGTGGTTTACCGTTATTTGCTCAGCGGATATTTGAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAAt                                                                                  <  2:464903/250‑1 (MQ=255)
                                                                                                                                                                    ttAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGAGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGTGGAACGAGTAATGTCg                                                                              >  1:189817/1‑250 (MQ=255)
                                                                                                                                                                      aGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATg                                                                           <  2:273909/251‑1 (MQ=255)
                                                                                                                                                                       gTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCTACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGTGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGGAATGTCGATGa                                                                          <  2:561420/251‑1 (MQ=255)
                                                                                                                                                                           gTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCaa                                                                      <  2:593316/251‑1 (MQ=255)
                                                                                                                                                                           gTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCaa                                                                      <  1:629731/251‑1 (MQ=255)
                                                                                                                                                                                      tCGACCAGCGGATATGTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGTCTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTg                                                            <  1:827489/250‑1 (MQ=255)
                                                                                                                                                                                       cGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGa                                                          <  2:499119/251‑1 (MQ=255)
                                                                                                                                                                                       cGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCGTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTAGTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGa                                                          <  2:223581/251‑1 (MQ=255)
                                                                                                                                                                                        gACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGAt                                                         <  1:54826/251‑1 (MQ=255)
                                                                                                                                                                                            aGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTaaa                                                     <  2:442944/251‑1 (MQ=255)
                                                                                                                                                                                                gATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTGACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACt                                                 <  2:299163/251‑1 (MQ=255)
                                                                                                                                                                                                 atatTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTg                                                >  2:303453/1‑251 (MQ=255)
                                                                                                                                                                                                     tttAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGtt                                            >  1:324977/1‑251 (MQ=255)
                                                                                                                                                                                                       tAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTc                                           <  1:441439/250‑1 (MQ=255)
                                                                                                                                                                                                         aTTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGa                                         >  1:770979/1‑250 (MQ=255)
                                                                                                                                                                                                                 gagaaaaaaTACAGTTGAAGGAATGGCTACAACATTGAAACACAGGACGGAGTTAATCAGGTTAAGCCGAAATTCACCGTCTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACGGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGatagca                                 <  2:206626/242‑4 (MQ=255)
                                                                                                                                                                                                                      aacaTACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTgg                            >  1:509120/1‑250 (MQ=255)
                                                                                                                                                                                                                       acaTACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGAc                          <  2:472989/251‑1 (MQ=255)
                                                                                                                                                                                                                       acaTACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGAc                          <  1:320369/251‑1 (MQ=255)
                                                                                                                                                                                                                        cattcAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACg                         <  1:416940/247‑1 (MQ=255)
                                                                                                                                                                                                                         aTACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGc                        >  1:613640/1‑251 (MQ=255)
                                                                                                                                                                                                                           aCAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCt                       <  1:119954/250‑1 (MQ=255)
                                                                                                                                                                                                                            cAGTCGAATGAATGGCTACAACATTGACACACAGGCCGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTCTTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTc                      >  2:443862/1‑250 (MQ=255)
                                                                                                                                                                                                                             aGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCTGGAAATTACAACTCACCGTTACTCAGCAACt                                                                                                                                               <  2:283961/128‑1 (MQ=255)
                                                                                                                                                                                                                             aGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCTGGAAATTACAACTCACCGTTACTCAGCAACt                                                                                                                                               >  1:283961/1‑128 (MQ=255)
                                                                                                                                                                                                                                    tgaatgGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTCTTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGCGTAATGTCGGTGAGACATTCCAGGCTGGATTAAATACTGCGTACGAAGATATCAACTGTACCCTCCGCTATAGc             >  2:333833/1‑251 (MQ=255)
                                                                                                                                                                                                                                     gaatgTCTCCAACATTGAAACACAGGTCAGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCAGGAAATTACAACCCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCtatatc             <  1:419449/250‑3 (MQ=255)
                                                                                                                                                                                                                                       atgGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTg          <  1:239794/251‑1 (MQ=255)
                                                                                                                                                                                                                                        tgGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAACTAATGCGTTCGAAGATATCAACTGGCCGCTCAGCTATAGCCTg          >  1:312000/1‑250 (MQ=255)
                                                                                                                                                                                                                                         ggCTACTACATGGAAACACAGGTCGGAGTTATTCAGGGTAAGCCGAAATTCTCCGACTATTACAACCTCTCTTATAGCAAACTCGTGAAATTACTACGCTCCGTTACTCAGCAACTCGGGCGCACATCAACACGATATTTGAGTGGTAGCCATCAAACTTATGGGTGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCTAAGATATCAACTGGACGCCCAGCTATAGCCtta         <  2:346046/250‑3 (MQ=255)
                                                                                                                                                                                                                                         ggCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATACCAACAGGACGCTCAGCTATAgaatgac        >  2:580613/1‑245 (MQ=255)
                                                                                                                                                                                                                                            tACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGa      >  1:497816/1‑250 (MQ=255)
                                                                                                                                                                                                                                             acaacaTTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGaa     <  2:497816/250‑1 (MQ=255)
                                                                                                                                                                                                                                              caacaTTGAAACACAGGCCGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGaaa    >  2:17879/1‑250 (MQ=255)
                                                                                                                                                                                                                                               aacaTTGAAACACAGGACGGAGTTATTCAGGTTAATCCGAAATTCACCGACGATAACAGTCCCGCATAGAACAAACGCTGGAAATGACAACTCACGGCTTCTCAGCAACTCTGGCGCTCATCAGCACTGTATTTGAGCGGTATCGAGCAAACTGATTGGGGAACGAGTAAGGTCGATGAGCAAATCCAGGCTGGAGTAAATACTGCGTTCGAAGATAGCAACTGGACGCTCAGCTATGGCCTGACGaaaaa  <  2:867791/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
GTGCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAA  >  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: