Predicted mutation
evidence seq id position mutation annotation gene description
RA AM406671 1,935,405 (T)7→6 pseudogene (1748/2210 nt) comEC ← pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM4066711,935,3990T.96.6% 4214.9 / 120.8 947pseudogene (1754/2210 nt)comECpseudogene
Reads supporting (aligned to +/- strand):  ref base T (30/2);  new base . (461/452);  total (493/454)
Fisher's exact test for biased strand distribution p-value = 2.73e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 1031 total aligned reads displayed.

GAAATTTTGTATGGTAAAAGTTTCCATAAGTCACTATGGAATCATTATTATCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAGGATAAGAAATCACCAACATGGTCAGCATCTGTATGAGTCAGAATAAGTTGGTCAATATGA  >  AM406671/1935251‑1935543
                                                                                                                                                    |                                                                                                                                                
gAAATTTTGTATGGTAAAAGTTTCCATAAGTCACTATGGAATCATTATTATCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCt                                                                                                                                                  <  2:3399131/149‑1 (MQ=255)
     tttGTATGGTAAAAGTTTCCATAAGTCACTATGGAATCAGTATTATCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCtttttt                                                                                                                                             <  1:6976062/149‑1 (MQ=255)
      ttGTATGGTAAAAGTTTCCATAAGTCACTATGGAATCATTATTATCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCtttttta                                                                                                                                            <  1:6889589/149‑2 (MQ=255)
       tGTATGGTAAAAGTTTCCATAAGTCACTATGGAATCATTATTATCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCttttttaa                                                                                                                                           <  2:2051719/149‑3 (MQ=255)
         tATGGTAAAAGTTTCCATAAGTCACTATGGAATCATTATTATCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCttttttaatt                                                                                                                                         <  1:1386105/149‑5 (MQ=255)
            ggTAAAAGTTTCCATAAGTCACTATGGAATCATTATTATCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAAttttt                                                                                                                                     <  1:1214923/149‑1 (MQ=255)
            ggTAAAAGTTTCCATAAGTCACTATGGAATCATTATTATCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAAttttt                                                                                                                                     <  1:2448272/149‑1 (MQ=255)
             gTAAAAGTTTCCATAAGTCACTATGGAATCATTATTATCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTc                                                                                                                                    <  1:9405163/149‑1 (MQ=255)
              tAAAAGTTTCCATAAGTCACTATGGAATCATTATTATCATCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTcc                                                                                                                                   <  2:33000/149‑1 (MQ=255)
               aaaaCTTTCCATAAGTCACTATGGAATCATTATTATCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAAAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCca                                                                                                                                  <  2:5863172/149‑1 (MQ=255)
                aaaGTTTCCATAAGTCACTATGGAATCATTATTATCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAAGCGGAATTTTTGCTTTTTTAATTTTTCcac                                                                                                                                 <  1:1689832/149‑1 (MQ=255)
                  aGTTTCCATAAGTCACTATGGAATCATTATTATCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACaa                                                                                                                               <  2:4258519/149‑1 (MQ=255)
                  aGTTTCCATAAGTCACTATGGAATCATTATTATCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACaa                                                                                                                               <  1:3415554/149‑1 (MQ=255)
                   gTTTCCATAAGTCACTATGGAATCATTATTATCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACaaa                                                                                                                              <  2:4071499/149‑1 (MQ=255)
                     ttCCATAAGTCACTATGGAATCATTATTATCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAac                                                                                                                            <  2:644616/149‑1 (MQ=255)
                     ttCCATAAGTCACTATGGAATCATTATTATCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAac                                                                                                                            >  2:8946643/1‑149 (MQ=255)
                       ccATAAGTCACTATGGAATCATTATTATCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAacta                                                                                                                          >  2:2897334/1‑149 (MQ=255)
                        cATAAGTCACTATGGAATCATTATTATCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAactac                                                                                                                         >  2:6082065/1‑149 (MQ=255)
                          tAAGTCACTATGGAATCATTATTATCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAAAAAAAAAAAAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTg                                                                                                                       <  2:9455064/149‑1 (MQ=255)
                           aaGTCACTATGGAATCATTATTATCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTGCTAACATGAATCGGAATTTTTGCTTTTTTAAATTTTCCACAAAACTACTGt                                                                                                                      >  1:4684074/1‑149 (MQ=255)
                           aaGTCACTATGGAATCATTATTATCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGt                                                                                                                      >  2:9273136/1‑149 (MQ=255)
                            aGTCACTATGGAATCATTATTATCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGTTTTTTTATTTTTTCCACAAAACTACTGtt                                                                                                                     >  2:9087482/1‑149 (MQ=255)
                               cACTATGGAATCATTATTATCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGt                                                                                                                  <  1:4190207/149‑1 (MQ=255)
                                  tATGGAATCATTATTATCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTaa                                                                                                               <  1:6681640/149‑1 (MQ=255)
                                       aaTCATTATTATCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATtctc                                                                                                          <  1:7368689/149‑1 (MQ=255)
                                        ttCATTATTATCACCTTTTTTAGTGTAACTATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTcc                                                                                                         <  2:5099898/148‑1 (MQ=255)
                                        aTCATTATTATCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTcc                                                                                                         >  1:650806/1‑149 (MQ=255)
                                         tCATTATTATCACCTTTTCTAGTGTAACCATTTGATAACACTTGCAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATATTTCCAAAAAACTACTGTTGGTTAATTCTccc                                                                                                        >  2:739141/1‑149 (MQ=255)
                                            ttattaTCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGa                                                                                                     <  2:6663093/149‑1 (MQ=255)
                                             tattatCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGAc                                                                                                    >  1:3923969/1‑149 (MQ=255)
                                             tattatCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGAc                                                                                                    >  1:351229/1‑149 (MQ=255)
                                               ttatCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTg                                                                                                  >  1:7498877/1‑149 (MQ=255)
                                               ttatCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTg                                                                                                  <  2:4416619/149‑1 (MQ=255)
                                                tatCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGa                                                                                                 >  1:8754122/1‑149 (MQ=255)
                                                tatCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGa                                                                                                 >  1:858137/1‑149 (MQ=255)
                                                tatCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGa                                                                                                 <  1:5093655/149‑1 (MQ=255)
                                                tatCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTAACAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGa                                                                                                 <  2:4720161/149‑1 (MQ=255)
                                                  tCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGAcc                                                                                               <  1:6412762/149‑1 (MQ=255)
                                                  tCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGAcc                                                                                               >  2:2740958/1‑149 (MQ=255)
                                                     ccTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAg                                                                                            <  1:3225740/149‑1 (MQ=255)
                                                      cTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGa                                                                                           <  2:3009856/149‑1 (MQ=255)
                                                       ttttttAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGAt                                                                                          <  1:846053/149‑1 (MQ=255)
                                                       ttttttAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGAt                                                                                          <  2:173043/149‑1 (MQ=255)
                                                         ttttAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGAttt                                                                                        >  1:3336205/1‑149 (MQ=255)
                                                          tttAGTGTAACCATTTGACAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTCAATTTTTCCACAAAACTACTGTTGATTAATTCTCCCGGACTGACCCAGATTtc                                                                                       <  2:5603834/149‑1 (MQ=255)
                                                           ttAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAAATTtct                                                                                      >  1:9236235/1‑149 (MQ=255)
                                                            tAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTtctc                                                                                     >  2:2527780/1‑149 (MQ=255)
                                                             aGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTctct                                                                                    <  2:7939117/149‑1 (MQ=255)
                                                             aGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTctct                                                                                    >  1:9558622/1‑149 (MQ=255)
                                                             aGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTctct                                                                                    <  2:872769/149‑1 (MQ=255)
                                                               tgtAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGCCCCAGATTTCTCtta                                                                                  >  2:319381/1‑149 (MQ=255)
                                                                  aaCCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCttattt                                                                               >  1:2680769/1‑149 (MQ=255)
                                                                    ccATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCttattttt                                                                             <  1:548648/149‑1 (MQ=255)
                                                                     cATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCttattttta                                                                            <  1:5217945/149‑1 (MQ=255)
                                                                            tAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATc                                                                     >  2:4957609/1‑149 (MQ=255)
                                                                            tAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATc                                                                     >  2:2278909/1‑149 (MQ=255)
                                                                             aaaaCTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCa                                                                    <  1:301092/149‑1 (MQ=255)
                                                                              aaaCTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAg                                                                   <  2:8325019/149‑1 (MQ=255)
                                                                                aCTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCt                                                                 <  1:8260023/149‑1 (MQ=255)
                                                                                 cTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCGACTTTACTAACATGAATCGGAATTTTCGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTa                                                                <  2:5365815/149‑1 (MQ=255)
                                                                                 cTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTa                                                                >  1:7081538/1‑149 (MQ=255)
                                                                                  ttGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTaa                                                               <  1:4345485/149‑1 (MQ=255)
                                                                                   tGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGCATTTTTGCTTTTTTACTTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAg                                                              >  1:1036137/1‑149 (MQ=255)
                                                                                    gaagaaAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAgg                                                             >  2:85691/1‑149 (MQ=255)
                                                                                     aagaaAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAGGa                                                            >  1:1796727/1‑149 (MQ=255)
                                                                                     aagaaAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGTTTTTTTTTTTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAGGa                                                            >  2:1534772/1‑149 (MQ=255)
                                                                                       gaaAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAGGATa                                                          >  1:4476743/1‑149 (MQ=255)
                                                                                       gaaAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAGGATa                                                          >  2:525205/1‑149 (MQ=255)
                                                                                       aaaaaCTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAGGATa                                                          <  2:7711400/148‑1 (MQ=255)
                                                                                         aaaCTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAGGATAAg                                                        <  1:2231344/149‑1 (MQ=255)
                                                                                          aaCTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAGGATAAGa                                                       >  1:7781376/1‑149 (MQ=255)
                                                                                          aaCTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAGGATAAGa                                                       >  1:7226853/1‑149 (MQ=255)
                                                                                             tATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAGGATAAGAAAt                                                    <  2:1319951/149‑1 (MQ=255)
                                                                                             tATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAGGATAAGAAAt                                                    <  1:1523427/149‑1 (MQ=255)
                                                                                                cAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAGGATAAGAAATCAc                                                 >  1:4891160/1‑149 (MQ=255)
                                                                                                  aaaaTCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAGGATAAGAAATCACCa                                               >  2:8179150/1‑149 (MQ=255)
                                                                                                    aaTCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAGGATAAGAAATCACCAAc                                             >  2:4634403/1‑149 (MQ=255)
                                                                                                     aTCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAGGATAAGAAATCACCAACa                                            >  1:6023358/1‑149 (MQ=255)
                                                                                                     aTCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAGGATAAGAAATCACCAACa                                            <  2:351229/149‑1 (MQ=255)
                                                                                                           aTTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAGGATAAGAAATCACCAACATGGTCa                                      >  1:808825/1‑149 (MQ=255)
                                                                                                            ttttATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAGGATAAGAAATCACCAACATGGTCAg                                     <  1:1961483/149‑1 (MQ=255)
                                                                                                              ttATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAGGATAAGAAATCACCAACATGGTCAGCa                                   >  2:1158267/1‑149 (MQ=255)
                                                                                                               tATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCTAGATTTCTCTTATTTTTATTTTATCAGCTAAGGATAAGAAATCACCAACATGGTCAGCAt                                  <  1:5670046/149‑1 (MQ=255)
                                                                                                               tATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAGGATAAGAAATCACCAACATGGTCAGCAt                                  <  1:38327/149‑1 (MQ=255)
                                                                                                               tATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAGGATAAGAAATCACCAACATGGTCAGCAt                                  <  1:3122835/149‑1 (MQ=255)
                                                                                                                aTCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAGGATAAGAAATCACCAACATGGTCAGCATc                                 >  1:8966564/1‑149 (MQ=255)
                                                                                                                 tCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAGGATAAGAAATCACCAACATGGTCAGCATCt                                <  1:623972/149‑1 (MQ=255)
                                                                                                                   cccAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAGGATAAGAAATCACCAACATGGTCAGCATCTGt                              <  2:570405/149‑1 (MQ=255)
                                                                                                                    ccAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTTCACAAAACTACTGTTGGTAAATTCTCCCGAACTGACCCAGAATTCTCTTATTTTTATTTAATTAGCTAAGGATAAGAAATCACCAACATGGTCAGCATCTgca                             >  2:6945687/1‑147 (MQ=255)
                                                                                                                          ttACTAACATGAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAGGATAAGAAATCACCAACATGGTCAGCATCTGTATGAGTc                       <  2:6819790/149‑1 (MQ=255)
                                                                                                                            aCTAACATGAATCGGAATTTTTGCTTTTTTTTTTTTTTTACAAAACTACTGTTGGTTAATTTTCCTGGTCTGACCTAGATTTCTCTTATTTTTATTTTATCAGCTTAGGATAAGAAATCACCAACATGGTCTGCATCTGTATGTTTCAg                     >  2:3710857/1‑149 (MQ=255)
                                                                                                                               aaCATGAATCGGAATTTTTGCTTTTTTAATTTTTCTACAAAACTACTGTTGGTTAATTCTCCAGGACTGACCCAGATTTCTCTTATTTTTATTTTATCTGCTATGGATACGAAATCACCAACATGGTGAGCATATGTATGTGTCAGAAt                  >  1:2850589/1‑149 (MQ=255)
                                                                                                                                aCATGAATCGGAATTTTTTCTTTTTTTTTTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTTTCTGCTAAGGATAAGAAATCACCATTATGTTCTGCATCTATATGAGTCAGATTa                 >  2:2097549/1‑149 (MQ=255)
                                                                                                                                    gAATCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAGGATAAGAAATCACCAACATGGTCAGCATCTGTATGAGTCAGAATAAGtt             <  2:7279446/149‑1 (MQ=255)
                                                                                                                                      aTCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAGGATAAGAAATCACCAACATGGTCAGCATCTGTATGAGTCAGAATAAGTTgg           <  2:1870784/149‑1 (MQ=255)
                                                                                                                                      aTCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAGGATAAGAAATCACCAACATGGTCAGCATCTGTATGAGTCAGAATAAGTTgg           <  2:3922926/149‑1 (MQ=255)
                                                                                                                                      aTCGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAGGATAAGAAATCACCAACATGGTCAGCATCTGTATGAGTCAGAATAAGTTgg           <  2:6293068/149‑1 (MQ=255)
                                                                                                                                        cGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAGGATAAGAAATCACCAACATGGTCAGCATCTGTATGAGTCAGAATAAGTTGGTc         >  2:1068874/1‑149 (MQ=255)
                                                                                                                                        cGGAATTTTTGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAGGATAAGAAATCACCAACATGGTCAGCATCTGTATGAGTCAGAATAAGTTGGTc         >  2:7005047/1‑149 (MQ=255)
                                                                                                                                             tttttGCTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAGGATAAGAAATCACCAACATGGTCAGCATCTGTATGAGTCAGAATAAGTTGGTCAatat    <  1:5817227/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GAAATTTTGTATGGTAAAAGTTTCCATAAGTCACTATGGAATCATTATTATCACCTTTTTTAGTGTAACCATTTGATAAAACTTGAAGAAAACTATCAAAAATCGGAATTTTATCCCCAACTTTACTAACATGAATCGGAATTTTTGCTTTTTTTAATTTTTCCACAAAACTACTGTTGGTTAATTCTCCCGGACTGACCCAGATTTCTCTTATTTTTATTTTATCAGCTAAGGATAAGAAATCACCAACATGGTCAGCATCTGTATGAGTCAGAATAAGTTGGTCAATATGA  >  AM406671/1935251‑1935543

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: