Predicted mutation
evidence seq id position mutation annotation gene description
RA AM406671 945,239 +C pseudogene (229/761 nt) tnp904 → pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM406671945,2391.C99.8% 1018.8 / NA 441pseudogene (229/761 nt)tnp904pseudogene
Reads supporting (aligned to +/- strand):  ref base . (1/0);  major base C (227/212);  minor base A (1/0);  total (229/212)
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 475 total aligned reads displayed.

GGGACGCTATGGTGCGGTTAGAATTACCAAGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGCAACAGCTCCTAATAAAGTATGGCTGGGAGACATGACCTATATCCCTACCAAAGAAGGCACCTTATACTTA  >  AM406671/945091‑945351
                                                                                                                                                      |                                                                                                                
gggACGCTATGGTGCGGTTAGAATTACCAAGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAg                                                                                                                   >  1:7842432/1‑149 (MQ=38)
  gACGCTATGGTGCGGTTAGAATTACCAAGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAgct                                                                                                                 <  1:7337079/149‑3 (MQ=37)
    cGCTATGGTGCGGTTAGAATTACCAAGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAgcttc                                                                                                               <  2:3592802/149‑5 (MQ=25)
     gCTATGGTGCGGTTAGAATTACCAAGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAgcttcg                                                                                                              >  1:1901225/1‑144 (MQ=25)
      cTATGGTGCGGTTAGAATTACCAAGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGc                                                                                                             >  1:1723324/1‑149 (MQ=21)
       tATGGTGCGGTTAGAATTACCAAGGTACTTCATAATACTGGTATTATAACCAACCCGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCt                                                                                                            >  1:1049676/1‑149 (MQ=12)
        aTGGTGCGGTTAGAATTACCAAGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCATGGGAAGCCGTTATAAATATAGAGATTATAACAGAAAAGGAGCTTCGCtt                                                                                                           >  2:1990275/1‑149 (MQ=11)
        aTGGTGCGGTTAGAATTACCAAGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCtt                                                                                                           <  1:4177609/149‑1 (MQ=21)
        aTGGTGCGGTTAGAATTACCAAGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCtt                                                                                                           <  1:2952194/149‑1 (MQ=21)
        aTGGTGCGGTTAGAATTACCAAGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCtt                                                                                                           <  1:769195/149‑1 (MQ=21)
         tGGTGCGGTTAGAATTACCAAGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCttt                                                                                                          <  1:456346/149‑1 (MQ=17)
         tGGTGCGGTTAGAATTACCAAGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCttt                                                                                                          >  2:6568752/1‑149 (MQ=18)
         tGGTGCGGTTAGAATTACCAAGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCttt                                                                                                          <  2:5877165/149‑1 (MQ=18)
         gggTGCGGTTAGAATTACCAAGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCttt                                                                                                          <  1:1569981/148‑1 (MQ=21)
          ggTGCGGTTAGAATTACCAAGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTc                                                                                                         >  1:5962296/1‑149 (MQ=21)
           gTGCGGTTAGAATTACCAAGGTACTTCATAATACTTGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCa                                                                                                        <  1:3681116/149‑1 (MQ=12)
           gTGCGGTTAGAATTACCAAGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCa                                                                                                        <  2:4122953/149‑1 (MQ=17)
            tGCGGTTAGAATTACCAAGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCaa                                                                                                       >  1:4821762/1‑149 (MQ=21)
             gCGGTTAGAATTACCAAGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAg                                                                                                      <  2:8840155/149‑1 (MQ=21)
               ggTTAGAATTACCAAGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTAGATAAGACTTTATGCCAAGGGAAGCCCTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGAc                                                                                                    <  2:5518903/149‑1 (MQ=11)
                gTTAGAATTACCAAGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGAcc                                                                                                   >  1:7225/1‑149 (MQ=21)
                 ttAGAATTACCAAGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGAccc                                                                                                  <  2:7705338/149‑1 (MQ=21)
                 ttAGAATTACCAAGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGAccc                                                                                                  >  2:7884131/1‑149 (MQ=21)
                 ttAGAATTACCAAGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGAccc                                                                                                  >  1:5574783/1‑149 (MQ=21)
                  tAGAATTACCAAGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCa                                                                                                 <  1:8153621/149‑1 (MQ=18)
                     aaTTACCAAGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAAtt                                                                                              <  2:3334588/149‑1 (MQ=21)
                     aaTTACCAAGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAAtt                                                                                              >  2:1564781/1‑149 (MQ=17)
                      aTTACCAAGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGAAGTTATGCCAAGGGGCGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGCCCCAAtct                                                                                             >  2:2239143/1‑147 (MQ=11)
                       ttACCAAGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATtta                                                                                            >  1:7599933/1‑149 (MQ=18)
                        tACCAAGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATttaa                                                                                           >  2:4820553/1‑149 (MQ=21)
                        tACCAAGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATttaa                                                                                           <  1:6591838/149‑1 (MQ=21)
                         aCCAAGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATttaat                                                                                          <  1:5348528/149‑1 (MQ=21)
                           cAAGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATttaatta                                                                                        >  1:8837863/1‑149 (MQ=21)
                            aaGGTACTTCATAATACTGGTATTATAACCAACACGCAACGTGTGGGGAAACTGATGCACTTGATGGGCCTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATttaattaa                                                                                       >  1:8765658/1‑149 (MQ=255)
                            aaGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATttaattaa                                                                                       <  2:7011119/149‑1 (MQ=18)
                             aGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTtaattaat                                                                                      <  2:1708816/149‑1 (MQ=21)
                             aGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATCtaattaat                                                                                      >  2:3456634/1‑149 (MQ=255)
                              ggTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGAccc                                                                                                  <  1:7918634/136‑1 (MQ=21)
                              ggTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATc                                                                                     <  1:966330/149‑1 (MQ=18)
                               gTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCa                                                                                    >  1:8043351/1‑149 (MQ=21)
                                tACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAg                                                                                   <  1:7035631/149‑1 (MQ=21)
                                tACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAg                                                                                   <  2:2316071/149‑1 (MQ=21)
                                tACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAg                                                                                   >  1:8265981/1‑149 (MQ=21)
                                 aCTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGa                                                                                  <  2:8629366/149‑1 (MQ=21)
                                  cTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGAt                                                                                 <  1:6387498/149‑1 (MQ=21)
                                  cTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGAt                                                                                 >  2:3018095/1‑149 (MQ=17)
                                  cTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGAt                                                                                 <  1:2528214/149‑1 (MQ=21)
                                   ttCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATc                                                                                >  2:4332625/1‑149 (MQ=18)
                                   atCATAATACTGGTATTATAACCACCACGTAACGTGTGGGGAAACTGAAGCACTTGATGGAGCTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACCGAACAGGACCTTTGCTTTAAACACCCACTTTACTTAATCAGCTc                                                                                >  2:5984430/2‑149 (MQ=11)
                                    tCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCt                                                                               <  1:1299683/149‑1 (MQ=17)
                                     cATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCtt                                                                              >  2:8516580/1‑149 (MQ=18)
                                      ataataCTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCttt                                                                             <  2:2842585/149‑1 (MQ=18)
                                      ataataATGGTATTATAAACAACACGAAACGTGGGGGGAAACTGATGCACTTGATGGGACTTTATGGCAAGGGAAGCCGTTATAAATATACACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCttt                                                                             >  2:4493634/1‑149 (MQ=11)
                                       taataCTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTa                                                                            <  2:7254491/149‑1 (MQ=21)
                                       taataCTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTa                                                                            <  2:5072336/149‑1 (MQ=17)
                                       taataCTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTa                                                                            >  1:226100/1‑149 (MQ=21)
                                        aataCTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTaa                                                                           <  1:7453537/149‑1 (MQ=21)
                                        aataCTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTaa                                                                           >  1:3466376/1‑149 (MQ=21)
                                        aataCTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTAGATGGGACTGTATGCGAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGTCCCAATTTAATTAATCCGATCTTTaa                                                                           <  2:8945278/149‑1 (MQ=11)
                                         ataCTGGTATTATAACCAACACGAAACGTGTGGGGAAAGTGTTTCACTTGATGGGACTATATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTaaa                                                                          >  2:4645524/1‑149 (MQ=11)
                                         ataCTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTaaa                                                                          <  2:2630294/149‑1 (MQ=17)
                                          taCTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGTCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAg                                                                         >  2:8689759/1‑149 (MQ=16)
                                          taCTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAg                                                                         >  2:3907828/1‑149 (MQ=21)
                                           aCTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGc                                                                        <  1:915029/149‑1 (MQ=18)
                                           aCTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGc                                                                        >  1:329906/1‑149 (MQ=21)
                                            cTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGCa                                                                       >  1:1219988/1‑149 (MQ=21)
                                            cTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGCa                                                                       >  1:1836055/1‑149 (MQ=21)
                                             tGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGCaa                                                                      <  2:4186186/149‑1 (MQ=21)
                                             tGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGCaa                                                                      <  1:4517110/149‑1 (MQ=21)
                                             tGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATAAGATCTTTAAAGCaa                                                                      >  2:1830654/1‑149 (MQ=12)
                                              ggTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGCAAc                                                                     >  1:6818171/1‑149 (MQ=21)
                                              ggTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGCAAc                                                                     >  1:3079294/1‑149 (MQ=18)
                                               gTATTATAACCAACACGATACGTGTGGGGTAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCACTTTAATTAATCAGATCTTTAAAGCAACa                                                                    <  2:5374734/149‑1 (MQ=11)
                                               gTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTAGAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGCAACa                                                                    >  2:5745014/1‑149 (MQ=255)
                                                tattatAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGCAACAg                                                                   <  2:9133677/149‑1 (MQ=21)
                                                 attatAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGCCACAGc                                                                  <  2:797027/149‑1 (MQ=16)
                                                 attatAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGCAACAGc                                                                  >  1:3848897/1‑149 (MQ=18)
                                                   tataACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGCAACAGCTc                                                                >  1:4361486/1‑149 (MQ=21)
                                                   tataACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGCAACAGCTc                                                                <  2:6708538/149‑1 (MQ=17)
                                                      aaCCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGCAACAGCTCCta                                                             >  1:8212756/1‑149 (MQ=14)
                                                        ccAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGCAACAGCTCCtaat                                                           >  1:2042508/1‑149 (MQ=11)
                                                         cAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATTACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGCAACAGCTCCtaata                                                          <  1:5763383/149‑1 (MQ=2)
                                                         cAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGCAACAGCTCCtaata                                                          >  2:8321420/1‑149 (MQ=14)
                                                          aaCACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGCAACAGCTCCtaataa                                                         >  1:218895/1‑149 (MQ=1)
                                                          aaCACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGCAACAGCTCCtaataa                                                         >  2:754157/1‑149 (MQ=14)
                                                            cacGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGCAACAGCTCCTAATAAAg                                                       >  2:2508942/1‑149 (MQ=14)
                                                              cGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCATATCTTTAAAGCAACAGCTCCTAATAAAGTa                                                     <  2:112130/149‑1 (MQ=0)
                                                               gAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGCAACAGCTCCTAATAAAGTAt                                                    >  1:6902644/1‑149 (MQ=11)
                                                               gAAACGTGTGATGACACTAATGCACTTGACTGGACTTTATGCCAAGGGAAGCCGTTATAAATATACTCATCCTAACAGAAAAGGAGCTTCGCTCGCGAGTCCCCATTTAACTAATCAGCTCTTGAAAGCAACACCTCCTAATAGAGtct                                                    <  2:4966042/149‑3 (MQ=0)
                                                                aaaCGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGCAACAGCTCCTAATAAAGTATg                                                   >  1:50159/1‑149 (MQ=11)
                                                                 aaCGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGCAACAGCTCCTAATAAAGTATgg                                                  >  2:6847648/1‑149 (MQ=14)
                                                                 aaCGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGCAACAGCTCCTAATAAAGTATgg                                                  >  2:3761350/1‑149 (MQ=14)
                                                                 aaCGTGAGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGAACTTTAAAGCAACAGCTCCTAATAAAGTATgg                                                  >  2:824596/1‑149 (MQ=2)
                                                                  aCGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGCAACAGCTCCTAATAAAGTATGGc                                                 >  2:8103074/1‑149 (MQ=11)
                                                                  aCGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGCAACAGCTCCTAATAAAGTATGGc                                                 <  1:2775512/149‑1 (MQ=11)
                                                                  aCGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATAGAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGCAACAGCTCCTAATAAAGTATGGc                                                 >  1:5390539/1‑149 (MQ=9)
                                                                   cGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGCAACAGCTCCTAATAAAGTATGGCt                                                >  2:6293108/1‑149 (MQ=11)
                                                                    gtgtGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGCAACAGCTCCTAATAAAGTATGGCTg                                               <  1:4636006/149‑1 (MQ=14)
                                                                    gtgtGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGCAACAGCTCCTAATAAAGTATGGCTg                                               <  2:927735/149‑1 (MQ=11)
                                                                     tgtgGGGAAACTGATGCACTTGATGGGACTTTATGCAAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGCTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGCAACAGCTCCTAATAAAGTATGGCTgg                                              >  2:6154163/1‑149 (MQ=9)
                                                                                                                                                      |                                                                                                                
GGGACGCTATGGTGCGGTTAGAATTACCAAGGTACTTCATAATACTGGTATTATAACCAACACGAAACGTGTGGGGAAACTGATGCACTTGATGGGACTTTATGCCAAGGGAAGCCGTTATAAATATAAACATTATAACAGAAAAGGAGTTCGCTTTCAAGACCCAATTTAATTAATCAGATCTTTAAAGCAACAGCTCCTAATAAAGTATGGCTGGGAGACATGACCTATATCCCTACCAAAGAAGGCACCTTATACTTA  >  AM406671/945091‑945351

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