Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,875,525 G→T 15.6% V30F (GTT→TTT)  yidX → putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,875,5250GT15.6% 159.0 / 5.8 69V30F (GTT→TTT) yidXputative lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (25/32);  new base T (9/3);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 6.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

ATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTTACTACTACTGGAATAACGAAG  >  NC_000913/3875242‑3875813
                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                
atgaCGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTc                                                                                                                                                                                                                                                                                              >  1:707158/1‑288 (MQ=255)
         gcgTTTGGCGTACATGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGAGCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCa                                                                                                                                                                                                                                                                                    <  2:773180/289‑1 (MQ=255)
              taccGGGCGTGTAATGTAGTGCAAACTTCTAGTGTTGTACTACAATTTAGATCACATAAAGAACAATGAGTAATCAAATGACATGAGTCGAGCAGAAATCTGAAAAGGGATATCAGGCGGTAAATAGAAGGGAAAGAAGTGTATGCTTTCAACGGCTTAGCTAAGCGTTTAAAGGATTAATCATGAAGTTGAATTTTAATGGATAGTTTAAGGCTGACGGTTTATACCCACTGGCGCTGATGCTTTAAGGCTGTATCTCGTATGCTGTGGTTTCCCATAGCGCaggtgg                                                                                                                                                                                                                                                                               <  2:516491/287‑6 (MQ=255)
                                 taCAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTAACaa                                                                                                                                                                                                                                                            >  1:325918/1‑289 (MQ=255)
                                 taCAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTGGTTTGCTAGACAGAAAAGGGTAGTTCAGGAAAGTATCaa                                                                                                                                                                                                                                                            >  1:450636/1‑289 (MQ=255)
                                 taCAAACTTAGAATGTTGTACTACAATGTAGATCACAAAAAGAACAATGCATAAAAAAATTACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATTAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCa                                                                                                                                                                                                                                                                                    <  2:962413/265‑1 (MQ=255)
                                  aCAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCaa                                                                                                                                                                                                                                                            >  1:130659/1‑288 (MQ=255)
                                    aaatctataTTGTTGGACTACAATTTAAATACAAAAAATGACAATGCATACAAAAATGACATGGGTAGAGAAGAAAACTGAAAATGGATATCAGGCGCTAACAAAGAGGAAAAGAAGAGTATGCTTACAACGGCTTAGCTACTCGTTTAAAAGATTAAGCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTACCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCg                                                                                                                                                                                                                                                         <  2:567405/284‑1 (MQ=255)
                                      aCTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTGACAAAAAAGAGAAGTTCAGGACAGTATCAATCGCa                                                                                                                                                                                                                                                       >  1:648837/1‑289 (MQ=255)
                                           tatTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTCGTTTCCCATACCGCAAAGGGTAGtt                                                                                                                                                                                                                                                                          >  1:962413/1‑265 (MQ=255)
                                                ttgtACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTACCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGTTTTCCTATTCCGCAAAGGTTAGTTCAGTAAAGTATCAATCGCAGTCAGacac                                                                                                                                                                                                                                              >  2:795230/1‑288 (MQ=255)
                                                      tacAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCa                                                                                                                                                                                                                                                                                    <  2:613405/244‑1 (MQ=255)
                                                      tacAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCa                                                                                                                                                                                                                                                                                    >  1:613405/1‑244 (MQ=255)
                                                          aTTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTGTTCTGCTATAACGAAAAGAGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTgg                                                                                                                                                                                                                                     >  1:371350/1‑287 (MQ=255)
                                                           tttAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGATGAATTTGAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGAGCTGATGCTTTCAGGCTGTATCTCGTATGCCTTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCACTCGAAGTCGAACACCATCACTGGGCt                                                                                                                                                                                                                                  >  2:886056/1‑289 (MQ=255)
                                                                     aaaaaGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGAAAAGGGTAGGTCAGGAAATTATCAATCGGAGTCAGCCATCACTCTTGGGTTATCGCAGGcc                                                                                                                                                                                                                        >  2:821402/1‑288 (MQ=255)
                                                                        aaGAACAATGCAGAAAAAAAGGACATGCGTCGGGCAGAAAGCTGGAAAGGGATATCAGGCGCTAAACAGGACGGAAAGAATAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATGTTGTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCaaaa                                                                                                                                                                                                                     <  2:707158/289‑1 (MQ=255)
                                                                               aTGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTTTGCCTTACCAAAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTa                                                                                                                                                                                                              >  1:62900/1‑289 (MQ=255)
                                                                               aTGCATAAAAAAATGACATGCGTCGGACAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGATGAGTATGCTTTCAACGGCTTAGCTACTGGTTTAAAGGATTAATAATGAAGTTGAATATTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCACATCCCGCAAAGGGTAGTTCAGGAAAGTATCAAGAGCAGTCCTACAACATCATCTGGCTATCTCAGGCCAAATATAga                                                                                                                                                                                                               >  2:978860/1‑287 (MQ=255)
                                                                                        aaaaTGACATGCGTCGGGAAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGGATGCTTTCAACGGCTTAGCAACTCGTTTAAAGGATTAATCATGAAGTGGAATTTTAAGGGATTTTTGAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACaa                                                                                                                                                                                                     <  2:288244/289‑1 (MQ=255)
                                                                                          aaTGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTGGTTTCCTATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGAAACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACaaaa                                                                                                                                                                                                   >  1:758684/1‑289 (MQ=255)
                                                                                                                  tGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATc                                                                                                                                                                            >  2:329329/1‑288 (MQ=255)
                                                                                                                    aaaaGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCg                                                                                                                                                                           <  1:741319/287‑1 (MQ=255)
                                                                                                                       aGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGa                                                                                                                                                                      >  1:294115/1‑289 (MQ=255)
                                                                                                                        tggtgatCAGGCGCTAAACAGGAGGGAAAGAATTGTATGCTTTCAACGGCTTAGCTACTCGTTTAAATGATTAATCATGAAGTTGAATTTTAAGTGATTTTTTAAGGCTGCATGTTTATTCCCACTGGCGCTGATGCTTTCAAGCTTTATCTCGTATGCTCTGTTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCAACACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGAt                                                                                                                                                                     <  2:758684/284‑1 (MQ=255)
                                                                                                                           atatCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTa                                                                                                                                                                   <  2:765863/288‑1 (MQ=255)
                                                                                                                                   cgcTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATc                                                                                                                                                            >  1:653851/1‑287 (MQ=255)
                                                                                                                                      tAAACAGGGGGGAAAGAAGAGGTTGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAAATTTAAGTGATTTTTTAAGGTTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCATATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGa                                                                                                                                                       <  1:821402/289‑1 (MQ=255)
                                                                                                                                      tAAACAGGGAGGAAATAAGTGTATGCTTTCAACGGCTTAGCTCCTCGTTGAAGGGATTAATAATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCAGTTTATACCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGa                                                                                                                                                       <  1:390606/289‑1 (MQ=255)
                                                                                                                                             gaggGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTCTTCTGCCTTAACGCAAAGAGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGcc                                                                                                                                                >  1:953757/1‑289 (MQ=255)
                                                                                                                                                  aaagaagatttatgtttCAACTGCTTAGCGACTCGTTTAAAGGATTTATCATGAAGTTGAATTTTAATGGATTTTTTAAGTCTTACGGTTTATTCCAAAAGGCGCTGATGCTTTCAGGCTGTATCTCGTAGGCTCTGGTTTACCATACAGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGAGAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGa                                                                                                                                           <  2:526931/275‑1 (MQ=255)
                                                                                                                                                  aaaatagagTATGGTTTCAACGGCTGAGCTACTCTTTTAAAGGATTAATCATGAATTTGAATTTTAAGGGATTTTTATAGTGTACCGTTTTATTCCCACAGGCGCTGATGCGTTCATGCTGTATCGCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGa                                                                                                                                           <  2:450636/284‑1 (MQ=255)
                                                                                                                                                   aatgagggTATGCTTTACACGGCTTAGCTACTCGTTTAAATGATTAATCATGAAGTTGAATTATAAGGGATTTTTTAAGGTTGCCTGTTTATTCCGACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCGTCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGAc                                                                                                                                          <  2:439995/282‑1 (MQ=255)
                                                                                                                                                   aataagaGTATGCTTTCAACTGCTTAGCTATTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGTATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGTCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGa                                                                                                                                           <  2:953757/285‑1 (MQ=255)
                                                                                                                                                   aagaagTGGATGCTTTCAAATGGTTAGCGACTCGTGTAAAGGATTAATCATGAAGTTTAATTTTAAGGGATTTTTTAAGGCTGCCTGTTTATTCCCACTGGAGATGATGCTTTCAGGCTATATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGa                                                                                                                                           <  2:963677/288‑1 (MQ=255)
                                                                                                                                                   aagaagAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAAGTTTAAGGGATTTTTTAAGGTTGCCGGTTTATTCCAACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGa                                                                                                                                           <  1:443859/288‑1 (MQ=255)
                                                                                                                                                           tATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAGGGTAGTTCAGCAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTGTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTaa                                                                                                                                  >  1:632992/1‑289 (MQ=255)
                                                                                                                                                            aTGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGa                                                                                                                                                                                                                   <  1:23802/207‑1 (MQ=255)
                                                                                                                                                            aTGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGa                                                                                                                                                                                                                   >  2:23802/1‑207 (MQ=255)
                                                                                                                                                            aTGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGt                                                                                                                                                   >  1:781853/1‑271 (MQ=255)
                                                                                                                                                            aTGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGt                                                                                                                                                   <  2:781853/271‑1 (MQ=255)
                                                                                                                                                                attaaaCGGCTTAGCTATTCGTTGAAAGGATTGATCGTTGATGTGAATTTTAAGGGATTTTTTTAGGCTTCCTGTTTATTCCAACTGGCGATGATGCTTTAAGGCTGTATCTCGTATGGTGTGGTGTCCCAGACTGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGTTAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATg                                                                                                                              <  2:607476/284‑1 (MQ=255)
                                                                                                                                                                 ttCAACGGCTTATGTACTCGTTTAACGGAGTAATCAGTAAGTTAAATTTTAAGTGATTTTTGAAGACTGCCGGTATATTCCGACTGGCGCTGATGCTTTCAGGCTGTATCACGTATGCTCTGGTTGCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGTGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCt                                                                                                                            <  2:325918/289‑1 (MQ=255)
                                                                                                                                                                   cAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGt                                                                                                                                                                                           >  1:420801/1‑224 (MQ=255)
                                                                                                                                                                   cAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGt                                                                                                                                                                                           <  2:420801/224‑1 (MQ=255)
                                                                                                                                                                                   cGTTTAAAGTATTACTCATGAATTTTAATTTTAAGGGATTGTTTGAGGTTGCTGGTTTATTCCCTCTGGCGCTGATGCTTTCAGACTGTATCTCATATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTTGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTaa                                                                                                          <  1:886056/289‑1 (MQ=255)
                                                                                                                                                                                      ttAAAGAATTAATCATTAAATTGAATTTTAAGGGATTTTTTAGAGCTTCCGGTTTATTCCCACTTGCGCTGATGCTGTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCg                                                                                                       <  2:294115/289‑1 (MQ=255)
                                                                                                                                                                                             aTTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGcaca                                                                                                  <  1:728537/287‑1 (MQ=255)
                                                                                                                                                                                              ttAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAACGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACCTCGTTAAGATGCTCAATGATCCAGCACTTAATCGGCACAata                                                                                               >  2:646190/1‑289 (MQ=255)
                                                                                                                                                                                              ttAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCCTGGTTTGCCATAACGAAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAata                                                                                               >  1:35346/1‑289 (MQ=255)
                                                                                                                                                                                                      gAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCCCTGGGCTATCGCAGGCACAAGATAGTAATGGAACTAAAGGCTATGTTTTTGTAGGGGAATCGCTGGTTTACCTTATCACTGATGGTGCGGAGGAATTCCTTACGCTGCTCAATGATCCCGCCTTAAACCGGACAAATACTCAGGtt                                                                                       >  2:541544/1‑289 (MQ=255)
                                                                                                                                                                                                                  tAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCGTCTCCTATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCa                                                                           >  1:583110/1‑289 (MQ=255)
                                                                                                                                                                                                                    aGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTGTTTTGCCAGAAAGAAAAGAGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGTTGACATCGTTAAGATGATCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGACGATGACGCAAg                                                                         >  1:748463/1‑289 (MQ=255)
                                                                                                                                                                                                                         ttttttAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCTTCTGCCATAAAGAAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTg                                                                    >  1:447712/1‑289 (MQ=255)
                                                                                                                                                                                                                                   cTGCCGGTTTATACCCACTGGCGCTGATTCTTTCAGGCTGTATCTCGTATGCTCTGTTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGTTCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATg                                                           <  1:795230/288‑1 (MQ=255)
                                                                                                                                                                                                                                      ccGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCggg                                                       >  1:787443/1‑289 (MQ=255)
                                                                                                                                                                                                                                      ccGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTGTTTTGCCTGACAGAAAAGGGTAGTTCAGGCAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGAAAAAGATCGTAATGGAACCAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGAACCAGCCCTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCggt                                                       >  1:673082/1‑288 (MQ=255)
                                                                                                                                                                                                                                      ccGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCCGTCGTCTGCTTTACCACAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCggg                                                       >  1:451662/1‑289 (MQ=255)
                                                                                                                                                                                                                                       cGGTTTATGCACACTGGCGCTGATGTTTTGAGGCTGAATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATAAATAGAAGTAAGAAACCATCACTGGGATAAAGCAGGCAAAAGATAGTAATGGAACAAAAGGATATGTTTTTGTAGGTGAATAGCGGGATTACATTATCACTTATGGTGCAGATGACATCGTTAAGATGATCAATGATCCAGAACTTAACCGTTACAATATTCAGGTTGCCGATGGCGCAAGATTTGTTTTAAATGCggg                                                       <  2:62900/288‑1 (MQ=255)
                                                                                                                                                                                                                                                         gcTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCgg                                     <  1:329329/288‑1 (MQ=255)
                                                                                                                                                                                                                                                         ccttagcgtttCAGCCTGTATTTCGTATTCTCTGGTTTCCCAGACCGCAAAGGGTAGTTCAGGAAAGTATCAATGGCAGTCAGACACCCGGCCTGGGCTAACGCAGACAAAAGATAGTAATGGAGCAAAAGGCTATGTTTTTGAAGGAGAAGCGCTGGATTACCTTATAACAAATGGTGCCAGTGACATCGTTAAGATGCTCAATGCTCCAGCACTTAGCCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAGAGAAATTTACCgg                                     <  2:657580/280‑1 (MQ=255)
                                                                                                                                                                                                                                                           tGATGCTTTCAGGCTGTATCTCGTATGCTGTCGTCTGCCATACAGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCGGATGACGCAAGTTTTGTTTTAAATGCGGGGAAAAAGAAATTTACGGGcac                                  >  1:24878/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                ctatCATGCTGTATCTAGTATGCTGTGTTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCAGCACTGGGCTATAGCAGGCAAAAGGTAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGAAATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGTAAGATTTGTTTTAAATGCGTGGAAAAAGAAATTTACCGGCACAatat                             <  1:978860/286‑1 (MQ=255)
                                                                                                                                                                                                                                                                 tttAAGGCTTTATTTCTTCTGCTATGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACAACATCACTGGGCTATCGCAGGCAAAAGCTAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAatat                             <  2:653851/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                        cTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTAAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCtt                        <  1:541544/286‑1 (MQ=255)
                                                                                                                                                                                                                                                                           tATCTCGTATGCTCTGGTTTCCCATACCGCAAAggg                                                                                                                                                                                                                                                                               <  2:34634/36‑1 (MQ=255)
                                                                                                                                                                                                                                                                           tATCTCGTATGCTCTGGTTTCCCATACCGCAAAggg                                                                                                                                                                                                                                                                               >  1:34634/1‑36 (MQ=255)
                                                                                                                                                                                                                                                                           tATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTtact                    >  2:251928/1‑287 (MQ=255)
                                                                                                                                                                                                                                                                              ctcGTATGCACTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGGCAGACACCATCACTAGGCTATCGCAGGCAAAATATAGTAATGGAACAAAAGGCTATGTTGTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTtactactac               <  2:201210/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                          ggTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCAAAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTTACTACTACTGGAATAACgaa   <  1:646190/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                           gTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTTACTACTACTGGAATAACgaag  >  1:573973/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                
