Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,123,930 A→T V176E (GTA→GAA)  cytR ← DNA‑binding transcriptional repressor CytR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,123,9300AT98.9% 334.1 / ‑4.5 91V176E (GTA→GAA) cytRDNA‑binding transcriptional repressor CytR
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (32/58);  total (32/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

AGCGCCATCACATCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTACTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCC‑TTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCAGATATCCGTGATTTG  >  NC_000913/4123684‑4124176
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                          
aGCGCCATCACATCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGc                                                                                                                                                                                                                                                        <  1:664860/251‑1 (MQ=255)
   gcCATCACATCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATc                                                                                                                                                                                                                                                     >  1:257383/1‑251 (MQ=255)
         acaTCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATgc                                                                                                                                                                                                                                               <  1:406977/251‑1 (MQ=255)
             cGCTATGGCAGAAGACAGCAGTAGGAGGCTGTGGAAGATCAAGCCGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGGCACAGCGGCCTCCCTTC‑GGGACCGGCTATACAGCCAATCCCTTTTTGCCCTTTCTCCTTTAAATAAATTTCTGCATCAAAAgcggcg                                                                                                                                                                                                                                           >  2:726559/1‑251 (MQ=255)
              gCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACCGCGGCATCTCTTC‑GGGACCGGCTATCCAGCCAAACCGTTTATGCCCTTTCTCATATAAATAATTTTTTTCATCACATGcggcgg                                                                                                                                                                                                                                          >  2:184338/1‑251 (MQ=255)
               ggatggcgaaagtcggcagTAGGCGGATGTGGAACATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGGAGGTGAAGTCGCCACGGGCGATGTAGTGCTGATCATCCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCGTGCCGGCGGTAGTGACACAGCGGCACCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATTCGCTTGCTCATATAAATTATTTTCTGCATCAAATGCGGCGG‑t                                                                                                                                                                                                                                        <  2:67508/236‑1 (MQ=255)
                 aTGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCCATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑Tcc                                                                                                                                                                                                                                      >  1:126626/1‑250 (MQ=255)
                 aTGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCa                                                                                                                                                                                                                                      <  1:105734/251‑1 (MQ=255)
                       gaagaCAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGACGCCACGGGCGCTGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACCCAGCGGCATTTCTTCCGGGACCG‑CTATACAGCCAATACGTTTTTGCCCCTTCTCATATAAACAATTTTCTGCCTCCAAAGCGGCGG‑TCAgctctg                                                                                                                                                                                                                                >  2:385200/1‑246 (MQ=255)
                        aagaCAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCAGCAAATGCTGCGG‑TCAGATTGTc                                                                                                                                                                                                                               <  2:111195/251‑1 (MQ=255)
                         tgcccttcttgcatgttagtgGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGGGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTc                                                                                                                                                                                                                               <  1:454910/232‑1 (MQ=255)
                            cagcagTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGGGGCACAGCGGCATCTCTTC‑GGGACCGGCGATACAGCCCATCCGTTTATGCCCCTGCTCATATAAATAATTTTTTTCATCAAATGcggcgggctctatgttgcg                                                                                                                                                                                                                             >  2:430904/1‑237 (MQ=255)
                                      cGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAg                                                                                                                                                                                                                <  1:118676/251‑1 (MQ=255)
                                       ggTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAAGGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAgg                                                                                                                                                                                                               >  1:15508/1‑251 (MQ=255)
                                           gtgGAAGAGCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGcagc                                                                                                                                                                                                           <  2:483375/251‑1 (MQ=255)
                                                aGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACAAAGCGGCATCTCTTC‑GGGGCCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCCTCAAATGCGGCGG‑TCCGGTTGTCGata‑ttgaatgttagccaccccac                                                                                                                                                                                                      >  1:289952/1‑232 (MQ=255)
                                                       gcagcgtctgcATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGcagctcc                                                                                                                                                                                                        <  2:669707/235‑1 (MQ=255)
                                                       gcagcTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑GCAGATTGTCGCTA‑TGGAATGTAGGgagctcccgatcggg                                                                                                                                                                                                >  1:669707/1‑242 (MQ=255)
                                                        cagcTGCTGCATCGCTTTGCTTCCGGCTTCGCAGGTGAAGTCGCCACGGGCGATGTATTGCGGCTCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGGGACACAGCGGCCTCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCCTGCTCATATAAATAATTTTCTGCATCCAATGCGGCGGGCCACAT‑GTCGATA‑TGAACTGTTGGcaagccccgcccccggg                                                                                                                                                                                              >  2:670609/1‑236 (MQ=255)
                                                            ggctgcATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATGTTCGGCATCAAAGGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCaaa                                                                                                                                                                                         <  2:698656/250‑1 (MQ=255)
                                                             gctgcATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAAt                                                                                                                                                                                        <  2:55154/251‑1 (MQ=255)
                                                                gcATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGAGTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAAttcg                                                                                                                                                                                     <  2:325297/251‑1 (MQ=255)
                                                                gcATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGAGACAAAGCGGCCTCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGGCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCCGATTGTTCCTA‑TCAAATTGAGGCAGATCCAGCCCC‑GGGGCAAAttcc                                                                                                                                                                                     >  1:648078/1‑250 (MQ=255)
                                                                gcATCGCTTTGCTTCCAGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAAtt                                                                                                                                                                                       <  1:380074/249‑1 (MQ=255)
                                                                       tttGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCCATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑GCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCCGCTCC‑GGGGCAAATTCGTTTGccc                                                                                                                                                                              >  1:179209/1‑250 (MQ=255)
                                                                       tttGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGcc                                                                                                                                                                               <  2:126626/250‑1 (MQ=255)
                                                                        ttGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCCGCTCCAGCTCC‑GGTGCAAATTCGTTCGccct                                                                                                                                                                             >  1:647276/1‑249 (MQ=255)
                                                                        ttGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGcc                                                                                                                                                                               <  1:420808/249‑1 (MQ=255)
                                                                         tGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAAGGCCGCAGCGACGCAGCGCCTGACCATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCGTC‑GGGAGCGGCTATGCAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCAt                                                                                                                                                                             <  1:360687/250‑1 (MQ=255)
                                                                          gCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATc                                                                                                                                                                            <  2:219628/250‑1 (MQ=255)
                                                                            ttCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGAAACAGCGGCATCTCTTC‑GGGGCCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TTGACAGTAGGCAGCTCCAACACC‑GGGGCAAATTCGTTTGCCCTCAcc                                                                                                                                                                         >  1:349699/1‑251 (MQ=255)
                                                                             tCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCAc                                                                                                                                                                          <  2:123110/249‑1 (MQ=255)
                                                                               cGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTACTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATc                                                                                                                                                                      <  2:657877/251‑1 (MQ=255)
                                                                                 gCTTCGAAGGGGAAGTCGCCACTGTCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGAATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCTGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAAGTCGTTCGCCATCGCCATcgg                                                                                                                                                                    <  2:285156/251‑1 (MQ=255)
                                                                                   ttCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTTTCGATA‑TGAACTGTAGGCAGCTCCCGCTCC‑CGTGCAAATTCGTTCGCCATCAcccacggcg                                                                                                                                                                  >  1:168972/1‑244 (MQ=255)
                                                                                    tCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCGCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATcggcg                                                                                                                                                                  <  1:725663/250‑1 (MQ=255)
                                                                                        aGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGAACCACCGGCATCTCTTC‑GGGGCCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCCAATGCGGCGG‑CCACATTTTCGATA‑TGAACTGTTGGCAGATCCCGCCCCCGGGCCAAATT‑GTTCGCCCTCACCCTCggcggc                                                                                                                                                                >  1:688257/1‑248 (MQ=255)
                                                                                          gTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATcggc                                                                                                                                                                   <  1:438077/243‑1 (MQ=255)
                                                                                              aGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATCGCCTTGCAGGCGGTACTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGGTA‑TGAACTGTAGGCAGCTCCAGCCCCCGGGGCAA‑TTCGTTCGCCATCACCATCCGCGGCAGATTAcac                                                                                                                                                       >  1:648809/1‑249 (MQ=255)
                                                                                               gTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACt                                                                                                                                                                                                                    <  2:554642/190‑1 (MQ=255)
                                                                                                   ccACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTc                                                                                                                                                  <  1:488061/251‑1 (MQ=255)
                                                                                                      cggtcGATGTATTGCGGATCAACCATAAGGCCGCCGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTTACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCTAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAAGTCGTTCGGCATCACCATCGGCGGCAGATTACGCTGTTCCTc                                                                                                                                               <  2:289952/247‑1 (MQ=255)
                                                                                                            aTGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCGTCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCt                                                                                                                                         <  1:232023/251‑1 (MQ=255)
                                                                                                                aTTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAAGTCGTTCGCTATCGCCATCGGCGGCAGAGTACGCTGTTCCTCAAGGCTGGCa                                                                                                                                     <  2:314563/251‑1 (MQ=255)
                                                                                                                  tGCGGAACAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTGC‑GGGACCGGCGATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAAGGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGGCATCTCCAGCGGCGGCAGATTACGCTGTTCCTCAAGGCTGGCAt                                                                                                                                    <  2:81062/250‑1 (MQ=255)
                                                                                                                   gCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCCGCTCC‑GGGGCAAATTCGTTTGCCCTCACCCTCGGCGGCAGATTAAGCTGTTCCTCAATGCTGGCATc                                                                                                                                   >  1:592510/1‑250 (MQ=255)
                                                                                                                      gATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAAc                                                                                                                               <  1:297885/251‑1 (MQ=255)
                                                                                                                      gATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAAc                                                                                                                               <  1:688972/251‑1 (MQ=255)
                                                                                                                      gATCAACCAGAATGCCGCAGCGACGCAGCGCCTGAACAGAGCCTTGCAGGCGGTAGTGACACAGCGGCAGCTCTTC‑GGGACCGGCTATACAGCCAAGCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑CGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGGTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCGGGCATCAAAc                                                                                                                               <  2:322735/251‑1 (MQ=255)
                                                                                                                         cAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAAGGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGc                                                                                                                            <  2:15508/251‑1 (MQ=255)
                                                                                                                            ccATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGc                                                                                                                         <  1:430904/251‑1 (MQ=255)
                                                                                                                                aaTGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGGTTGTCGATA‑TTAACTGTAGGCAGCTCCAGCTCC‑CGGGCCAATTCGTTCGCCCTCACCCTCGGCGGCAGATTACCCTGTTCCTCAATGCTGGCCTCAAACGGCAGCC‑ttt                                                                                                                    >  1:136453/1‑250 (MQ=255)
                                                                                                                                aaTGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCAGCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCC‑tt                                                                                                                     <  1:670609/250‑1 (MQ=255)
                                                                                                                                   gCCGCATCGACGCAGCGCCTGAACAAAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGGTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAAGAATTTTCTGCATCAAATGCGGCGG‑GCAGATTGTCGATA‑TGAACTGTAGACAGCTCCAGCTCC‑GGTGCAAATTCGTGCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCC‑TTGAAc                                                                                                                 <  1:384867/251‑1 (MQ=255)
                                                                                                                                        aGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTg                                                                                                                                                     <  2:651614/211‑1 (MQ=255)
                                                                                                                                        aGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGAAACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCCGAATTTCGATA‑TTAACTGTAGGCAGCTCCCGCCCC‑CGGGCCAATTCGTTTGCCATCCCCATCGGCGGCCGATTTCGCTg                                                                                                                                                     >  1:651614/1‑211 (MQ=255)
                                                                                                                                          cGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCTTTTATGCCCGTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCGGATTGTCGATA‑TTAACTGTAGGCAGCTCCAGCTCC‑GTTGCAAATTCGTGCGCCATCACCATCGGCGGCATATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCC‑TTGAACCCAGcaa                                                                                                          <  1:733222/251‑1 (MQ=255)
                                                                                                                                           gACGCAGCGCCTGAACATAGCCTTGCCGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCC‑GTGAACCCAGcaac                                                                                                         <  2:688257/251‑1 (MQ=255)
                                                                                                                                                   gcCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCC‑TTGAACCCAGCAACAACATGcc                                                                                                 <  1:193603/251‑1 (MQ=255)
                                                                                                                                                     cTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAAATGTAGGCAGCTCCCGCTCC‑CGGGCAAAATTGTTCGCCCTCACCATCGGCGGCAGATTAAGCTGTTCCTCCATGCTGGGCTCCAACGGCAGCC‑TTTAACCCAGCAACAACATGCCAt                                                                                               >  1:253185/1‑251 (MQ=255)
                                                                                                                                                       gAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCGCAATGCTGGCATCAAACGGCAGCC‑TTGAACCCAGCAACAACATGCCATCa                                                                                             <  1:76187/251‑1 (MQ=255)
                                                                                                                                                          cATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCC‑TTGAACCCAGCAACAACATGCCATCAAtt                                                                                          <  1:567441/251‑1 (MQ=255)
                                                                                                                                                           aTAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCC‑TTGAACCCAGCAACAACATGCCATCAAttt                                                                                         <  1:726559/251‑1 (MQ=255)
                                                                                                                                                           aTAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCC‑TTGAACCCAGCAACAACATGCCATCAAttt                                                                                         <  2:349699/251‑1 (MQ=255)
                                                                                                                                                           aTAGCCTTGCAGGCGGCAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCCGATTGTCGATA‑TTAACTGTAGGCAGCTCCAGCTCC‑CGGGCAAATTCGTTTGCCATCACCATCCGCGGCCGATTACGCTGTTTCTCAAAGCTGGCATCAAAACGCAGCC‑CTGAACCCAGCAACAACATGCCATCCAttt                                                                                         >  1:255436/1‑251 (MQ=255)
                                                                                                                                                            tAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATAAAGGCAATCCGGTTATGCCCTTGCTCAAATAAATAATTTTCTGCATCCAAAGCGGAGG‑CCAGATTGTTGATAATGAAC‑GTAGGCAGCTCCCGCACC‑GGGGGAAAATTGTTTGCCCTCAACACCGGCGGGAGATTAAGCTGTTCCTCAAAGCCCGCGTCAAACGGGAGCCATGGAAACCAGCAACAAAAGGCCATAAAttt                                                                                         >  1:674622/1‑251 (MQ=255)
                                                                                                                                                                  ttcAGGCGGTAGTGACACAGCGGCATCTCTTC‑GTGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTTCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCC‑TTGAACCCAGCAACAACATGCCATCAATTTGCTTgg                                                                                   <  2:35978/248‑1 (MQ=255)
                                                                                                                                                                  tGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCC‑TTGAACCCAGCAACAACATGCCATCAATTTGCTTgg                                                                                   <  1:173209/250‑1 (MQ=255)
                                                                                                                                                                    cAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCC‑TTGAACCCAGCAACAACATGCCATCAATTTGCTTGGtga                                                                                <  2:178004/251‑1 (MQ=255)
                                                                                                                                                                      ggCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTCCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGAATTTCGATA‑TGCACTGTAGGCAGCTTCAGCTCC‑GGTGCAAATTTGTTCGCCATCACCATCCGCGGCAGGTTCCGCTGTTTCCCAATGATTGCATCAAAAGGCAGCC‑TTGAACCCCGCAACAACATTCCATCAATTTGCTTGGtgatg                                                                              >  1:459373/1‑251 (MQ=255)
                                                                                                                                                                         ggTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCC‑TTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATc                                                                           <  1:273611/251‑1 (MQ=255)
                                                                                                                                                                            aGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCC‑TTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCaa                                                                         <  2:168972/250‑1 (MQ=255)
                                                                                                                                                                             gTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCC‑TTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAAt                                                                       <  2:136453/251‑1 (MQ=255)
                                                                                                                                                                                      cGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATGCGTTCGCCATCGCCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCC‑TTGAACCCAGCAACAACATGCCATCAATCTGCTTGGGTATGATCAAATCGATAAAgg                                                              <  2:35015/251‑1 (MQ=255)
                                                                                                                                                                                       ggCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTGTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTTTTCCTCAATGCTGGCATCAAACGGCAGCC‑TTTAACCCAGCAACAACATGCCATCAATTTGCTTCGTGATGATCAAATCGATAAAGGt                                                             <  2:459373/251‑1 (MQ=255)
                                                                                                                                                                                        gCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCC‑TTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAAGCGATAAAGGtt                                                            <  2:661313/251‑1 (MQ=255)
                                                                                                                                                                                              ctTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCC‑TTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGttttt                                                         <  2:592510/248‑1 (MQ=255)
                                                                                                                                                                                                    ggACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCC‑TTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTCctgctg                                                 <  2:179209/251‑1 (MQ=255)
                                                                                                                                                                                                             ataCAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCC‑TTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTGCCTGCTGATTTTGAt                                         <  2:647276/250‑1 (MQ=255)
                                                                                                                                                                                                             ataCAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCGCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCAGGCATCAAACGGCAGCC‑TTGAACCCAGCAACAACATGCCATCAATTTGCTGGGTGATGATCAAATCGATAAAGGTTTTTTCCTGCTGATTTTGa                                          <  1:648760/249‑1 (MQ=255)
                                                                                                                                                                                                                 aGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCC‑TTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTCCTGCTGATTTTGATGCGca                                    <  1:385200/251‑1 (MQ=255)
                                                                                                                                                                                                                   ccAATCCGGTTATTCCCGTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCGGTTCCTCAAGGCTGGCATCAAACGGCAGCC‑TTGAACCCAGCAACAACATGCCATCAAGTTGCTTGGGGATGATCAAATCGATAAAGGTGTTTGCCTGCTTATAGGGATGCGcaca                                  <  2:344332/251‑1 (MQ=255)
                                                                                                                                                                                                                            ttATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCC‑TTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAAt                         >  1:706861/1‑251 (MQ=255)
                                                                                                                                                                                                                                gCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑ACAGATTGTCCATA‑TGAACTGTAGGCAGCTCCCGCTCC‑GGCGCAAATTCGTTCGCCATCACCATCGCCGGCAGATTACGCCGGTCCTCAATGCTGGCATCAAACAGCAGAA‑TTGAAACCAGCAACAACATGCCAACAATTTGCTTGGGGATGATCAAATCGATAAAGGTTTTTTCCTGCTGATTTTGATGGGCACAGACGCCAATCAGc                     >  1:456513/1‑251 (MQ=255)
                                                                                                                                                                                                                                 cccTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGTCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCAACATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCAGCAAACGGCAGCC‑TTGAACCCAGCAACAACATGCCATCAAATTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTCCTGCTGATTTTGATTCGCAAAGGCGCCAATCAGc                     >  1:202240/1‑250 (MQ=255)
                                                                                                                                                                                                                                   cTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCC‑TTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTCCTGCTGATTTTTATGCGCACAGTCGCCAATCAGCAcc                  >  1:566580/1‑251 (MQ=255)
                                                                                                                                                                                                                                    ttGCTCATATAAATAATTTTCTGCATCAAATGCGGCGG‑TCCGATTGTCGATA‑TGAACTGTCGGGCGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCCTCGGCGGCAGTTTTCGCTGTTCCTCCATGCTGGGCTCCAACGGCAGCC‑TTGAACCCAGACACCACATGCCATCTATTTGCTTGGGGATGATCAAATCGAAAAAGGTTTTTTCCTGCTGATTTTGATGCGCGCAGTCGCCAATCAGCAcc                  >  2:625236/1‑250 (MQ=255)
                                                                                                                                                                                                                                            ataaataaTTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGGTA‑TGAACTGTAGGCAGCTCCCGCTCC‑GGTGCAAATTCGTTCGCACTCACCATCGGCGGCAGATTACGCTGTTCCTCAAAGCTGGCAACAAACTGCAGCC‑TTGAACCCAGCAACAACATGCCATCAATTTGCTTGGGGATGATCAAATCGATAAAGGGTTTTTACTGCTGATTTTGATGCGCACAGTCGGCCATCAGCAACAGCTATCCg         >  1:314879/1‑251 (MQ=255)
                                                                                                                                                                                                                                              aaataaTTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCC‑TTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCAGATATCCGt        <  2:257383/250‑1 (MQ=255)
                                                                                                                                                                                                                                              aaataaTTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCC‑TTGAACCCAGCAAAAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGGTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCAGATATCCGTg       >  1:668095/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   atTTTCTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCC‑TTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCAGATATCCGTGATTTg  >  1:731974/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                          
AGCGCCATCACATCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTC‑GGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTACTGCATCAAATGCGGCGG‑TCAGATTGTCGATA‑TGAACTGTAGGCAGCTCCAGCTCC‑GGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCC‑TTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCAGATATCCGTGATTTG  >  NC_000913/4123684‑4124176

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