Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 834,352 A→T 100% Q428L (CAG→CTG)  dinG → ATP‑dependent DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913834,3520AT97.8% 294.7 / ‑2.1 90Q428L (CAG→CTG) dinGATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base T (33/55);  total (33/57)
Fisher's exact test for biased strand distribution p-value = 5.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

TTCGAGGCGCAAAGCAAACTCTGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTGTCACCTCCG  >  NC_000913/834255‑834447
                                                                                                 |                                                                                               
ttCGAGGCGCAAAGCAAACTCTGGCGGCTGGCTTCGCTGGCGCAATCTTCTGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGc                                                                                               <  2:1240526/100‑1 (MQ=255)
 tCGAGGCGCAAAGCAAACTCTGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGgctgct                                                                                              <  2:1841336/100‑1 (MQ=255)
 tCGAGGCGCAAAGCAAACTCTGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCt                                                                                              <  1:352471/100‑1 (MQ=255)
  cGAGGCGCAAAGCAAACTCTGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAAGGCTGCTa                                                                                             >  2:1139255/1‑100 (MQ=255)
  cGAGGCGCAAAGCAAACGCTGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTa                                                                                             <  2:2221117/100‑1 (MQ=255)
   gAGGCGCAAAGCAAACTCTGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTac                                                                                            <  1:2063273/100‑1 (MQ=255)
      gcgcAAAGCAAACTCTGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACAcc                                                                                         <  1:182791/100‑1 (MQ=255)
      gcgcAAAGCAAACTCTGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACAcc                                                                                         <  2:1258480/100‑1 (MQ=255)
      gcgcAAAGCAAACTCTGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACAcc                                                                                         <  2:1784733/100‑1 (MQ=255)
        ggaaagcaaaCTCTGGCGGCTGGCGTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAACGGGCGACGCGGGAAGAGCGCGAAGTGCTGCTACACctc                                                                                       <  2:1582275/98‑1 (MQ=255)
           aaccaaaCTCTGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTgg                                                                                    <  2:714507/97‑1 (MQ=255)
            agcaaaCTCTGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGt                                                                                   >  1:207195/1‑100 (MQ=255)
            agcaaaCTCTGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGGGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGt                                                                                   <  2:363313/100‑1 (MQ=255)
             gcaaaCTCTGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGtt                                                                                  >  1:586279/1‑100 (MQ=255)
             gcaaaCTCTGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGtt                                                                                  >  1:1928517/1‑100 (MQ=255)
             gcaaaCTCTGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGtt                                                                                  <  2:486075/100‑1 (MQ=255)
             gcaaaCTCTGGCGGCTGGCTTCGCTCGCGCAATCTTCCTGTGCACCGGTGACCAAATGGGCGACGCGGGCAGACAGCGAATGGCTGCTACACGTCTGGtt                                                                                  >  2:2436140/1‑100 (MQ=255)
                aaCTCTGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCa                                                                               >  2:120385/1‑100 (MQ=255)
                 cctctGTCGGCTGGCTTCGCTGGCGCACTCTTCCGGTGCACCGGTGACCAACTGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACCCCTCTGGTTTCAc                                                                              >  2:247334/2‑100 (MQ=255)
                  cactGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTCTCACt                                                                             >  1:1727816/3‑100 (MQ=255)
                     tGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCg                                                                          <  2:2359047/100‑1 (MQ=255)
                       gcggcTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTg                                                                        <  2:2152626/100‑1 (MQ=255)
                       gcggcTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTg                                                                        <  2:765825/100‑1 (MQ=255)
                       gcggcTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTg                                                                        <  2:1710454/100‑1 (MQ=255)
                       gcggcTGGCTCCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTg                                                                        <  2:1681761/100‑1 (MQ=255)
                         ggctggctTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTggg                                                                      <  2:1847693/100‑1 (MQ=255)
                         ggctggctTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTgg                                                                       <  1:605183/99‑1 (MQ=255)
                           ctggctTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGaa                                                                    <  2:1550746/100‑1 (MQ=255)
                             ggctTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAAt                                                                   >  2:2207166/1‑99 (MQ=255)
                              gctTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATAc                                                                 <  1:261448/100‑1 (MQ=255)
                               ctTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACg                                                                >  1:342170/1‑100 (MQ=255)
                               ctTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACg                                                                <  1:480845/100‑1 (MQ=255)
                               ctTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACg                                                                <  1:62446/100‑1 (MQ=255)
                                ttCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACgt                                                               >  1:615362/1‑100 (MQ=255)
                                 tCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACgtg                                                              <  1:779654/100‑1 (MQ=255)
                                 tCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACgtg                                                              <  1:1192637/100‑1 (MQ=255)
                                  cGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACgtgt                                                             >  1:1530504/1‑100 (MQ=255)
                                  cGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACgtgt                                                             >  1:2203743/1‑100 (MQ=255)
                                   gCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTc                                                            <  2:1845145/100‑1 (MQ=255)
                                   gCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTc                                                            <  2:406857/100‑1 (MQ=255)
                                    cTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCa                                                           >  1:123407/1‑100 (MQ=255)
                                      ggCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGc                                                         <  2:769921/100‑1 (MQ=255)
                                       gcgcAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCg                                                        >  1:340709/1‑100 (MQ=255)
                                        cgcAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGa                                                       >  1:1703650/1‑100 (MQ=255)
                                        cgcAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGa                                                       >  1:278507/1‑100 (MQ=255)
                                        cgcAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGa                                                       <  1:1448453/100‑1 (MQ=255)
                                          cAATCTGCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATc                                                     <  2:1454589/100‑1 (MQ=255)
                                           aaTCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCa                                                    >  1:221796/1‑100 (MQ=255)
                                               ttCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTg                                                <  2:1748948/100‑1 (MQ=255)
                                                 ccGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGa                                              <  2:223657/100‑1 (MQ=255)
                                                 ccGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGa                                              <  1:1889695/100‑1 (MQ=255)
                                                   ggTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGaaa                                            >  1:642164/1‑100 (MQ=255)
                                                   ggTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGa                                              <  2:1782865/98‑1 (MQ=255)
                                                        aCCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGgctg                                       <  2:664027/100‑1 (MQ=255)
                                                         ccGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGgctgc                                      <  1:2126999/100‑1 (MQ=255)
                                                         ccGGTGACCAAATGGGCGACGCGGGAAGAGCACGAAGAGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGgctgc                                      <  2:1726019/100‑1 (MQ=255)
                                                           ggTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGctgctg                                    >  1:1900678/1‑100 (MQ=255)
                                                           ggTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGctgctg                                    >  1:1583772/1‑100 (MQ=255)
                                                              gACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTgg                                 <  1:2448722/100‑1 (MQ=255)
                                                              gACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTgg                                 <  2:1121043/100‑1 (MQ=255)
                                                               aCCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGgc                                <  1:2433215/100‑1 (MQ=255)
                                                               aCCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGgc                                >  2:2531962/1‑100 (MQ=255)
                                                                ccAAATTGGCGACGCGGGAAGGGCGCGAAGGGGTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGgcg                               <  2:1094210/100‑1 (MQ=255)
                                                                ccAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGgcg                               >  1:584734/1‑100 (MQ=255)
                                                                 cAAATGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGgcgc                              >  2:108635/1‑100 (MQ=255)
                                                                  aaaTGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCGGCGATCAGCTGGAAAGGCTGCTGTGGCGCa                             >  2:256804/1‑100 (MQ=255)
                                                                  aaaTGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCa                             <  1:878031/100‑1 (MQ=255)
                                                                   aaTGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCa                             <  2:1842270/99‑1 (MQ=255)
                                                                   aaTGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAg                            <  2:1428254/100‑1 (MQ=255)
                                                                     tGGGCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTa                          >  1:1137652/1‑100 (MQ=255)
                                                                       ggCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTAtt                        >  2:2519332/1‑100 (MQ=255)
                                                                       ggCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTAtt                        >  1:1012551/1‑100 (MQ=255)
                                                                        gCGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTc                       <  2:950250/100‑1 (MQ=255)
                                                                         cGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTcc                      <  2:2079446/100‑1 (MQ=255)
                                                                         cGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTcc                      >  1:445797/1‑100 (MQ=255)
                                                                         cGACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTcc                      <  1:1944775/100‑1 (MQ=255)
                                                                          gACGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCg                     <  1:1613593/100‑1 (MQ=255)
                                                                           aCGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGc                    <  1:475386/100‑1 (MQ=255)
                                                                           aCGCGGGAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGc                    <  1:1623752/100‑1 (MQ=255)
                                                                             gcgGGAAGAGCGCGAAGGGCTGCTACACTTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGcac                  <  1:604672/100‑1 (MQ=255)
                                                                               gggAAGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACat                >  2:1923255/1‑100 (MQ=255)
                                                                                  aaGAGCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACattat             <  2:650995/100‑1 (MQ=255)
                                                                                   agagCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGACAGGCTGCTGTGGCGCAGTATTCCGCACattatt            >  2:2478529/1‑100 (MQ=255)
                                                                                    gagCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTg           <  2:1080364/100‑1 (MQ=255)
                                                                                     agCGCGAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTGt          >  1:413355/1‑100 (MQ=255)
                                                                                         cgAAGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTGTCAcc      >  1:833153/1‑100 (MQ=255)
                                                                                          gAAGGGCTGCTACACCTCTGGTTTTACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTGTCAcc      <  1:1282788/99‑1 (MQ=255)
                                                                                           aaGGGCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTGTCACCTc    <  1:602929/100‑1 (MQ=255)
                                                                                             gggCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTGTCACCTCCg  <  1:1074084/100‑1 (MQ=255)
                                                                                             gggCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTGTCACCTCCg  <  1:1336055/100‑1 (MQ=255)
                                                                                             gggCTGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTGTCACCTCCg  <  2:1940602/100‑1 (MQ=255)
                                                                                                 |                                                                                               
TTCGAGGCGCAAAGCAAACTCTGGCGGCTGGCTTCGCTGGCGCAATCTTCCGGTGCACCGGTGACCAAATGGGCGACGCGGGAAGAGCGCGAAGGGCAGCTACACCTCTGGTTTCACTGCGTGGGAATACGTGTCAGCGATCAGCTGGAAAGGCTGCTGTGGCGCAGTATTCCGCACATTATTGTCACCTCCG  >  NC_000913/834255‑834447

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