Predicted mutation
evidence seq id position mutation annotation gene description
RA AM406671 1,681,037 T→C E332G (GAA→GGA)  pepV ← Dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM4066711,681,0370TC99.8% 4425.3 / NA 1062E332G (GAA→GGA) pepVDipeptidase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  major base C (527/531);  minor base G (1/1);  total (530/532)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 1074 total aligned reads displayed.

GTTCAACTTCAACAACCCCATCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTTCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACAAAGTAAGATAAGTTGCTCCATT  >  AM406671/1680891‑1681182
                                                                                                                                                  |                                                                                                                                                 
gTTCAACTTCAACAACCCCATCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTccc                                                                                                                                                 <  2:8793871/149‑1 (MQ=255)
  tcaacttcaacAACCCCATCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCAt                                                                                                                                               <  1:387109/149‑1 (MQ=255)
  tcaacttcaacAACCCCATCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCAt                                                                                                                                               >  1:2032673/1‑149 (MQ=255)
     acttcaacAACCCCATCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCaa                                                                                                                                            >  2:10088861/1‑149 (MQ=255)
     acttcaacAACCCCATCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCaa                                                                                                                                            <  1:8882077/149‑1 (MQ=255)
       ttccacAACCCCATCAAGTTTAGCAAGAATTTCTTGCATACGCTCTTGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCCAATGACCAAACACCAGCATTCATTTAGGTATTTCCCATCAAtt                                                                                                                                          <  2:65602/149‑1 (MQ=255)
       ttcaacAACCCCATCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAAtt                                                                                                                                          <  2:10527389/149‑1 (MQ=255)
        ccaacatCCCCATCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGATGTTCCCTCGTGGGAAACGAAAATTGAGGGCGATGTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATtc                                                                                                                                         <  2:7836525/148‑1 (MQ=255)
         caacaaCCCCATCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATtca                                                                                                                                        >  1:7624073/1‑149 (MQ=255)
           acaacCCCATCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATtcatc                                                                                                                                      <  1:8139806/149‑1 (MQ=255)
            cacccccATCAAGTTTAGCCAGCATTTTTTGCATCCGCTCTGAGCTGTTCCCTTGTGGGAAACGGAAATTGAGTGCGATTTTACCTTCCCCATTTTCATCAAATGACCAATCACCAGCATTCATTGAGGTATTTCGCATCAATtcatca                                                                                                                                     <  2:9359860/146‑1 (MQ=255)
              acCCCATCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATcaac                                                                                                                                   <  1:861308/149‑1 (MQ=255)
               ccccATCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATcaaca                                                                                                                                  >  1:11072603/1‑149 (MQ=255)
               ccccATCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATcaaca                                                                                                                                  >  2:6730824/1‑149 (MQ=255)
                cccATCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATcaacaa                                                                                                                                 >  2:7902110/1‑149 (MQ=255)
                 ccATCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACaaa                                                                                                                                <  1:1317329/149‑1 (MQ=255)
                 ccATCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACaaa                                                                                                                                >  1:10836400/1‑149 (MQ=255)
                   aTCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAg                                                                                                                              <  1:7787716/149‑1 (MQ=255)
                     cAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCt                                                                                                                            >  1:9312840/1‑149 (MQ=255)
                         tttAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTAc                                                                                                                        >  1:6908193/1‑149 (MQ=255)
                          ttAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTAcc                                                                                                                       <  2:2286173/149‑1 (MQ=255)
                           tAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCt                                                                                                                      >  1:11273729/1‑149 (MQ=255)
                           tAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCt                                                                                                                      >  2:8093804/1‑149 (MQ=255)
                            aGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTa                                                                                                                     >  2:8487734/1‑149 (MQ=255)
                               aaGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAAtt                                                                                                                  <  1:6126046/149‑1 (MQ=255)
                                 gAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAAtttt                                                                                                                >  1:4517198/1‑149 (MQ=255)
                                    tttCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCa                                                                                                             <  1:2210830/149‑1 (MQ=255)
                                     ttCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCAc                                                                                                            <  1:5386447/149‑1 (MQ=255)
                                        ttGCATACGCTCTGGGCTGTTCCCTCGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCtt                                                                                                         >  2:6274228/1‑149 (MQ=255)
                                             tACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGa                                                                                                    >  2:1709513/1‑149 (MQ=255)
                                             tACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGa                                                                                                    <  1:10558345/149‑1 (MQ=255)
                                              aCGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGAt                                                                                                   >  1:9123904/1‑149 (MQ=255)
                                                  tctGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTc                                                                                               >  2:1422431/1‑149 (MQ=255)
                                                     gggCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCaag                                                                                            <  1:4728688/149‑1 (MQ=255)
                                                     gggCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCaag                                                                                            <  1:519319/149‑1 (MQ=255)
                                                      ggCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCaaga                                                                                           <  1:1640224/149‑1 (MQ=255)
                                                       gCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCaagaa                                                                                          <  1:4301897/149‑1 (MQ=255)
                                                         tGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGt                                                                                        >  1:1208644/1‑149 (MQ=255)
                                                         tGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGt                                                                                        >  1:9623168/1‑149 (MQ=255)
                                                           ttCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGttt                                                                                      >  2:10687945/1‑149 (MQ=255)
                                                             cccTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGttttt                                                                                    <  2:2967266/149‑1 (MQ=255)
                                                              ccTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTc                                                                                   >  1:7088565/1‑149 (MQ=255)
                                                              ccTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTc                                                                                   >  1:7919176/1‑149 (MQ=255)
                                                              ccTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTCTTCACCTTCATGATCTTCAAGAAGTTTTTc                                                                                   >  1:5569316/1‑149 (MQ=255)
                                                                  gtgGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGcc                                                                               <  2:8932587/149‑1 (MQ=255)
                                                                   tgGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGccc                                                                              <  2:865718/149‑1 (MQ=255)
                                                                      gAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCaa                                                                           >  1:9924251/1‑149 (MQ=255)
                                                                        aaCGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACt                                                                         <  2:4265150/149‑1 (MQ=255)
                                                                        aaCGGAAATTGAGGGCGATTTTACCTACACCATTTTCATCAAATGACCAAACCCCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATNATCTTCAAGAAGTTTTTCAGCCCCAACt                                                                         <  1:6358963/149‑1 (MQ=255)
                                                                          cGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACttt                                                                       <  2:5743078/149‑1 (MQ=255)
                                                                           ggAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTa                                                                      >  1:1851146/1‑149 (MQ=255)
                                                                             aaaTTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAAt                                                                    >  2:7742988/1‑149 (MQ=255)
                                                                                  gAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAgg                                                               <  1:6553046/149‑1 (MQ=255)
                                                                                   aGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGgc                                                              <  2:10861842/149‑1 (MQ=255)
                                                                                    gggCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGgca                                                             >  2:4709986/1‑149 (MQ=255)
                                                                                     ggCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGgcag                                                            <  2:6067017/149‑1 (MQ=255)
                                                                                     ggCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGgcag                                                            >  1:681000/1‑149 (MQ=255)
                                                                                     ggCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTCTTTCAGCCCCAACTTTAATGAAGgcag                                                            >  2:2804925/1‑149 (MQ=255)
                                                                                       cGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGgcagca                                                          >  2:3993276/1‑149 (MQ=255)
                                                                                       cGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGgcagca                                                          <  1:3750877/149‑1 (MQ=255)
                                                                                       cGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGgcagca                                                          >  1:5798309/1‑149 (MQ=255)
                                                                                         aTTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGgcagcagc                                                        >  1:10360662/1‑149 (MQ=255)
                                                                                         aTTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGgcagcagc                                                        >  1:4065557/1‑149 (MQ=255)
                                                                                         aTTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGgcagcagc                                                        >  1:48246/1‑149 (MQ=255)
                                                                                            ttACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCAcc                                                     >  2:10358631/1‑149 (MQ=255)
                                                                                             tACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCg                                                    >  1:3234267/1‑149 (MQ=255)
                                                                                                cTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCa                                                 >  2:5144135/1‑149 (MQ=255)
                                                                                                 ttCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAg                                                >  2:3199524/1‑149 (MQ=255)
                                                                                                    aCCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCaa                                             <  2:4941801/149‑1 (MQ=255)
                                                                                                     ccATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCTTGATCTTCAAGAAGTTTTTCAGCACCAACTTTAATGAAGGCAGCAGCACCGTCAGCaaa                                            >  2:4439026/1‑149 (MQ=255)
                                                                                                     ccATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCaaa                                            >  2:9646868/1‑149 (MQ=255)
                                                                                                      cATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAg                                           >  1:7255627/1‑149 (MQ=255)
                                                                                                       aTTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGt                                          >  1:10015130/1‑149 (MQ=255)
                                                                                                       aTTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGt                                          <  1:2505280/149‑1 (MQ=255)
                                                                                                          ttCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCaa                                       >  2:2627226/1‑149 (MQ=255)
                                                                                                           tcatcaAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCaaa                                      <  2:6913015/149‑1 (MQ=255)
                                                                                                            catcaAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAAt                                     >  2:11071067/1‑149 (MQ=255)
                                                                                                             atcaAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAAtt                                    <  1:2720175/149‑1 (MQ=255)
                                                                                                               caAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGa                                  >  2:1131365/1‑149 (MQ=255)
                                                                                                               caAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGa                                  >  2:5219216/1‑149 (MQ=255)
                                                                                                                  aTGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCACCAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTc                               >  2:5844555/1‑149 (MQ=255)
                                                                                                                      ccAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGa                           >  2:356069/1‑149 (MQ=255)
                                                                                                                      ccAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGa                           >  2:1263296/1‑149 (MQ=255)
                                                                                                                        aaaCACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAAc                         >  1:8447584/1‑149 (MQ=255)
                                                                                                                          acacCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACaa                       <  2:2499303/149‑1 (MQ=255)
                                                                                                                           cacCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACaaa                      >  2:7544091/1‑149 (MQ=255)
                                                                                                                           cacCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTACATGATCATCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACaaa                      >  2:7158863/1‑149 (MQ=255)
                                                                                                                              cAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACAAAGTa                   >  2:9196152/1‑149 (MQ=255)
                                                                                                                              cAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGGACCTAATATTTCACCTTCATGATCTTCAAGAAGTTTTTCAGGCCCAACTTTAATGAAGGGAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACAAAgca                   >  2:8269272/1‑147 (MQ=255)
                                                                                                                               aGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACAAAGTaa                  <  1:10172402/149‑1 (MQ=255)
                                                                                                                                  attcattGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACAAAGTAAGAt               >  1:3478204/1‑149 (MQ=255)
                                                                                                                                   ttcattGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCTCTTTTATTATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACAAAGTAAGATa              >  2:6418508/1‑149 (MQ=255)
                                                                                                                                    tcattGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACAAAGTAAGATaa             <  1:2938240/149‑1 (MQ=255)
                                                                                                                                    tcattGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGAGCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACAAAGTAAGATaa             >  1:5940336/1‑149 (MQ=255)
                                                                                                                                      attGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACAAAGTAAGATACGt           >  2:1937909/1‑149 (MQ=255)
                                                                                                                                      attGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACAAAGTAAGATAAGt           >  2:3822054/1‑149 (MQ=255)
                                                                                                                                       ttGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACAAAGTAAGATAAGtt          <  1:7419220/149‑1 (MQ=255)
                                                                                                                                       ttGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACAAAGTAAGATAAGtt          <  2:5372668/149‑1 (MQ=255)
                                                                                                                                         gAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACAAAGTAAGATAAGTTGc        >  1:829514/1‑149 (MQ=255)
                                                                                                                                              aTTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACAAAGTAAGATAAGTTGCTCCAt   >  1:5043484/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GTTCAACTTCAACAACCCCATCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTTCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACAAAGTAAGATAAGTTGCTCCATT  >  AM406671/1680891‑1681182

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