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.7% 3842.1 / NA 927E332G (GAA→GGA) pepVDipeptidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (424/500);  minor base A (3/0);  total (427/500)
Fisher's exact test for biased strand distribution p-value = 9.74e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 935 total aligned reads displayed.

GTTCAACTTCAACAACCCCATCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTTCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACAAAGTAAGATAAGTTGCTCCATT  >  AM406671/1680891‑1681182
                                                                                                                                                  |                                                                                                                                                 
gTTCAACTTCAACAACCCCATCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTccc                                                                                                                                                 >  2:3330925/1‑149 (MQ=255)
gTTCAACTTCAACAACCCCATCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTccc                                                                                                                                                 >  2:577692/1‑149 (MQ=255)
  tcaacttcaacAACCCCATCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCAt                                                                                                                                               >  1:6134605/1‑149 (MQ=255)
      cttcaacAACCCCATCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAAt                                                                                                                                           >  1:1891221/1‑149 (MQ=255)
       ttcaacAACCCCATCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAAtt                                                                                                                                          <  1:7617121/149‑1 (MQ=255)
        tcaacAACCCCATCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATtc                                                                                                                                         >  1:4668446/1‑149 (MQ=255)
        tcaacAACCCCATCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATtc                                                                                                                                         <  2:2875907/149‑1 (MQ=255)
           acaacCCCATCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATtcatc                                                                                                                                      <  2:9437979/149‑1 (MQ=255)
           acaacCCCATCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATtcatc                                                                                                                                      >  1:4067683/1‑149 (MQ=255)
              acCCCATCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATcaac                                                                                                                                   >  1:7198619/1‑149 (MQ=255)
              acCCCATCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATcaac                                                                                                                                   >  2:1409658/1‑149 (MQ=255)
               ccccATCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATcaaca                                                                                                                                  >  2:3918122/1‑149 (MQ=255)
                 ccATCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACaaa                                                                                                                                >  1:7399158/1‑149 (MQ=255)
                  cATCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACaaaa                                                                                                                               <  2:2314818/149‑1 (MQ=255)
                   aTCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAg                                                                                                                              <  1:999139/149‑1 (MQ=255)
                    tCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGc                                                                                                                             >  1:2284315/1‑149 (MQ=255)
                    tCAAGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGc                                                                                                                             >  2:184708/1‑149 (MQ=255)
                       aGTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGt                                                                                                                          <  1:1670474/149‑1 (MQ=255)
                        gTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTa                                                                                                                         <  1:3398412/149‑1 (MQ=255)
                        gTTTAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTa                                                                                                                         >  1:3496165/1‑149 (MQ=255)
                         tttAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTAc                                                                                                                        <  2:7745626/149‑1 (MQ=255)
                          ttAGCAAGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTAcc                                                                                                                       <  1:7445660/149‑1 (MQ=255)
                               aaGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAAtt                                                                                                                  >  2:4134884/1‑149 (MQ=255)
                                aGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAAttt                                                                                                                 <  2:5109166/149‑1 (MQ=255)
                                aGAATTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAAttt                                                                                                                 <  1:8674261/149‑1 (MQ=255)
                                aGAATTTCGTGCAGACGCTCTGGGCGTTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCACATGACCAAACACCAGGATTCATTGAGGTATTTCCCGTCAATTCATCAACAAAAGCTGTACCTAAttt                                                                                                                 <  2:455641/149‑1 (MQ=255)
                                  aaTTTCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAAttttt                                                                                                               >  2:5898801/1‑149 (MQ=255)
                                     ttCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCAc                                                                                                            <  1:3241672/149‑1 (MQ=255)
                                      tCTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCAcc                                                                                                           <  1:7813620/149‑1 (MQ=255)
                                       tttGCATTCGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCt                                                                                                          <  2:6470535/148‑1 (MQ=255)
                                       cTTTCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCAGCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCt                                                                                                          <  2:1961851/149‑1 (MQ=255)
                                       cTTGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCt                                                                                                          >  1:8943609/1‑149 (MQ=255)
                                         tGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTc                                                                                                        <  1:5670124/149‑1 (MQ=255)
                                         tGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTc                                                                                                        >  2:75627/1‑149 (MQ=255)
                                         tGCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTc                                                                                                        <  1:2428282/149‑1 (MQ=255)
                                          gCATACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTTAGAGCGATTTTACCTTCACCGTTTTCATCAAACGACCAAACACCAGCATTCATTGAGGTATTCCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCa                                                                                                       <  2:6933105/149‑1 (MQ=255)
                                             tACGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGa                                                                                                    >  2:2654512/1‑149 (MQ=255)
                                              aCGCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGAt                                                                                                   >  1:1129812/1‑149 (MQ=255)
                                                gCTCTGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCt                                                                                                 >  2:6272597/1‑149 (MQ=255)
                                                 ctctGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCtt                                                                                                <  1:6301240/149‑1 (MQ=255)
                                                 ctctGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCtt                                                                                                <  2:5596686/149‑1 (MQ=255)
                                                  tctGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTc                                                                                               <  1:2046574/149‑1 (MQ=255)
                                                   ctGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCa                                                                                              <  1:4362752/149‑1 (MQ=255)
                                                    tGGGCTGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCaa                                                                                             >  2:9118182/1‑149 (MQ=255)
                                                    tGGGCTGCTCCCTTGTGGGATACGGAAATTGAGGGCGATGTTACCTTCCCCATGTTCATCAAATGACCGAACACCAGCATTCACTCCGGCATTTCCCATCAATTCATCAACAAAAACTGTACCTAATTTTTCACCTTCATGATCTTCaa                                                                                             >  2:2496818/1‑149 (MQ=255)
                                                         tGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGt                                                                                        >  2:8801470/1‑149 (MQ=255)
                                                         tGTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGt                                                                                        <  2:3722394/149‑1 (MQ=255)
                                                          gTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGtt                                                                                       <  1:3132361/149‑1 (MQ=255)
                                                          gTTCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGtt                                                                                       >  1:8171135/1‑149 (MQ=255)
                                                           ttCCCTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGttt                                                                                      >  2:6742693/1‑149 (MQ=255)
                                                           ttCCCTTGTGGGAAACGGAAATTGAGGGCAATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAAGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGttt                                                                                      <  1:150226/149‑1 (MQ=255)
                                                             cccTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGttttt                                                                                    <  1:4540025/149‑1 (MQ=255)
                                                              ccTTGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTc                                                                                   >  2:547892/1‑149 (MQ=255)
                                                                ttGTGGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAg                                                                                 <  2:4665695/149‑1 (MQ=255)
                                                                 tgtgGGAAACGGAAATTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGc                                                                                <  1:1244487/149‑1 (MQ=255)
                                                                       aaaCGGAAATTGAGGGCGATTTTGCCTTCACCATTTTCATCAAATGACCGAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGCACGTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAAc                                                                          >  2:5271474/1‑149 (MQ=255)
                                                                              aaTTGAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATg                                                                   >  1:8804328/1‑149 (MQ=255)
                                                                                  gAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAgg                                                               >  1:5479792/1‑149 (MQ=255)
                                                                                  gAGGGCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAgg                                                               <  2:8236820/149‑1 (MQ=255)
                                                                                     ggCGATTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCCTCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCACCTTTAATGAAGgcag                                                            >  2:4864739/1‑149 (MQ=255)
                                                                                         aTTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGgcagcagc                                                        <  2:7467771/149‑1 (MQ=255)
                                                                                         aTTTTACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGgcagcagc                                                        <  2:7919110/149‑1 (MQ=255)
                                                                                            ttACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCAcc                                                     <  2:1551630/149‑1 (MQ=255)
                                                                                            ttACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCAcc                                                     <  2:7304844/149‑1 (MQ=255)
                                                                                             tACCTTCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCg                                                    >  1:6400630/1‑149 (MQ=255)
                                                                                                 ttCACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAg                                                <  2:845199/149‑1 (MQ=255)
                                                                                                   cACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCa                                              >  2:6534773/1‑149 (MQ=255)
                                                                                                   cACCATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCa                                              <  1:5057819/149‑1 (MQ=255)
                                                                                                     ccATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCaaa                                            >  1:2981602/1‑149 (MQ=255)
                                                                                                     ccATTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCaaa                                            <  1:6814056/149‑1 (MQ=255)
                                                                                                       aTTTTCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGt                                          >  2:9208349/1‑149 (MQ=255)
                                                                                                         tttCATCAAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCa                                        <  1:1340916/149‑1 (MQ=255)
                                                                                                           tcatcaAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGACGGCAGCAGCACCGTCAGCAAAGTCaaa                                      <  2:4326982/149‑1 (MQ=255)
                                                                                                           tcatcaAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCaaa                                      <  1:4215795/149‑1 (MQ=255)
                                                                                                           tcatcaAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCaaa                                      >  1:8336950/1‑149 (MQ=255)
                                                                                                           tcatcaAATGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCaaa                                      <  1:9118857/149‑1 (MQ=255)
                                                                                                             atgaaaTGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAAtt                                    <  2:3174712/146‑1 (MQ=255)
                                                                                                                  aTGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTc                               >  2:3005965/1‑149 (MQ=255)
                                                                                                                  aTGACCAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTc                               >  1:5168724/1‑149 (MQ=255)
                                                                                                                       cAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGaa                          <  2:1131731/149‑1 (MQ=255)
                                                                                                                       cAAACACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGaa                          <  1:7027890/149‑1 (MQ=255)
                                                                                                                        aaaCACCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAAc                         <  2:6134605/149‑1 (MQ=255)
                                                                                                                          acacCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACaa                       <  2:3592046/149‑1 (MQ=255)
                                                                                                                          acacCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACaa                       >  1:6608926/1‑149 (MQ=255)
                                                                                                                          acacCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACaa                       <  1:5391550/149‑1 (MQ=255)
                                                                                                                            acCAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACAAAg                     <  1:5925872/149‑1 (MQ=255)
                                                                                                                             ccAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACaaa                      >  2:2195145/1‑147 (MQ=255)
                                                                                                                              cAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACAAAGTa                   <  2:8928342/149‑1 (MQ=255)
                                                                                                                              cAGCATTCATTGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACAAAGTa                   >  2:7064158/1‑149 (MQ=255)
                                                                                                                                    tcattGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACAAATTAAGATaa             >  2:8700209/1‑149 (MQ=255)
                                                                                                                                    tcattGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACAAAGTAAGATaa             <  2:1271684/149‑1 (MQ=255)
                                                                                                                                    tcattGAGGGATTTCCCATCAATGAATCAACAAAAGCTGTGCCTAATTTTCCACCTTCATGATCTTCAAGAAGTGTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACAACGTAAGATaa             >  2:3431371/1‑149 (MQ=255)
                                                                                                                                     cattGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACAAAGTAAGATAAg            <  2:3500138/149‑1 (MQ=255)
                                                                                                                                     cattGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACAAAGTAAGATAAg            <  1:3604791/149‑1 (MQ=255)
                                                                                                                                      attGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACAAAGTAAGATAAGt           <  1:9427552/149‑1 (MQ=255)
                                                                                                                                      attGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACAAAGTAAGATAAGt           >  1:2636686/1‑149 (MQ=255)
                                                                                                                                       ttGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACAAAGTAAGATAAGtt          >  1:3736971/1‑149 (MQ=255)
                                                                                                                                        tGAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACAAAGTAAGATAAGTTg         >  1:1000697/1‑149 (MQ=255)
                                                                                                                                         gAGGTATTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGAAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACAAAGTAAGATAAGTTGc        <  1:488159/149‑1 (MQ=255)
                                                                                                                                              aTTTCCCATCAATTCATCAACAAAAGCTGTACCTAATTTTTCACCTTCATGATCTTCAAGAAGTTTTTCAGCCCCAACTTTAATGTAGGCAGCAGCACCGTCAGCAAAGTCAAATTGATTCAAGAACAAAGTAAGATAAGTTGCTCCAt   <  1:5780052/149‑1 (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: