Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,007,486 A→T 7.8% K541N (AAA→AAT) ‡ recQ → ATP‑dependent DNA helicase RecQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,007,4860AT 7.8% 306.0 / 13.4 102K541N (AAA→AAT) ‡recQATP‑dependent DNA helicase RecQ
Reads supporting (aligned to +/- strand):  ref base A (47/47);  new base T (4/4);  total (51/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-01

Only 100 of 103 total aligned reads displayed.

TACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTT  >  NC_000913/4007338‑4007635
                                                                                                                                                    |                                                                                                                                                     
tACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAAc                                                                                                                                                      <  1:928719/150‑1 (MQ=255)
tACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAAc                                                                                                                                                      <  1:934535/150‑1 (MQ=255)
tACAACTGACAGAGGCCGCGCGCCCGGTGCTCAGCGGCGAATCCTCATTGCAACTTGCCGTGCCGAGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATAGCAAACTGTTCGCCAAATTACGCATAc                                                                                                                                                      >  1:955546/1‑150 (MQ=255)
  cAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTg                                                                                                                                                    >  2:284766/1‑150 (MQ=255)
       gACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTaa                                                                                                                                               >  1:1494651/1‑150 (MQ=255)
       gACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTaa                                                                                                                                               >  2:1171190/1‑150 (MQ=255)
          agagGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATc                                                                                                                                            <  2:275713/150‑1 (MQ=255)
          agagGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATc                                                                                                                                            >  2:883230/1‑150 (MQ=255)
          agagGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATc                                                                                                                                            <  1:355769/150‑1 (MQ=255)
          agagGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATc                                                                                                                                            <  2:292591/150‑1 (MQ=255)
          agagGCCGCGCGCCAGGTGCTGCGCGGCTAATCCTCTTTGCAACATGCCGTGCCGCGTATCTTGTCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATc                                                                                                                                            <  1:482670/150‑1 (MQ=255)
             ggCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAGTTTGCGTaatttggc                                                                                                                                         <  2:216795/150‑7 (MQ=255)
              gCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATa                                                                                                                                        <  1:284766/150‑1 (MQ=255)
                cgcgcgCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGc                                                                                                                                      >  2:959561/1‑150 (MQ=255)
                 gcgcgcCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGcc                                                                                                                                     <  1:978608/150‑1 (MQ=255)
                       ccGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGaa                                                                                                                               <  1:220219/150‑1 (MQ=255)
                         ggTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAg                                                                                                                             <  1:193489/150‑1 (MQ=255)
                         ggTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAg                                                                                                                             <  2:741456/150‑1 (MQ=255)
                            gctgcGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAGTTTGCGTaatttgg                                                                                                                                          <  2:761802/134‑6 (MQ=255)
                            gctgcGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAGTTTGCGTaatttgg                                                                                                                                          >  1:761802/1‑129 (MQ=255)
                            gctgcGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTaa                                                                                                                          <  1:122193/150‑1 (MQ=255)
                            gctgcGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTaa                                                                                                                          <  1:1636150/150‑1 (MQ=255)
                                 gcggcgAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTcc                                                                                                                     >  1:1281688/1‑150 (MQ=255)
                                   ggcgAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCccg                                                                                                                   >  2:27864/1‑150 (MQ=255)
                                    gcgAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCccgc                                                                                                                  >  1:377657/1‑150 (MQ=255)
                                      gAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAGTTTGCGTaatttgg                                                                                                                                          <  1:1023598/124‑6 (MQ=255)
                                      gAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAGTTTGCGTaatttgg                                                                                                                                          >  2:1023598/1‑119 (MQ=255)
                                        aTCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTa                                                                                                              >  2:1343049/1‑150 (MQ=255)
                                          ccTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACg                                                                                                            <  2:1270360/150‑1 (MQ=255)
                                          ccTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACg                                                                                                            >  2:517381/1‑150 (MQ=255)
                                           ctctTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACgt                                                                                                           <  2:1281982/150‑1 (MQ=255)
                                              tttGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATTCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGGAAACTGCGTAAATCGATAGCCGATGAAAGTATTGTGCCGCCGTACgtggt                                                                                                        <  2:1655172/150‑1 (MQ=255)
                                               ttGCAACTTGCCGTGCCGCATATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACgtggtg                                                                                                       >  1:675854/1‑150 (MQ=255)
                                                tGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGgtgt                                                                                                      >  1:1516022/1‑150 (MQ=255)
                                                 gCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGtt                                                                                                     <  1:846494/150‑1 (MQ=255)
                                                  cAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGttt                                                                                                    <  2:257193/150‑1 (MQ=255)
                                                  cAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGttt                                                                                                    <  2:1032665/150‑1 (MQ=255)
                                                   aaCTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTa                                                                                                   <  2:185012/150‑1 (MQ=255)
                                                    aCTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACgtgg                                                                                                         >  2:29024/1‑143 (MQ=255)
                                                    aCTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACgtgg                                                                                                         <  1:29024/143‑1 (MQ=255)
                                                      ttGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAacg                                                                                                <  2:286861/150‑1 (MQ=255)
                                                       tgccgtgccgCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCTGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAacga                                                                                               >  2:784666/1‑150 (MQ=255)
                                                       tgccgtgccgCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAacgc                                                                                               >  2:1486046/1‑149 (MQ=255)
                                                        gccgtgccgCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACgtgg                                                                                                         >  2:952103/1‑139 (MQ=255)
                                                        gccgtgccgCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACgtgg                                                                                                         <  1:952103/139‑1 (MQ=255)
                                                           gtgccgCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCa                                                                                           <  1:520644/150‑1 (MQ=255)
                                                               cgcgTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCt                                                                                       <  1:1580268/150‑1 (MQ=255)
                                                                gcgTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCtt                                                                                      >  1:96987/1‑150 (MQ=255)
                                                                   tATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCttgat                                                                                   <  2:807988/150‑1 (MQ=255)
                                                                    aTCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCttgatt                                                                                  >  2:1247773/1‑150 (MQ=255)
                                                                     tCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCttgattg                                                                                 <  2:1541985/150‑1 (MQ=255)
                                                                     tCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCttgattg                                                                                 <  2:1281688/150‑1 (MQ=255)
                                                                      cGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCttgattga                                                                                >  2:234340/1‑150 (MQ=255)
                                                                       gTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTgag                                                                               >  2:822305/1‑150 (MQ=255)
                                                                          gcgcTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCATACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATg                                                                            <  2:132809/150‑1 (MQ=255)
                                                                           cgcTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGTCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATgg                                                                           >  1:1589376/1‑150 (MQ=255)
                                                                            gcTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGc                                                                          >  2:84562/1‑150 (MQ=255)
                                                                            gcTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGGGTTTAACGACGCAACCTTGATTGAGATGGc                                                                          <  2:1093583/150‑1 (MQ=255)
                                                                               cAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGa                                                                       <  1:1170235/150‑1 (MQ=255)
                                                                                 aaCCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAAc                                                                     >  2:732579/1‑150 (MQ=255)
                                                                                 aaCCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAAc                                                                     <  1:959561/150‑1 (MQ=255)
                                                                                   ccGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAg                                                                   >  1:112491/1‑150 (MQ=255)
                                                                                     gAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGAt                                                                 >  2:879218/1‑150 (MQ=255)
                                                                                         gCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCg                                                             >  1:1272954/1‑150 (MQ=255)
                                                                                          cGATGCAGAAATCGTTCGGAGGCAACTATGATAGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGa                                                            >  1:583686/1‑150 (MQ=255)
                                                                                           gATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGAt                                                           <  1:517381/150‑1 (MQ=255)
                                                                                           gATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGAt                                                           >  1:961481/1‑150 (MQ=255)
                                                                                            aTGCAGAAATCGTTCGGCGTCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATc                                                          <  2:1045197/150‑1 (MQ=255)
                                                                                            aTGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATc                                                          <  2:1142827/150‑1 (MQ=255)
                                                                                              gCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCAc                                                        >  2:1226129/1‑150 (MQ=255)
                                                                                               cAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCAcc                                                       >  2:1221863/1‑150 (MQ=255)
                                                                                                 gAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGc                                                     >  2:1201305/1‑150 (MQ=255)
                                                                                                      cGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCg                                                >  1:304387/1‑150 (MQ=255)
                                                                                                           ggcggcAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATg                                           >  2:1545784/1‑150 (MQ=255)
                                                                                                            gcggcAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGc                                          <  1:879218/150‑1 (MQ=255)
                                                                                                             cggcAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCt                                         >  1:916647/1‑150 (MQ=255)
                                                                                                                     aTGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTa                                 <  1:784666/150‑1 (MQ=255)
                                                                                                                       gATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAAc                               <  1:622301/150‑1 (MQ=255)
                                                                                                                         tCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACgg                             >  2:1087446/1‑150 (MQ=255)
                                                                                                                          cGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGc                            >  1:307684/1‑150 (MQ=255)
                                                                                                                              aaCTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTg                        <  1:1504280/150‑1 (MQ=255)
                                                                                                                               aCTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTgg                       <  2:1516022/150‑1 (MQ=255)
                                                                                                                               aCTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTgg                       <  1:1201305/150‑1 (MQ=255)
                                                                                                                               aCTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTgg                       <  1:883230/150‑1 (MQ=255)
                                                                                                                                cTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTggg                      <  2:377657/150‑1 (MQ=255)
                                                                                                                                cTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTggg                      <  2:1494651/150‑1 (MQ=255)
                                                                                                                                 tGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGa                     <  1:1640912/150‑1 (MQ=255)
                                                                                                                                      gCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTATAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATgcgc                >  1:1250817/1‑150 (MQ=255)
                                                                                                                                       ccAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCa               >  1:1669353/1‑150 (MQ=255)
                                                                                                                                       ccAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACATGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCa               >  1:1269044/1‑150 (MQ=255)
                                                                                                                                          atttACGCAGTTTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAg             >  2:1307793/3‑149 (MQ=255)
                                                                                                                                          atttACGCAGTTTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAg             <  1:1307793/147‑1 (MQ=255)
                                                                                                                                            ttACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTg          >  1:803669/1‑150 (MQ=255)
                                                                                                                                              aCGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGa        >  2:465116/1‑150 (MQ=255)
                                                                                                                                               cGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATTTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGaa       <  2:307684/150‑1 (MQ=255)
                                                                                                                                                 cAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACg     >  1:863434/1‑150 (MQ=255)
                                                                                                                                                 cAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACg     >  2:1079906/1‑150 (MQ=255)
                                                                                                                                                 cAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACg     >  1:42037/1‑150 (MQ=255)
                                                                                                                                                 cAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACg     >  2:252144/1‑150 (MQ=255)
                                                                                                                                                    aCTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCtt  <  2:42037/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCGAATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAACTATGATCGCAAACTGTTCGCCAAATTACGCAAACTGCGTAAATCGATAGCCGATGAAAGTAATGTCCCGCCGTACGTGGTGTTTAACGACGCAACCTTGATTGAGATGGCTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTT  >  NC_000913/4007338‑4007635

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