ATGACGTTCGCGTTTGCCGTGCGTGTAATGTAGTACAAACTTATATTGTTGTACTACAATTTAGATCACAAAAAGAACAATGCATAAAAAAATGACATGCGTCGGGCAGAAATCTGAAAAGGGATATCAGGCGCTAAACAGGAGGGAAAGAAGAGTATGCTTTCAACGGCTTAGCTACTCGTTTAAAGGATTAATCATGAAGTTGAATTTTAAGGGATTTTTTAAGGCTGCCGGTTTATTCCCACTGGCGCTGATGCTTTCAGGCTGTATCTCGTATGCTCTGGTTTCCCATACCGCAAAGGGTAGTTCAGGAAAGTATCAATCGCAGTCAGACACCATCACTGGGCTATCGCAGGCAAAAGATAGTAATGGAACAAAAGGCTATGTTTTTGTAGGGGAATCGCTGGATTACCTTATCACTGATGGTGCCGATGACATCGTTAAGATGCTCAATGATCCAGCACTTAACCGGCACAATATTCAGGTTGCCGATGACGCAAGATTTGTTTTAAATGCGGGGAAAAAGAAATTTACCGGCACAATATCGCTTTACTACTACTGGAATAACGAAG  >  NC_000913/3875242‑3875813

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: