Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,626,422 T→A 6.4% Q320L (CAG→CTG)  dcp ← dipeptidyl carboxypeptidase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,626,4220TA 6.4% 260.7 / 4.2 91Q320L (CAG→CTG) dcpdipeptidyl carboxypeptidase II
Reads supporting (aligned to +/- strand):  ref base T (43/40);  new base A (5/2);  total (49/42)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

CTTCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAAGTAGTGTTGCCTGTT  >  NC_000913/1626276‑1626569
                                                                                                                                                  |                                                                                                                                                   
cTTCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCTCCTCATCAAGATCAAATTTCTCCCGCCGTACATGTTCCGCATAAAATGCCCAGTCCCCCGGCTGCGCGCTAAACCCGCCCTGCTGGTTATCGATAACCGCCTg                                                                                                                                                    >  6:376819/1‑148 (MQ=255)
cTTCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTg                                                                                                                                                    >  6:149520/1‑148 (MQ=255)
cTTCAATTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATAAAATTTCTCACGCCGTACATGTTACGCATAAAATGCCCAGTCCAACGGCTGCGCGCTAAACACGCCCTTCTGCTTATCGATAACCGCATg                                                                                                                                                    >  5:80574/1‑148 (MQ=255)
 ttCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTa                                                                                                                                                  >  5:141146/1‑149 (MQ=255)
         acacCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGc                                                                                                                                           <  5:149520/148‑1 (MQ=255)
            ccGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAAt                                                                                                                                       >  3:279145/1‑149 (MQ=255)
              gtgtTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAAtt                                                                                                                                      >  8:125026/1‑148 (MQ=255)
                     aTTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGc                                                                                                                                           >  7:163202/1‑136 (MQ=255)
                     aTTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGc                                                                                                                                           <  8:163202/136‑1 (MQ=255)
                     aTTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCAtcgctc                                                                                                                              >  7:93618/1‑149 (MQ=255)
                        cAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAGATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCa                                                                                                                           >  2:59743/1‑149 (MQ=255)
                        cAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCAGTATGGAGGCTAATTCATCGCTCGCa                                                                                                                           >  8:156053/1‑149 (MQ=255)
                             tatGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCAGTATGGAGGCTAAt                                                                                                                                       <  3:161485/132‑1 (MQ=255)
                             tatGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCAGTATGGAGGCTAAt                                                                                                                                       <  3:167419/132‑1 (MQ=255)
                             tatGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCAGTATGGAGGCTAAt                                                                                                                                       >  4:167419/1‑132 (MQ=255)
                             tatGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCAGTATGGAGGCTAAt                                                                                                                                       >  4:161485/1‑132 (MQ=255)
                                   ttGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCCCGCTAAACCCCCCCTGCTGCTTATCGATAACCGCCTGAATGGAGGCTAAATCAACGCACGCACCTTGACgcgc                                                                                                                >  5:44059/1‑149 (MQ=255)
                                           cgcCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTgg                                                                                                          <  6:271070/147‑1 (MQ=255)
                                           cgcCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGaa                                                                                                        <  8:93618/149‑1 (MQ=255)
                                             ccTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAAc                                                                                                       >  6:124228/1‑148 (MQ=255)
                                              cTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCCCGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCAGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACaa                                                                                                     >  8:216656/1‑149 (MQ=255)
                                              cTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCAtcgctcg                                                                                                                             <  8:366947/125‑1 (MQ=255)
                                              cTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCAtcgctcg                                                                                                                             >  7:366947/1‑125 (MQ=255)
                                               tcatcaAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCATATCGATAACCGCCTGTATGGAGGCTAAATCATCGCTCGCACGTAGACGCGCCGCTGGAACAAt                                                                                                    >  3:357920/1‑149 (MQ=255)
                                                  tcaAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTTTGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTc                                                                                                 <  6:141146/149‑1 (MQ=255)
                                                    aaGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTcc                                                                                                <  7:156053/148‑1 (MQ=255)
                                                      gATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCg                                                                                              >  4:110374/1‑148 (MQ=255)
                                                          aaaTTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCCTCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATa                                                                                          >  5:241129/1‑148 (MQ=255)
                                                          aaaTTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATa                                                                                          >  8:197764/1‑148 (MQ=255)
                                                          aaaTTTCTCCCGCCGTACCTGTTCCGAATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAAACGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCaga                                                                                          >  5:235156/1‑146 (MQ=255)
                                                           aaTTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATaaa                                                                                        <  7:125026/149‑1 (MQ=255)
                                                               tctcCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCTCGTTGGCGCGCCGCTGGAACAATTTCCCGCATaa                                                                                         <  2:32166/144‑1 (MQ=255)
                                                               tctcCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGGCGCGCCGCTGGAACAATTTCCCGCATaa                                                                                         >  1:32166/1‑144 (MQ=255)
                                                               tctcCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGGCGCGCCGCTGGAACAATTTCCCGCATaa                                                                                         <  2:38577/144‑1 (MQ=255)
                                                               tctcCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGGCGCGCCGCTGGAACAATTTCCCGCATaa                                                                                         >  1:38577/1‑144 (MQ=255)
                                                                   cctccgTACCTGTTCCGCATAAAATGCCCAGTCCCTCGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTaa                                                                                   <  6:44059/143‑1 (MQ=255)
                                                                    cgccgTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAGCCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAAATCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAATtgc                                                                               >  7:56250/1‑149 (MQ=255)
                                                                      ccgTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGtgctg                                                                             <  3:110374/149‑1 (MQ=255)
                                                                        gTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGTAACAATTTCCCGCATAAAGTTAAGTgctgct                                                                           >  4:75479/1‑149 (MQ=255)
                                                                           ccTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTc                                                                         <  4:279145/148‑1 (MQ=255)
                                                                           ccTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCa                                                                        <  8:95276/149‑1 (MQ=255)
                                                                            cTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAg                                                                       <  8:354780/149‑1 (MQ=255)
                                                                             tGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAgg                                                                      >  6:184566/1‑149 (MQ=255)
                                                                              gTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGgt                                                                     >  7:964/1‑149 (MQ=255)
                                                                                   gCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTaaat                                                                                                                                      >  2:373413/1‑77 (MQ=255)
                                                                                   gCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAAtt                                                                                                                                      <  1:373413/79‑1 (MQ=255)
                                                                                   gCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGtgc                                                                               >  5:409283/1‑134 (MQ=255)
                                                                                   gCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGtgc                                                                               <  6:409283/134‑1 (MQ=255)
                                                                                    cATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGttttt                                                               <  5:184566/149‑1 (MQ=255)
                                                                                       aaaaTGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGc                                                             >  7:153975/1‑148 (MQ=255)
                                                                                              ccAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCTCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGAt                                                     >  6:368077/1‑149 (MQ=255)
                                                                                              ccAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGAt                                                     <  1:222167/149‑1 (MQ=255)
                                                                                                aGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTCGAACAATTTCCCGCATAAAGTTAAGTGCTGCATCAGGTGTTTTTGCCATCTGATCg                                                   >  1:107573/1‑149 (MQ=255)
                                                                                                 gTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCgg                                                  >  2:164654/1‑149 (MQ=255)
                                                                                                   cccACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCGTCAGGTGTTTTTGCCATCTGATCGGc                                                 >  3:225891/1‑148 (MQ=255)
                                                                                                     cACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGAt                                              >  1:421323/1‑149 (MQ=255)
                                                                                                      aCGGCTGCGCGCTAACCCCGCCCTGCTGCTTATCGATAACCGCCGGTATGGAGGCTAATTCATCGCTCGCCCGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGGGTGTTTGCCAGCTGATCGGCGAtt                                             >  7:434547/1‑149 (MQ=255)
                                                                                                      aCGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTaa                                                                                   <  5:382366/111‑1 (MQ=255)
                                                                                                      aCGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTaa                                                                                   >  6:382366/1‑111 (MQ=255)
                                                                                                      aCGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGAtt                                             <  6:77796/149‑1 (MQ=255)
                                                                                                      aCGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCAGTATGGAGTCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTTATGCTTCAGGTGTTTTTGACATCTTATCGACGAtt                                             >  4:22126/1‑149 (MQ=255)
                                                                                                        ggCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGAtttt                                           >  5:47825/1‑149 (MQ=255)
                                                                                                          cTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTc                                          >  4:237637/1‑148 (MQ=255)
                                                                                                             cgcgcTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCa                                        <  3:237637/147‑1 (MQ=255)
                                                                                                             cgcgcTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCa                                        <  7:189792/147‑1 (MQ=255)
                                                                                                                gcTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCgg                                   >  4:41969/1‑149 (MQ=255)
                                                                                                                 cTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGc                                  <  8:374449/149‑1 (MQ=255)
                                                                                                                 cTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGc                                  <  7:216656/149‑1 (MQ=255)
                                                                                                                  tAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCa                                 >  8:62235/1‑149 (MQ=255)
                                                                                                                     aCCCGACCTGCTGCTTATCGTTAACCGCCTGTATCTTGGGTTATTCATCGATCGCTCTTTTACGCTTCGCTGGAACAATTTCCCGCATAATGTTAAGTGCTGCTTCATGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCTTaa                              <  8:404139/149‑1 (MQ=255)
                                                                                                                      cccGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGGTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAAt                             <  1:59743/149‑1 (MQ=255)
                                                                                                                        cGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGa                           <  8:964/149‑1 (MQ=255)
                                                                                                                        cGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGa                           <  3:75479/149‑1 (MQ=255)
                                                                                                                          cccTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAgg                         <  2:421323/149‑1 (MQ=255)
                                                                                                                            ctgctgCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGa                        >  2:231035/1‑148 (MQ=255)
                                                                                                                             tgctgcTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGaa                       >  3:4095/1‑148 (MQ=255)
                                                                                                                              gctgctTATCGATAACCGCCTGTATGGAGGCTAATGCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCGTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGaaaa                     <  6:241129/149‑1 (MQ=255)
                                                                                                                               cttctTATCGATAACCGCCTGTATTGAGGCTAATTCATAGCTCGCACGTTTACGCGCCGCTGGAACAATTTCCCTCATAAAGTTAATTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGaaaa                     <  6:235156/145‑1 (MQ=255)
                                                                                                                               ctccctATCGATAACCCCCTGTATGGAGGCTAATTCATCGCTCGCACGTGGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAACTGCTGATTCACGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGaaaa                     <  6:121659/143‑1 (MQ=255)
                                                                                                                                  ctTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCaa                 <  5:376819/149‑1 (MQ=255)
                                                                                                                                     aTCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAag                >  7:369066/1‑147 (MQ=255)
                                                                                                                                     aTCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAag                >  1:186888/1‑147 (MQ=255)
                                                                                                                                      tCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAagta              <  2:186888/148‑1 (MQ=255)
                                                                                                                                        gATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTcc                                         <  4:38414/119‑1 (MQ=255)
                                                                                                                                         aTAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAAGTAgtgt          >  4:125022/1‑149 (MQ=255)
                                                                                                                                                ccTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTaa                                                                                   <  5:95992/69‑1 (MQ=255)
                                                                                                                                                ccTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTaa                                                                                   >  6:95992/1‑69 (MQ=255)
                                                                                                                                                ccTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAAGTAGTGTTGCCtgt   >  2:228361/1‑149 (MQ=255)
                                                                                                                                                 cTGTATGGAGGCTAATTCATCTCTAGCACGTTGACGCGCCGATGGAACAATTTCCAGCATAAAGTTAAGTGCTGCTTCTGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCTTAATGAGGAAAACCAAGTAGTGTTGCCtgtt  <  3:22126/149‑1 (MQ=255)
                                                                                                                                                 cTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAAGTAGTGTTGCCtgtt  <  1:13177/149‑1 (MQ=255)
                                                                                                                                                 cTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAAGTAGTGTTGCCtgtt  <  1:234483/149‑1 (MQ=255)
                                                                                                                                                  tGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAAGTAgtgt          >  5:31699/1‑140 (MQ=255)
                                                                                                                                                  tGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAAGTAgtgt          <  6:31699/140‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CTTCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGATAACCGCCTGTATGGAGGCTAATTCATCGCTCGCACGTTGACGCGCCGCTGGAACAATTTCCCGCATAAAGTTAAGTGCTGCTTCAGGTGTTTTTGCCATCTGATCGGCGATTTTCCATGCGGCATAATGAGGAAAACCAAGTAGTGTTGCCTGTT  >  NC_000913/1626276‑1626569

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