Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,027,435 T→G 3.9% intergenic (+208/+173) fre → / ← fadA NAD(P)H‑flavin reductase/3‑ketoacyl‑CoA thiolase (thiolase I)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,027,4350TG 3.9% 320.7 / 4.5 115intergenic (+208/+173)fre/fadANAD(P)H‑flavin reductase/3‑ketoacyl‑CoA thiolase (thiolase I)
Reads supporting (aligned to +/- strand):  ref base T (61/47);  new base G (2/4);  total (64/51)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.03e-01

Only 100 of 116 total aligned reads displayed.

TGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAA  >  NC_000913/4027209‑4027672
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               
tgggttttgtttttttttgtaaataAAAACACCCCGCCCCGACAGGGGGGAAGAACGGCAACTAAACTGTTATTCAGGGGCATTTAGATCTATGACGTAACTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATc                                                                                                                                                                                                                                      <  1:701341/216‑1 (MQ=18)
 ggttttgaatttttcgagaataaAAAAAACCCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCg                                                                                                                                                                                                                                     <  2:872966/219‑1 (MQ=21)
 ggcgttgactttgcctgaaaaaaaaaaacaCCCGCCCCGGGCAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATTTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGGCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCg                                                                                                                                                                                                                                     <  2:672019/215‑1 (MQ=21)
          ttttccggaaataaAAAAAACCCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACaaa                                                                                                                                                                                                                           <  2:133659/229‑1 (MQ=255)
          tttATCATAGATAAAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACaa                                                                                                                                                                                                                            <  2:309252/238‑1 (MQ=255)
            tATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCCTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCCCATCCGGCAATTGGGGCATGATGCCTGATGCGACGCTGTGCCGTCTTATCGTGCCTACACATa                                                                                                                                                                                                                         >  2:78128/1‑239 (MQ=255)
                  agaTATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGa                                                                                                                                                                                                                 >  2:148825/1‑239 (MQ=255)
                      aaaaaaaacCCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCg                                                                                                                                                                                                             <  2:406762/237‑1 (MQ=255)
                          aaaacCCCGCCCCTGACAGGCGGGAAGAACGGCAACGAAACTTTTAGTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAgg                                                                                                                                                                                                         <  2:664475/239‑1 (MQ=255)
                              aCCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCCATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCgg                                                                                                                                                                                                     >  2:365503/1‑239 (MQ=255)
                              aCCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTGGCGGCTTATCGTGCCTAACAAAAGTCCGAAACGTAGGCCgg                                                                                                                                                                                                     >  2:303764/1‑239 (MQ=255)
                                   ccccTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCGAACCATGCCCGATGCGCCGCTGTCGCGTCTTATCGTGCCTACAAATAATCCGAACCGGAGGACGGATAAGGCGTTTTCGCCGCATCCGCCAATTGGTGCATGAAGCCTGATGCGAACCCTGTCCGCTCTAACCGGCCCACAAAAAGTCCGAAACCGAGGCCgggcaac                                                                                                                                                                                                >  1:388200/1‑234 (MQ=21)
                                   ccccTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAAGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTATGCCGGATaaa                                                                                                                                                                                                >  2:357251/1‑238 (MQ=255)
                                       tGACAGGCGGGAAGAACGGGAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGt                                                                                                                                                                                            <  2:445649/239‑1 (MQ=255)
                                              cGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCTTAGGCCGGGTAAGCGGTTTAccc                                                                                                                                                                                      >  2:962334/1‑236 (MQ=255)
                                               gggAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCAAGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTAcgcc                                                                                                                                                                                     <  2:846402/238‑1 (MQ=255)
                                                                     ttATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAACAAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCCCATCCGGCCATTGGTGCatga                                                                                                                                                               >  2:494386/1‑238 (MQ=255)
                                                                       aTTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCCTGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGCAGGCCGGAAAAGGGGTTTACCCCCCACCCGGACATTtgtgaaaagac                                                                                                                                                             >  1:387616/1‑231 (MQ=21)
                                                                       aTTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGc                                                                                                                                                                                                        <  1:435226/195‑1 (MQ=255)
                                                                       aTTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTCGGCCGCATCCGCCCATTGGTGCataatg                                                                                                                                                             >  2:374789/1‑238 (MQ=255)
                                                                             tGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGTTGCCTGATGCGACGCTGTCGCGTCTTATCGTGcc                                                                                                                                                                                                                                 >  1:253226/1‑166 (MQ=255)
                                                                             tGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGGGCTTACAAAAATCCCGAACCGTAGGCCGGGTAAGGCCTTTACCCCCCCACCGCCAAttcgggctaaacgccgcttc                                                                                                                                                      >  1:306494/1‑221 (MQ=18)
                                                                             tGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCCCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGCCTTATCGTCCCTCCAAATAGTCCGAACCTTAGGCCGGGAAAGGCTTTTACCCCCACCCCGGCAATTgggccataatcccgaat                                                                                                                                                       >  1:462506/1‑223 (MQ=18)
                                                                               gCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGc                                                                                                                                                     >  2:47552/1‑238 (MQ=255)
                                                                                 aTTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGGGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGttaaacccgctcaaacac                                                                                                                                                    <  2:911202/238‑18 (MQ=14)
                                                                                  tttAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATACGGCGTTCACGCCGCATCCGGCAAGTGGttaaacccgctcaaacacc                                                                                                                                                   >  2:490141/1‑220 (MQ=14)
                                                                                      attttttggcgtttttGGGAACCCATGCCTGATGCGGCGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTgg                                                                                                                                             <  1:47552/224‑1 (MQ=18)
                                                                                       aTCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGTGTGCAGGAGTCCTGGTGCCACCCTGgc                                                                                                                                            >  1:355728/1‑239 (MQ=21)
                                                                                         cTATGCCGTATCGGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGgcgc                                                                                                                                          <  2:410659/239‑1 (MQ=255)
                                                                                           aTGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCCGATGCGACGCTGgcgc                                                                                                                                          >  2:264163/1‑237 (MQ=255)
                                                                                            tGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGCGGTTTACGCCGCATCCGGCAATTGGGCCATGATGCCTGATGCGACGCTgggcggtg                                                                                                                                       >  2:581808/1‑233 (MQ=255)
                                                                                                gtttCTGGCAAACCAGGCCCGATGCTACGTTGTCCCGTCTTATCGATCCTACAAATAGGCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                   <  2:599431/236‑1 (MQ=25)
                                                                                                gttgtttggaaaaCATGCCCGTTGGCACGCTGTCGCGTCTTATCCGGCCTAAAAATAGCCCGACCCTTAGGCCGGATAAGGGGTTTCCGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                   <  1:746817/230‑1 (MQ=18)
                                                                                                gTATCTGGCAAACCATGCCCGATGCGACGCTTGCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                   <  2:316581/239‑1 (MQ=255)
                                                                                                 tATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGTCCGGGTCAGGCGTTTACGCCGCATCCGGCAATTGTTGCATGATGCCTGATGGGAGGCTGGCGCGTCTTTTc                                                                                                                                  >  2:1001391/1‑239 (MQ=255)
                                                                                                 tATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATc                                                                                                                                  >  2:287237/1‑239 (MQ=255)
                                                                                                 tATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGCCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATAGTGCCTGAGGCGACGCGGgcgggcgttata                                                                                                                                  >  1:953818/1‑230 (MQ=18)
                                                                                                   tCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGACAAGGCGTTTACGCCCCATCCGGCCATTGGGGCATGTTGCCTAAACCGACGCTGCCCCGCCTTacctg                                                                                                                                >  2:530904/1‑235 (MQ=18)
                                                                                                           tcctttcccGGTGCGCGGCTGTCGGGTTTTTTCGTGGCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGTGGTGCATGATGCCTGATGCGGCGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAt                                                                                                                        <  1:287237/233‑1 (MQ=21)
                                                                                                                 cccGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGGGCATGAGGCCTGATTCGCCGCTGGCCCGCCTTATAAGGCCTACATTgggg                                                                                                                   >  2:510791/1‑236 (MQ=255)
                                                                                                                   cGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCg                                                                                                                <  2:796149/239‑1 (MQ=255)
                                                                                                                      tGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATATTCCGGACCGTACGCCGGCTAAGGGGTTTACCCCGCATCCGGACATTGGTGCATGATGCCTGGTGCGACGCCTGCGCGTCTTTCCAAGCCTACAtagtggcgggac                                                                                                             >  1:231628/1‑228 (MQ=18)
                                                                                                                      tGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACCCTGGCCCGGCTTTTCAGGCCTCCATTGGGGCCggtt                                                                                                             >  2:958521/1‑237 (MQ=255)
                                                                                                                         gaaccTGTCGCGTCTTATCGTGCCTACAAATTGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCg                                                                                                           <  2:1018670/234‑1 (MQ=255)
                                                                                                                            gttgttgctttgtagtgTGCCTAAAAATAGGTCAACACGGAGGCCGGAATAGGGGTTTACGGCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAg                                                                                                       <  1:724391/223‑1 (MQ=18)
                                                                                                                            gctgtcggttttttgcgggCCTAAAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAg                                                                                                       <  1:365503/221‑1 (MQ=255)
                                                                                                                            cgttgcgggttttttggTGCGTAAAAAGAGGCCGGACCGTAGGGCGGATAAGGCGTTTGCGGCCGCTCCGGAAATTGGTGCAGGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAg                                                                                                       <  1:78128/223‑1 (MQ=18)
                                                                                                                              tGTCGCGTCTTATCGTGCCTTCAAATAGACCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGa                                                                                                      <  2:842131/238‑1 (MQ=255)
                                                                                                                                tCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAGTAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGGCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCg                                                                                                   >  2:141174/1‑239 (MQ=255)
                                                                                                                                tCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGGGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCCGGGCTACATGGGTGCGAGACGGGTAGAACg                                                                                                   >  1:658829/1‑239 (MQ=21)
                                                                                                                                               ccTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGGAGGCCGGAAAAGGCGTTTACGCCGCATCCGGCAATTGGGGCATAATGCCCGATGCGAGGCTGGCGACGCTTATCCGGCCAAAACTGTGGCCGGCTCGGGAACCCGGAAAGGGGGTTCagg                                                                                    >  1:849085/1‑237 (MQ=18)
                                                                                                                                               ccTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGCCCGAACCGTAGGCCGGCAAAGGCGGTTACGCCGCATCCGGCCATTGGGGCATGATGCCCGATGCAACGCTGTCCGGTCTTATCAGCCCTCAATTGGGGACGGATGGGGCGACCGGGCAAgggggttccc                                                                                    >  1:904489/1‑231 (MQ=18)
                                                                                                                                                    aaaTAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCAt                                                                                                                                                                                 <  2:426087/141‑1 (MQ=255)
                                                                                                                                                          tCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGAc                                                                                                     <  1:510374/211‑1 (MQ=255)
                                                                                                                                                          tCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGc                                                                         >  2:823863/1‑239 (MQ=255)
                                                                                                                                                          tCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATAGGTAGACCGGATAAGGCGTTCACGCCCCATCCGGc                                                                         >  2:359868/1‑239 (MQ=255)
                                                                                                                                                          tCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATCATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGCGCCGGACCGGTAGACCGGATAAGGCGTTACCTCCGCATCCGGc                                                                         >  2:466081/1‑239 (MQ=255)
                                                                                                                                                          tCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGGGCCTACAATTAGTCCGAACCGTAGGCCGGAGAAGGCGTTTACGCCGCATCCGGCCAATGGGGCATGATGCCCGATGCGACGCTGGCGCAGCCTATCACGCCTCCATTGTGGCCGGTCCGGTGGACCGGGAAAgggtttacggccccacccgcc                                                                         >  1:361880/1‑220 (MQ=18)
                                                                                                                                                            cGACCCGTTGGGCGGTTAGGGTGTTTCCCCGCACTCGGTCAATTGGGGCAGTAGGCCTGAGGGGACGGTGGCGCTGCTTATCGGGCCTAAAAATAGTCCGACCCGTAGGCCGGATAAGGCGTTTACGGCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCaa                                                                       <  1:141174/239‑1 (MQ=17)
                                                                                                                                                            cGAACCGTAGGCCGGATAAGGCGTTTCCGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCaa                                                                       <  2:788740/239‑1 (MQ=255)
                                                                                                                                                            cGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAg                                                                                                       >  2:408999/1‑207 (MQ=255)
                                                                                                                                                            cGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCaa                                                                       <  2:904489/239‑1 (MQ=255)
                                                                                                                                                            cGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCaa                                                                       <  2:231628/239‑1 (MQ=255)
                                                                                                                                                            cGAACCGTAGGCCGGATAAGGCGTTTACCCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCaa                                                                       <  2:50526/239‑1 (MQ=255)
                                                                                                                                                              aaCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGGGGTTACGCCCCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGCCTTTCCAGGCCCACATTGGTGCCGGAACGGTAGACCGGGAAAGGGGTTCACCCCACATCCGGCAAgg                                                                     >  1:788740/1‑238 (MQ=21)
                                                                                                                                                              aaCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTg                                                                                                                                                                       >  1:426087/1‑141 (MQ=255)
                                                                                                                                                              aaCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGt                                                                     >  2:361411/1‑239 (MQ=255)
                                                                                                                                                                   tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGCCGCGCCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGCTAAGGCGTTAACTCCCCATCCGGCAAGTGGTTa                                                                >  2:240922/1‑239 (MQ=255)
                                                                                                                                                                    aGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATaat                                                                                                                                                                                                                       >  1:660095/1‑87 (MQ=255)
                                                                                                                                                                    aGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACa                                                                                                                         >  2:601542/1‑181 (MQ=255)
                                                                                                                                                                           aTAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                   >  1:1005297/1‑164 (MQ=255)
                                                                                                                                                                           aTAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                   <  2:1005297/164‑1 (MQ=255)
                                                                                                                                                                           aTAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAg                                                                                                       >  2:713160/1‑192 (MQ=255)
                                                                                                                                                                           aTAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAg                                                                                                       <  1:713160/192‑1 (MQ=255)
                                                                                                                                                                           aTAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGACCGGTAGACCGGATAAGGCGTTCACGCCGCTTCCGGCAAGTGGTTAAACCCGcc                                                        >  2:416168/1‑238 (MQ=255)
                                                                                                                                                                             aaGGCGTTTACGCCGCATCCGGCAATGGGTGCATGATGCCTGATGCGACGCTGTCGCGTTTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCa                                                      <  2:658829/239‑1 (MQ=255)
                                                                                                                                                                                   tttCCGCCGCATCCGGCAATTGGTGCATGATGCCTGATGGGACGTTGTCGCGTCTTAGCGGGCCTACAAATAGTCCGAACCGTGGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAacac                                                 <  1:510791/238‑1 (MQ=255)
                                                                                                                                                                                   tttACGCCGCATCCGGCAATTGGTGCATGATGCTTGATGGGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAacac                                                 <  2:202559/238‑1 (MQ=255)
                                                                                                                                                                                    ttACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGTGTCCGGATCGGTCGACCGGAAAAGGCGTTCACCCCCCCTCCGGCCAATGGGTAAACCCCCCCAAacaac                                                >  1:138842/1‑236 (MQ=18)
                                                                                                                                                                                    ttACGCCGCATCCGGCAATGGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCg                                               <  2:323282/239‑1 (MQ=255)
                                                                                                                                                                                     tACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGAATGATGCCTGATGCGACGCTGGCGCGCCTTCTCCGGCCTACATTCGTGCCGGACCGGTCGACCGGACAAGGGGTTCAACCCCCAGCCGGCCACTGGTTCAAACCGCCCAAACAccc                                               >  1:219860/1‑237 (MQ=18)
                                                                                                                                                                                      atacctaacctCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCCACTCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGACTTATCAGGCCGACATTGTGGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCACCCCAAacaaca                                               >  1:202559/11‑236 (MQ=18)
                                                                                                                                                                                           cgcATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGCCCTACATTGGTGCCGGATCGGCAGACCGGATAAGGGGTTCCCGCCGCATCCGGCAAGGGGTTAAACCCGCCCCAAAACCATCGCaa                                         >  1:917612/1‑239 (MQ=21)
                                                                                                                                                                                           cccattCCGGCTATTGGGTATGTATGCCTGATGCGACGTTGGCCGGTCTTATCGGGCCTACAAATGGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGCCCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCa                                          <  1:303764/234‑1 (MQ=18)
                                                                                                                                                                                                tCCGGCAATTGGTGAATGTGGCCTGATGCGACCGTGTCGCCTCTTTGCGTGCCTAAAAAAAGTCCGAACCGTAGGCCGGATAAGGCGTTTCCGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATAcc                                     <  1:264163/239‑1 (MQ=21)
                                                                                                                                                                                                tCCGGAATTTGGTGCATGGTGCCGTATGGGGCTCTGGTCGGTCTTATCGGGCCTACAAATAGTCCGCACCGGAGGCCGGATAAGGCGTTTCCGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATAcc                                     <  1:840206/239‑1 (MQ=21)
                                                                                                                                                                                                        ttGGTGCATGATGCCTGATGCGACGCGGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCa                             <  2:953818/239‑1 (MQ=255)
                                                                                                                                                                                                          ggTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGa                           <  2:430810/239‑1 (MQ=255)
                                                                                                                                                                                                            tGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGAc                          <  2:917612/238‑1 (MQ=255)
                                                                                                                                                                                                                   atgCCTGATGCGACGCTGTCGCGTCTTATCGGGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                                                                                                                                                             <  1:649887/85‑1 (MQ=18)
                                                                                                                                                                                                                     gCCTGATGCGACGTCGTCGCGTCTTATCTTGCCTAAAAATAGTCCGAACCGTAGGCCGGAAAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATAcc                                     <  2:34203/219‑1 (MQ=255)
                                                                                                                                                                                                                       cTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATTCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGGGTTCACCCCGCATCCGGCAAGGGGTTAAACCCGCCCAAACACCCTCGCAAAACCCTGACCCAAACCGAAacaaaca               >  1:967742/1‑235 (MQ=21)
                                                                                                                                                                                                                           tGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGACCGGTAGACCGGCTAAGGCGTTCACGCCGCATCCGGCAAGGGGTTAAAACCGCTCAAACACCGTCGCAATACCCTGCACCAACCCGATACACAtcgtc            >  2:36988/1‑238 (MQ=21)
                                                                                                                                                                                                                                aCGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCCAACCCCTAACCCAGACCGATACACAtcgcgcccg        >  2:86377/1‑232 (MQ=255)
                                                                                                                                                                                                                                aCGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGGGGTTAAACCCGCTCAAACACCCCCCCAATACCCTGACCCAGACCGAAACACAccgcccccca       >  1:128313/1‑232 (MQ=21)
                                                                                                                                                                                                                                aCGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGAAAACCGGATAAGGCGGTCACGCCGCATCCCGCCAGTGGTTAAACCCGCCCACAAACCTCCCCCATACCCTGACCCAGTcccgtcccttaacaaaacat       >  1:34203/1‑221 (MQ=18)
                                                                                                                                                                                                                                 cGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGCTAAGGCGTTCACGCCGCATCCGGCAAGGGGTTAAACCCGCCCAAACACCGTCGCCATACCCTGACCCAGACCGATACCAATCGCCGCCAg       >  1:188131/1‑238 (MQ=255)
                                                                                                                                                                                                                                 cGCCGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGAc                                                                                                     >  2:686128/1‑144 (MQ=255)
                                                                                                                                                                                                                                 cGCCGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGAc                                                                                                     <  1:686128/144‑1 (MQ=255)
                                                                                                                                                                                                                                       tCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACCTCGTCGCCAGACCaa  >  2:537897/1‑238 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               
TGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAA  >  NC_000913/4027209‑4027672

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