Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,649,460 T→A 77.7% K65* (AAA→TAA)  hdeB ← acid‑resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,649,4600TA77.7% 513.1 / 143.0 240K65* (AAA→TAA) hdeBacid‑resistance protein
Reads supporting (aligned to +/- strand):  ref base T (26/27);  new base A (95/90);  total (123/117)
Fisher's exact test for biased strand distribution p-value = 8.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 250 total aligned reads displayed.

TACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCAGTGACAAAGCCGCTACA  >  CP009273/3649313‑3649608
                                                                                                                                                   |                                                                                                                                                    
tACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTaat                                                                                                                                                    <  2:4119363/150‑1 (MQ=255)
   ttGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTaatata                                                                                                                                                 >  2:4478263/1‑150 (MQ=255)
    tGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCTGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATAc                                                                                                                                                >  1:2496037/1‑150 (MQ=255)
      cACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTg                                                                                                                                              >  1:3933153/1‑150 (MQ=255)
       aCCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGt                                                                                                                                             >  2:939243/1‑150 (MQ=255)
         cTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGttt                                                                                                                                           >  2:2401506/1‑150 (MQ=255)
           cATTTATTCGGCAAGTCATTAGATGCTTTATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTttct                                                                                                                                         <  1:4822991/150‑1 (MQ=255)
            aTTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTttctt                                                                                                                                        >  1:4783942/1‑150 (MQ=255)
            aTTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTCATATACTGTttctt                                                                                                                                        >  2:3015394/1‑150 (MQ=255)
             ttAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTttcttc                                                                                                                                       <  1:558353/150‑1 (MQ=255)
                  tCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATgca                                                                                                                                  >  1:2894311/1‑150 (MQ=255)
                  tCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATgca                                                                                                                                  >  1:4652911/1‑150 (MQ=255)
                   cGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATgcag                                                                                                                                 >  1:513426/1‑150 (MQ=255)
                    ggCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTATCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATgcagc                                                                                                                                >  2:1648762/1‑150 (MQ=255)
                     gCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATgcagca                                                                                                                               >  1:2268072/1‑150 (MQ=255)
                      cAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATgcagcag                                                                                                                              >  2:2867064/1‑149 (MQ=255)
                      cAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCAt                                                                                                                              >  2:53458/1‑150 (MQ=255)
                       aaGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATc                                                                                                                             <  1:4012847/150‑1 (MQ=255)
                        aGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATcc                                                                                                                            <  1:4245546/150‑1 (MQ=255)
                               aGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGc                                                                                                                     >  1:1518603/1‑150 (MQ=255)
                                   gCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAAcc                                                                                                                 <  2:2300013/150‑1 (MQ=255)
                                   gCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAAcc                                                                                                                 <  2:2813586/150‑1 (MQ=255)
                                   gCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAAcc                                                                                                                 <  1:3213894/150‑1 (MQ=255)
                                   gCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAAcc                                                                                                                 >  2:1087565/1‑150 (MQ=255)
                                   gCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAAcc                                                                                                                 <  1:4591857/150‑1 (MQ=255)
                                    cTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCg                                                                                                                <  2:3811028/150‑1 (MQ=255)
                                     ttGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCgg                                                                                                               <  2:3784962/150‑1 (MQ=255)
                                     ttGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCgg                                                                                                               <  1:202697/150‑1 (MQ=255)
                                      tGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCggg                                                                                                              >  1:835653/1‑150 (MQ=255)
                                      tGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCggg                                                                                                              >  2:2099100/1‑150 (MQ=255)
                                      tGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCggg                                                                                                              >  2:2553904/1‑150 (MQ=255)
                                            ttGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCAtt                                                                                                        >  1:1199411/1‑150 (MQ=255)
                                            ttGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCAtt                                                                                                        <  2:4035884/150‑1 (MQ=255)
                                            ttGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCAtt                                                                                                        <  1:4140640/150‑1 (MQ=255)
                                              gAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGc                                                                                                      >  1:5057042/1‑150 (MQ=255)
                                               aaGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCt                                                                                                     >  1:2771140/1‑150 (MQ=255)
                                               aaGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCt                                                                                                     >  2:3745009/1‑150 (MQ=255)
                                               aaGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCt                                                                                                     >  1:4759038/1‑150 (MQ=255)
                                                aGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCtt                                                                                                    >  2:2634354/1‑150 (MQ=255)
                                                aGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCtt                                                                                                    >  2:3326415/1‑150 (MQ=255)
                                                  gTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCtttt                                                                                                  <  1:852770/150‑1 (MQ=255)
                                                   tataCAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTg                                                                                                 >  1:764488/1‑150 (MQ=255)
                                                    ataCAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTgg                                                                                                <  2:504431/150‑1 (MQ=255)
                                                       cAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGAtt                                                                                             <  1:1258600/150‑1 (MQ=255)
                                                        aaaTTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTc                                                                                            <  1:338544/150‑1 (MQ=255)
                                                         aaTTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATAATGTTTCTTCATGAAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCa                                                                                           >  1:1309570/1‑150 (MQ=255)
                                                          aTTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAg                                                                                          >  1:1630647/1‑150 (MQ=255)
                                                          aTTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAg                                                                                          >  2:802572/1‑150 (MQ=255)
                                                          aTTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGAGTCAg                                                                                          >  2:1166733/1‑150 (MQ=255)
                                                           tttttCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGa                                                                                         <  2:13411/150‑1 (MQ=255)
                                                           tttttCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGa                                                                                         <  2:3895644/150‑1 (MQ=255)
                                                              ttCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCa                                                                                      >  1:3853810/1‑150 (MQ=255)
                                                                     ggTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAAtt                                                                               <  1:4047008/150‑1 (MQ=255)
                                                                        ttttCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCt                                                                            <  2:5141035/150‑1 (MQ=255)
                                                                          ttCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTgg                                                                          <  1:356908/150‑1 (MQ=255)
                                                                          ttCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTgg                                                                          >  2:1839258/1‑150 (MQ=255)
                                                                              tACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAgg                                                                      <  2:1679494/150‑1 (MQ=255)
                                                                               aCAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGt                                                                     <  1:3457373/150‑1 (MQ=255)
                                                                                     ttCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATc                                                               <  1:1836921/150‑1 (MQ=255)
                                                                                      tCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCt                                                              <  2:3183103/150‑1 (MQ=255)
                                                                                       cGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCtt                                                             >  2:1905034/1‑150 (MQ=255)
                                                                                       cGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCtt                                                             >  1:5018184/1‑150 (MQ=255)
                                                                                          tCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAg                                                          >  2:2988093/1‑150 (MQ=255)
                                                                                                tAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGAtt                                                    <  2:4533956/150‑1 (MQ=255)
                                                                                                    aaTTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCAtt                                                <  2:2924476/150‑1 (MQ=255)
                                                                                                     aTTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTg                                               <  2:2276347/150‑1 (MQ=255)
                                                                                                     aTTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTg                                               >  1:3320056/1‑150 (MQ=255)
                                                                                                           gtgagATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGcc                                         >  1:100999/1‑150 (MQ=255)
                                                                                                            tgagATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCa                                        >  2:3502906/1‑150 (MQ=255)
                                                                                                             gagaTCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCaa                                       >  1:2949288/1‑150 (MQ=255)
                                                                                                             gagaTCGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCaa                                       <  2:2425993/150‑1 (MQ=255)
                                                                                                                  cGGTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAg                                  >  2:3639915/1‑150 (MQ=255)
                                                                                                                   ggTTTCATTTAAAGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGa                                 >  1:2295042/1‑150 (MQ=255)
                                                                                                                          tttAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCg                          <  1:3270545/150‑1 (MQ=255)
                                                                                                                          tttAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCg                          <  2:1569777/150‑1 (MQ=255)
                                                                                                                           ttAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGt                         <  2:3876054/150‑1 (MQ=255)
                                                                                                                           ttAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGATGCCAACGCAGATTGTGCGt                         >  1:4407450/1‑150 (MQ=255)
                                                                                                                             aaaGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTc                       >  2:489044/1‑150 (MQ=255)
                                                                                                                              aaGTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCa                      >  1:2080291/1‑150 (MQ=255)
                                                                                                                                gTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCAcc                    >  2:579631/1‑150 (MQ=255)
                                                                                                                                gTAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCAcc                    >  1:327941/1‑150 (MQ=255)
                                                                                                                                 tAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCa                   <  1:3897846/150‑1 (MQ=255)
                                                                                                                                 tAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCa                   >  1:1044499/1‑150 (MQ=255)
                                                                                                                                 tAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCa                   <  1:1149404/150‑1 (MQ=255)
                                                                                                                                 tAACGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCa                   >  1:565885/1‑150 (MQ=255)
                                                                                                                                    cGGTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCAGTg                >  1:3807245/1‑150 (MQ=255)
                                                                                                                                     ggTATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCAGTGa               >  1:4265180/1‑150 (MQ=255)
                                                                                                                                       tATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCAGTGACa             <  2:2780200/150‑1 (MQ=255)
                                                                                                                                       tATCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCAGTGACa             <  2:4652911/150‑1 (MQ=255)
                                                                                                                                        aTCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCAGTGACaa            >  2:1428558/1‑150 (MQ=255)
                                                                                                                                         tCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACTCAGATTGTGCGTTCACCAGTGACaaa           <  1:1336149/150‑1 (MQ=255)
                                                                                                                                         tCGCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCAGTGACaaa           >  2:955113/1‑150 (MQ=255)
                                                                                                                                          cGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCAGTGACAAAg          >  1:2572867/1‑150 (MQ=255)
                                                                                                                                          cGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCAGTGACAAAg          >  1:1997839/1‑150 (MQ=255)
                                                                                                                                           gCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCAGTGACAAAGc         >  2:666280/1‑150 (MQ=255)
                                                                                                                                           gCCACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCAGTGACAAAGc         >  2:3252244/1‑150 (MQ=255)
                                                                                                                                            ccACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCAGTGACAAAGcc        <  2:1518603/150‑1 (MQ=255)
                                                                                                                                            ccACCTTAATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCAGTGACAAAGcc        >  1:3628644/1‑150 (MQ=255)
                                                                                                                                                  ttATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCAGTGACAAAGCCGCTACa  <  2:1986556/150‑1 (MQ=255)
                                                                                                                                                  ttATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCAGTGACAAAGCCGCTACa  <  1:3702857/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACGGTATCGCCACCTTTATATACTGTTTCTTCATGCAGCATCCACCATGCAACCGGGGTCATTGCTTTTGGATTCAGATCAATAAATTCCTGGCAGGTCATATCTTTAGCGGATTCATTGGCTGCCAACGCAGATTGTGCGTTCACCAGTGACAAAGCCGCTACA  >  CP009273/3649313‑3649608

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