Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,778,219 C→A 11.6% D53Y (GAC→TAC)  yibN ← putative rhodanese‑related sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,778,2190CA11.6% 948.6 / 66.2 266D53Y (GAC→TAC) yibNputative rhodanese‑related sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base C (134/101);  new base A (14/17);  total (148/118)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

Only 100 of 285 total aligned reads displayed.

ACCGTCTACCACGATAACCGGTTTGTCTTTGTGCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGATGACGGCCAACAAATTGCATAATTT  >  CP009273/3778070‑3778368
                                                                                                                                                     |                                                                                                                                                     
aCCGTCTACCACGATAACCGGTTTGTCTTTGTGCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGtc                                                                                                                                                       >  2:3278564/1‑150 (MQ=255)
   gTCTACCACGATAACCGGTTTGTCTTTGTGCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGtcttc                                                                                                                                                    <  1:4269817/150‑1 (MQ=255)
      tACCACGATAACCGGTTTGTCTTTGTGCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCttt                                                                                                                                                 <  2:1873305/150‑1 (MQ=255)
          aCGATAACCGGTTTGTCTTTGTGCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTttgttg                                                                                                                                             <  1:932523/150‑1 (MQ=255)
               aaCCGGTTTGTCTTTGTGCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAg                                                                                                                                        >  1:1196517/1‑150 (MQ=255)
                    gTTTGTCTTTGTGCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGAcgcg                                                                                                                                   >  1:1217602/1‑150 (MQ=255)
                    gTTTGTCTTTGTGCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGAcgcg                                                                                                                                   <  2:3508189/150‑1 (MQ=255)
                       tGTCTTTGTGCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAg                                                                                                                                <  1:1294188/150‑1 (MQ=255)
                       tGTCTTTGTGCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAg                                                                                                                                <  2:2554401/150‑1 (MQ=255)
                         tCTTTGTGCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCt                                                                                                                              >  2:2904441/1‑150 (MQ=255)
                           tttGTGCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTc                                                                                                                            <  1:4096369/150‑1 (MQ=255)
                            ttGTGCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTca                                                                                                                           <  2:1294219/150‑1 (MQ=255)
                               tgCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTcacca                                                                                                                        <  1:2353330/150‑1 (MQ=255)
                                  ttCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGa                                                                                                                     >  2:2101232/1‑150 (MQ=255)
                                   tctcAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAg                                                                                                                    >  2:2592310/1‑150 (MQ=255)
                                     tcAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTa                                                                                                                  >  1:2996921/1‑150 (MQ=255)
                                      cAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTaa                                                                                                                 <  2:2676697/150‑1 (MQ=255)
                                      cAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTaa                                                                                                                 >  2:244776/1‑150 (MQ=255)
                                       aaGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAAt                                                                                                                <  1:4534939/150‑1 (MQ=255)
                                        aGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATc                                                                                                               >  1:4523927/1‑150 (MQ=255)
                                        aGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTATTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATc                                                                                                               >  1:1514964/1‑150 (MQ=255)
                                         gCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCa                                                                                                              <  1:3863235/150‑1 (MQ=255)
                                          cTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCAc                                                                                                             <  2:1821503/150‑1 (MQ=255)
                                               cAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCa                                                                                                        >  2:1402918/1‑150 (MQ=255)
                                                aaCATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCAc                                                                                                       <  1:2479599/150‑1 (MQ=255)
                                                  cATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACtt                                                                                                     >  2:3570448/1‑150 (MQ=255)
                                                    ttgttgGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTc                                                                                                   <  2:992627/150‑1 (MQ=255)
                                                         ggCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGt                                                                                              >  2:4349176/1‑150 (MQ=255)
                                                         ggCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGt                                                                                              >  1:553380/1‑150 (MQ=255)
                                                         ggCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGt                                                                                              <  2:937377/150‑1 (MQ=255)
                                                          gCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTc                                                                                             >  1:3367716/1‑150 (MQ=255)
                                                          gCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTATTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTc                                                                                             <  2:4161604/150‑1 (MQ=255)
                                                          gCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTATTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTc                                                                                             >  1:3310434/1‑150 (MQ=255)
                                                           cTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCa                                                                                            <  2:3710969/150‑1 (MQ=255)
                                                             ttGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTc                                                                                                         <  1:1332093/135‑1 (MQ=255)
                                                                aTTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGctc                                                                                       >  2:580137/1‑150 (MQ=255)
                                                                aTTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGctc                                                                                       <  1:1040949/150‑1 (MQ=255)
                                                                aTTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGctc                                                                                       >  1:2130930/1‑150 (MQ=255)
                                                                 tttCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGctct                                                                                      >  2:1745891/1‑150 (MQ=255)
                                                                   tCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTa                                                                                    <  1:3118757/150‑1 (MQ=255)
                                                                      cTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAAc                                                                                 <  1:2076525/150‑1 (MQ=255)
                                                                      cTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAAc                                                                                 <  1:2424955/150‑1 (MQ=255)
                                                                       tCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACg                                                                                >  1:1215882/1‑150 (MQ=255)
                                                                       tCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACg                                                                                >  1:4218963/1‑150 (MQ=255)
                                                                       tCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACg                                                                                >  2:3419234/1‑150 (MQ=255)
                                                                       tCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACg                                                                                >  2:1095974/1‑150 (MQ=255)
                                                                         ggCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTa                                                                              >  2:930020/1‑150 (MQ=255)
                                                                         ggCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTa                                                                              >  2:1597340/1‑150 (MQ=255)
                                                                         ggCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTa                                                                              >  2:2349230/1‑150 (MQ=255)
                                                                         ggCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTa                                                                              >  1:2190420/1‑150 (MQ=255)
                                                                           cAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGt                                                                            <  1:633547/150‑1 (MQ=255)
                                                                             aCAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTca                                                                          >  1:2623527/1‑150 (MQ=255)
                                                                             aCAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTca                                                                          <  2:1766340/150‑1 (MQ=255)
                                                                              cAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTcac                                                                         <  2:310771/150‑1 (MQ=255)
                                                                               aGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTcaca                                                                        >  2:2211204/1‑150 (MQ=255)
                                                                               aGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTATTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTcaca                                                                        >  2:758279/1‑150 (MQ=255)
                                                                                gATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACaa                                                                       <  1:1015182/150‑1 (MQ=255)
                                                                                 aTTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTATTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAg                                                                      <  1:3803853/150‑1 (MQ=255)
                                                                                 aTTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTATTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAg                                                                      <  1:232147/150‑1 (MQ=255)
                                                                                    aaTAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAACCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAAc                                                                   >  2:877451/1‑150 (MQ=255)
                                                                                     aTAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAAcc                                                                  >  2:773864/1‑150 (MQ=255)
                                                                                           ccTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGccacc                                                            >  1:3784329/1‑150 (MQ=255)
                                                                                               cGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTa                                                        >  2:4006330/1‑150 (MQ=255)
                                                                                                gATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTaa                                                       <  2:4302343/150‑1 (MQ=255)
                                                                                                  tatGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAAcg                                                     >  2:595076/1‑150 (MQ=255)
                                                                                                   atGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAAcgc                                                    >  1:3074864/1‑150 (MQ=255)
                                                                                                      gCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGAt                                                 <  1:3472077/150‑1 (MQ=255)
                                                                                                       ccTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCATCTTTGTAGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGAAACGCGATc                                                >  2:3831726/1‑150 (MQ=255)
                                                                                                        cTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTATTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATcc                                               >  2:1943769/1‑150 (MQ=255)
                                                                                                         tttACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCa                                              <  1:2977350/150‑1 (MQ=255)
                                                                                                           tACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAgg                                            <  2:648358/150‑1 (MQ=255)
                                                                                                             cGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCg                                          >  1:2021213/1‑150 (MQ=255)
                                                                                                              ggAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGa                                         >  1:4035782/1‑150 (MQ=255)
                                                                                                                aaGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATa                                       <  2:1684693/150‑1 (MQ=255)
                                                                                                                 aGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATAc                                      >  2:170602/1‑150 (MQ=255)
                                                                                                                 aGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATAc                                      >  1:1397518/1‑150 (MQ=255)
                                                                                                                  gTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACt                                     >  2:2702107/1‑150 (MQ=255)
                                                                                                                    catcaCGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCa                                   >  2:3829259/1‑150 (MQ=255)
                                                                                                                     atcaCGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAg                                  >  2:1552176/1‑150 (MQ=255)
                                                                                                                     atcaCGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAg                                  <  1:1700510/150‑1 (MQ=255)
                                                                                                                     atcaCGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGAACCAGGCGATACTCAg                                  >  1:300946/1‑150 (MQ=255)
                                                                                                                      tcaCGCTGACGTAAATCCACAACCACAGCGTATTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGt                                 >  1:3031271/1‑150 (MQ=255)
                                                                                                                        aCGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTAt                               <  1:1941164/150‑1 (MQ=255)
                                                                                                                          gCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATgg                             <  2:4480282/150‑1 (MQ=255)
                                                                                                                             gACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGAt                          <  1:3197911/150‑1 (MQ=255)
                                                                                                                             gACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGAt                          >  2:3207466/1‑150 (MQ=255)
                                                                                                                             gACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGAt                          <  1:1782935/150‑1 (MQ=255)
                                                                                                                              aCGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGATg                         <  1:2258315/150‑1 (MQ=255)
                                                                                                                               cGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGATGa                        >  2:3887414/1‑150 (MQ=255)
                                                                                                                                gTAAATCCACAACCACAGCGTATTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGATGAc                       <  1:3558591/150‑1 (MQ=255)
                                                                                                                                   aaTCCACAACCACAGCGTATTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGATGACGGc                    >  2:84847/1‑150 (MQ=255)
                                                                                                                                    aTCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGATGACGGcc                   >  1:4604270/1‑150 (MQ=255)
                                                                                                                                    aTCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGATGACGGcc                   >  2:1563906/1‑150 (MQ=255)
                                                                                                                                      ccACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGATGACGGCcaa                 <  1:3278564/150‑1 (MQ=255)
                                                                                                                                       cacaACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACATACCCCCCCCCCCCCCCCCCCGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGATGACGGCcaac                <  1:2601289/150‑1 (MQ=255)
                                                                                                                                         caACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGATGACGGCcaacaa              <  2:3888768/150‑1 (MQ=255)
                                                                                                                                             cacaGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGATGACGGCCAACAAATTg          <  1:1479751/150‑1 (MQ=255)
                                                                                                                                                aGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGATGACGGCCAACAAATTGCAt       <  2:1720399/150‑1 (MQ=255)
                                                                                                                                                    tcttctTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGATGACGGCCAACAAATTGCATAAtt   >  2:513463/1‑150 (MQ=255)
                                                                                                                                                    tcttctTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGATGACGGCCAACAAATTGCATAAtt   <  2:1196563/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCGTCTACCACGATAACCGGTTTGTCTTTGTGCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGATGACGGCCAACAAATTGCATAATTT  >  CP009273/3778070‑3778368

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: