Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,549,932 T→A 9.0% intergenic (+222/+21) fimH → / ← gntP minor component of type 1 fimbriae/fructuronate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,549,9320TA 9.0% 666.2 / 43.7 199intergenic (+222/+21)fimH/gntPminor component of type 1 fimbriae/fructuronate transporter
Reads supporting (aligned to +/- strand):  ref base T (92/89);  new base A (9/9);  total (101/98)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.41e-03

Only 100 of 202 total aligned reads displayed.

ACTGCAGGGAATGTGCAATCGATTATTGGCGTGACTTTTGTTTATCAATAAAGAAATCACAGGACATTGCTAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCAGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGG  >  NC_000913/4549660‑4550203
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                               
aCTGCAGGGAATGTGCAATCGATTATTGGCGTGACTTTTGTTTATCAATAAAGAAATCACAGGACATTGCTAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCAtt                                                                                                                                                                                                                                                                                >  2:434656/1‑274 (MQ=255)
 cTGCAGGGAATGTGCAATCGATTATTGGCGTGACTTTTGTTTATCAATAAAGAAATCACAGGACATTGCTAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCCACCCGGAGGGGCTAGGGGACTGTCCGGGCACATTTTTAGACTTCGTCATAAGTCTGAGCCTGCCAtt                                                                                                                                                                                                                                                                                >  2:122418/1‑273 (MQ=255)
         aaTGTGCAATCGATTATTGGCGTGACTTTTGTTTATCAATAAAGAAATCACAGGACATTGCTAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGTATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGCCTTTGTCATCAGTCTGAGCCTGCAATTGGCAGGccc                                                                                                                                                                                                                                                                       >  4:349198/1‑272 (MQ=255)
         aaTGTGCAATCGATTATTGGCGTGACTTTTGTTTATCAATAAAGAAATCACAGGACATTGCTAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGAGTCAGCGTCGCATATGGCATAAACAAAGCGCCCTTTGCTGGTCTGTTCCCCTCACCCTACCCCTACCCCCGGTGGGGCGCTGGGACTGTCCGGGCACATTTTTAGCCTTTGTCCTAAGCCTGACCGTGCATTTGGCAGGctc                                                                                                                                                                                                                                                                       >  4:127963/1‑274 (MQ=255)
                       tattGGCGTGACTTTTGTTTATCAATAAAGAAATCACAGGACATTGCTAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCGTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACg                                                                                                                                                                                                                                                        <  4:232530/275‑1 (MQ=255)
                              gTGACTTTTGTTTATCAATAAAGAAATCCCAGGACATTGCTAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCa                                                                                                                                                                                                                                                  <  2:699475/274‑1 (MQ=255)
                                     ttGTTTATCAATAAAGAAATCAAAGGACATTGCTAATGCTGGTAGGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCa                                                                                                                                                                                                                                                  <  2:360243/267‑1 (MQ=255)
                                                         cacaGGACATTGCTAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGccc                                                                                                                                                                                                                      <  2:319820/275‑1 (MQ=255)
                                                         cacaGGACATTGCTAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTACGCCAGATAAGCCGCGTCAGCGTCGCATCTCGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGccc                                                                                                                                                                                                                      <  1:90035/275‑1 (MQ=255)
                                                                         tGCTGGGTCGAAAAATAACCTGAAGCTAAAAACCTGGACGTTAGCCCGGTGGAGGCCAGGTAAGACGCGTCAGCGGCGCATCTGGAAGAAACAAAGAGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGAcc                                                                                                                                                                                                      <  4:68619/275‑1 (MQ=255)
                                                                               tACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTcc                                                                                                                                                                                                >  2:306257/1‑275 (MQ=255)
                                                                               tACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTCGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAACCAGCTAAACCACGGAGTCGACCAGCTc                                                                                                                                                                                                 >  2:167453/1‑274 (MQ=255)
                                                                                          aCCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGAcccc                                                                                                                                                                                      <  4:306256/274‑1 (MQ=255)
                                                                                               aaGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGc                                                                                                                                                                                                                                                                         >  2:750892/1‑186 (MQ=255)
                                                                                                aGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCCTCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCCGTCTGAGCCTGCCATTGCCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATAAGCCCAACCCCGGCGTTGACCAGCTCCAGCAGACCCCAGGtttt                                                                                                                                                                               >  2:697396/1‑275 (MQ=255)
                                                                                                    aaaaaCCTGCACGTTCGCCCCTTGTAGGCCAGATACTACGCTTCAGCGTCGCAACTGGCATAAACAAAGCGGACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAAc                                                                                                                                                                           <  4:434654/275‑1 (MQ=255)
                                                                                                    aaaaaCCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAAc                                                                                                                                                                           <  2:349200/275‑1 (MQ=255)
                                                                                                    aaaaaCCGGCAAGTTAGACGTGTGGAGGCCGGATAAGTCGCGTCAGCGGCGCAGCTGGCATAAAAAAATCGCAGTTTGATGGTCTGTTACCCTCGCACTAACCCTCTCCCCGGAGGGGCGGGGGGACTGTCCGGGCACAGTGTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTGTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAAc                                                                                                                                                                           <  4:359261/275‑1 (MQ=255)
                                                                                                                          tGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGctct                                                                                                                                                                                                                                                                      >  2:46792/1‑162 (MQ=255)
                                                                                                                               gCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTgg                                                                                                                                                                                                                                                                              >  7:747967/1‑149 (MQ=255)
                                                                                                                                 cAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCa                                                                                                                                                                                                                                                                            >  8:327639/1‑149 (MQ=255)
                                                                                                                                   gATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAgg                                                                                                                                                                                                                                                                          >  7:183634/1‑149 (MQ=255)
                                                                                                                                    aTAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAg                                                                                                                                                                                                            >  2:48325/1‑210 (MQ=255)
                                                                                                                                       aGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCAATGGCAGGctc                                                                                                                                                                                                                                                                       <  6:700897/148‑1 (MQ=255)
                                                                                                                                          cgcgTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGgt                                                                                                                                                                                                                                                                   <  8:747967/149‑1 (MQ=255)
                                                                                                                                              tCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGAGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTg                                                                                                                                                                                                                                                                               <  8:202223/133‑1 (MQ=255)
                                                                                                                                               cAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCt                                                                                                                                                                                                                                                              <  8:222566/149‑1 (MQ=255)
                                                                                                                                                  cGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGgt                                                                                                                                                                                                                                                                   <  4:809698/141‑1 (MQ=255)
                                                                                                                                                   gTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTa                                                                                                                                                                                                                                                          <  7:327639/149‑1 (MQ=255)
                                                                                                                                                     gtcagtgggCATTAACAAGACGTCCTGTGCGGGTCTGTGCCCCTCACCCTTTCCCGCTCCCCGGAGGGTGGAGGGGGCTGTTCTGGCACATTTTTAGAGTTTGTCATCAGTCGGAGCTTGCCATTGGGAGGCTCTGGGGGCCTTTAACGCTACCATGCTAATAATCAGAACAATAATCAGACCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGagg                                                                                                                          <  4:697394/268‑3 (MQ=255)
                                                                                                                                                               accagcAAAGTGCGCTTTGTTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCt                                                                                                                                                                                                                                              >  8:668342/6‑149 (MQ=38)
                                                                                                                                                               aTAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCt                                                                                                                                                                                                                                              >  7:304418/1‑149 (MQ=255)
                                                                                                                                                               aTAAACAAAGCACACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCt                                                                                                                                                                                                                                              >  7:514057/1‑149 (MQ=255)
                                                                                                                                                                 cagcAAAGTGCGCTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTAGTCATCAGTCTGAGCCTGCCATAGGCAGGATCTCGTGTCCTTTTACGCTACCATGCtaa                                                                                                                                                                                                                                            >  8:560295/4‑149 (MQ=38)
                                                                                                                                                                             cTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGctct                                                                                                                                                                                                                                                                      >  5:880574/1‑111 (MQ=255)
                                                                                                                                                                              tttGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGctc                                                                                                                                                                                                                                                                       >  7:691158/1‑109 (MQ=255)
                                                                                                                                                                                   tGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCa                                                                                                                                                                                                                          <  7:847410/149‑1 (MQ=255)
                                                                                                                                                                                       cTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAgg                                                                                                                                                                                                                                                                          <  2:763734/97‑1 (MQ=255)
                                                                                                                                                                                          ttCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCCGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCa                                                                                       >  2:594022/1‑273 (MQ=255)
                                                                                                                                                                                           tCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAAc                                                                                                                                                                                                                   >  5:856607/1‑148 (MQ=255)
                                                                                                                                                                                              ccTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACg                                                                                                                                                                                                               >  7:549292/1‑149 (MQ=255)
                                                                                                                                                                                               cTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCTTCAGTCTGAGCCTGCCATTGGCAGGCCCTGGTGTCCTTTTACGCTACCATGCTAANAATCAGCNCAANAATCAGNCCAACCACGGAGtt                                                                                                                                                                                                          >  3:34625/1‑153 (MQ=255)
                                                                                                                                                                                               cTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACg                                                                                                                                                                                                               >  6:748435/1‑148 (MQ=255)
                                                                                                                                                                                                tCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGa                                                                                                                                                                                                             >  8:597060/1‑149 (MQ=255)
                                                                                                                                                                                                tCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGa                                                                                                                                                                                                             <  8:304418/149‑1 (MQ=255)
                                                                                                                                                                                                   cccTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCa                                                                                                                                                                                                                                      >  7:420681/1‑121 (MQ=255)
                                                                                                                                                                                                   cccTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCa                                                                                                                                                                                                                                      <  8:420681/121‑1 (MQ=255)
                                                                                                                                                                                                   cccTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGAcc                                                                                                                                                                                                      <  1:34626/153‑1 (MQ=255)
                                                                                                                                                                                                     cTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTg                                                                                                                                                                                                         >  5:377014/1‑148 (MQ=255)
                                                                                                                                                                                                       aaCCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGAc                                                                                                                                                                                                       >  5:523955/1‑148 (MQ=255)
                                                                                                                                                                                                           ctctCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGc                                                                                                                                                                                                   >  5:350310/1‑148 (MQ=255)
                                                                                                                                                                                                           ctctCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGc                                                                                                                                                                                                   >  6:779558/1‑148 (MQ=255)
                                                                                                                                                                                                           ctctCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCt                                                                                                                                                                                                  >  7:283/1‑149 (MQ=255)
                                                                                                                                                                                                             ctcCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCAATGGCAGGct                                                                                                                                                                                                                                                                        <  6:349589/77‑1 (MQ=255)
                                                                                                                                                                                                              tcCCCGGAGGGGCGAGGGGACTGTCCGGTCACACTTCTAGACTTTGTCATCAGTCTTAGCCTCCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGGTGACCAGCTCCAGCAGACCCCAGGTTTGCAACGTGTCTTTTACTGACAGGTCAAAGTAACCCATGAGCAACCAGAAAGATGCATCGCTAATGTGGGTGAGGGTGTTGGAACCGGCAACCGGCGCCAGtccaagatgcg                                                                 >  4:469127/1‑265 (MQ=255)
                                                                                                                                                                                                              tcCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACACTAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCNCGGTTTTCAACGTGTCTTTTACTTCCAGGTCa                                                                                                                                                      >  4:760844/1‑190 (MQ=255)
                                                                                                                                                                                                              tcCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGaaa                                                                                                                              >  4:576930/1‑214 (MQ=255)
                                                                                                                                                                                                                  cGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCa                                                                                                                                                                                                                                      <  2:563586/106‑1 (MQ=255)
                                                                                                                                                                                                                  cGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCcagca                                                                                                                                                                                            >  5:284780/1‑148 (MQ=255)
                                                                                                                                                                                                                  cGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGAc                                                                                                                                                            <  2:808909/180‑1 (MQ=255)
                                                                                                                                                                                                                  cGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTGACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGa                                                                                                                                                             <  4:647301/179‑1 (MQ=255)
                                                                                                                                                                                                                     aGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCAt                                                                                                                                                                                                                                                                                 >  2:90016/1‑60 (MQ=255)
                                                                                                                                                                                                                          cGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCCGCAGCGCCGGATTCACGc                                                     >  2:397276/1‑275 (MQ=255)
                                                                                                                                                                                                                            aGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCCACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGGCGCCAGTACCCGCCACCCCGGATTCCCGcc                                                    >  4:526522/1‑274 (MQ=255)
                                                                                                                                                                                                                                 aCTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCaa                                                                                                                                                                            <  8:514057/149‑1 (MQ=255)
                                                                                                                                                                                                                                 aCTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCaa                                                                                                                                                                            <  8:132417/149‑1 (MQ=255)
                                                                                                                                                                                                                                    gTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACgt                                                                                                                                                                         <  8:126400/149‑1 (MQ=255)
                                                                                                                                                                                                                                      ccGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACgtg                                                                                                                                                                        >  5:712296/1‑148 (MQ=255)
                                                                                                                                                                                                                                         ggCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGTTAATAATTAGCAAAATAAGTAGCCCAACCACGGACTTGACTACCTCCAGCAGACTCCAGGTGTTCAACGTGTCTTTTACGGACAGGTCAAAGTACCCTTTGGACAACCAGAAAGGTGCAGCGTCAATGTGGGAGAGGGTGTTGGAACCGGCAGCCGGCGCCAGGAACAGCAGGCCAGGATTAACACCAACCAGCTGACAAg                                      >  4:473628/1‑275 (MQ=255)
                                                                                                                                                                                                                                          gCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCt                                                                                                                                                                                                                                                      >  4:382000/1‑66 (MQ=255)
                                                                                                                                                                                                                                               tttttAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTg                                                                                                                                                                                                         >  6:83964/1‑106 (MQ=255)
                                                                                                                                                                                                                                               tttttAGACTTGTGCATCAGTCGTAGCATGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCAGGCTAATAATCAGCACAAGAATCAGACCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGAAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCAGTTGCt                                 <  2:198055/274‑1 (MQ=255)
                                                                                                                                                                                                                                                  ttAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGAc                                                                                                                                                            <  6:523955/148‑1 (MQ=255)
                                                                                                                                                                                                                                                      aCTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTATCAACGTGTCTTTTACTGACAGGTc                                                                                                                                                       >  7:589940/1‑149 (MQ=255)
                                                                                                                                                                                                                                                       cTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACt                                                                                                                                                               >  8:701795/1‑140 (MQ=255)
                                                                                                                                                                                                                                                       cTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCa                                                                                                                                                      >  7:215177/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         ttGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCtaataa                                                                                                                                                                                                                                         >  6:217753/1‑64 (MQ=255)
                                                                                                                                                                                                                                                         ttGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCaaa                                                                                                                                                    <  7:560295/149‑1 (MQ=255)
                                                                                                                                                                                                                                                         ttGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCaa                                                                                                                                                     >  5:840342/1‑148 (MQ=255)
                                                                                                                                                                                                                                                          tGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGAc                                                                                                                                                            >  3:75691/1‑140 (MQ=255)
                                                                                                                                                                                                                                                          tGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGaaa                                                                                                                              >  4:458311/1‑170 (MQ=255)
                                                                                                                                                                                                                                                            tcatcaGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAg                                                                                                                                                                                                    <  2:284111/98‑1 (MQ=255)
                                                                                                                                                                                                                                                              atcaGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCa                                                                                                                                                                                                                                                  <  4:776359/50‑1 (MQ=255)
                                                                                                                                                                                                                                                              atcaGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGccc                                                                                                                                                                                                                      >  7:128611/1‑78 (MQ=255)
                                                                                                                                                                                                                                                              atcaGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGATCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCAGTTGCt                                 <  2:526663/259‑1 (MQ=255)
                                                                                                                                                                                                                                                                   tCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCttt                                                                                                                                           <  6:428259/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                    tagAGCCTGCCAATGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTa                                                                                                                                                 >  6:725639/3‑141 (MQ=255)
                                                                                                                                                                                                                                                                    cagAGCCTGCCAATGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTa                                                                                                                                                 <  5:725639/139‑1 (MQ=255)
                                                                                                                                                                                                                                                                    cTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGa                                                                                                                                         <  7:6474/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                     agAGCCTGCCAATGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCaa                                                                                                                                                     >  6:846610/2‑136 (MQ=255)
                                                                                                                                                                                                                                                                      gAGCCTGCCAATGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCaa                                                                                                                                                                            <  6:642271/112‑1 (MQ=255)
                                                                                                                                                                                                                                                                       aGCCTGCCAATGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTg                                                                                                                                                                                                         <  8:380391/82‑1 (MQ=255)
                                                                                                                                                                                                                                                                       aGCCTGCCAATGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGaaca                                                                                                                                      >  7:361492/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                        gCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGaaca                                                                                                                                      <  6:856607/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                        gCCTGCCAATGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCaa                                                                                                                                                     <  5:112179/133‑1 (MQ=255)
                                                                                                                                                                                                                                                                          cTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGa                                                                                                                            >  4:530748/1‑156 (MQ=255)
                                                                                                                                                                                                                                                                            gCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCa                                                                                                                                  >  5:636095/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                             ccATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCAGTTGCTGGATCTAGGATAGCAGCACTGATAATACcgg  >  2:734498/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                             ccATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAAAGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATTCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCAAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCAGTTGCTGGAGCTAGGATTGCAGAAACGATAATCCcgg  >  1:159771/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                              cATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGa                                                                                                                                <  5:77834/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                               
ACTGCAGGGAATGTGCAATCGATTATTGGCGTGACTTTTGTTTATCAATAAAGAAATCACAGGACATTGCTAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCAGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGG  >  NC_000913/4549660‑4550203

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