Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,778,219 C→A 12.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,2190CA12.6% 853.7 / 66.8 246D53Y (GAC→TAC) yibNputative rhodanese‑related sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base C (126/89);  new base A (16/15);  total (142/104)
Fisher's exact test for biased strand distribution p-value = 5.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 269 total aligned reads displayed.

ACCGTCTACCACGATAACCGGTTTGTCTTTGTGCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGATGACGGCCAACAAATTGCATAATTT  >  CP009273/3778070‑3778368
                                                                                                                                                     |                                                                                                                                                     
 ccGTCTACCACGATAACCGGTTTGTCTTTGTGCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGtct                                                                                                                                                      >  1:1117084/1‑150 (MQ=255)
 ccGTCTACCACGATAACCGGTTTGTCTTTGTGCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGtct                                                                                                                                                      >  2:692676/1‑150 (MQ=255)
 ccGTCTACCACGATAACCGGTTTGTCTTTGTGCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGtct                                                                                                                                                      <  1:1672443/150‑1 (MQ=255)
  cGTCTACCACGATAACCGGTTTGTCTTTGTGCTTCTCAAGCTCACCAACATTGTCGGCATTGATTTCGCTCGGCAACAGAATAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGtctt                                                                                                                                                     <  2:3565128/150‑1 (MQ=255)
        ccACGATAACCGGTTTGTCTTTGTGCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTttgt                                                                                                                                               >  2:182542/1‑150 (MQ=255)
         cACGATAACCGGTTTGTCTTTGTGCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTttgtt                                                                                                                                              >  1:2361056/1‑150 (MQ=255)
            gATAACCGGTTTGTCTTTGTGCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGAt                                                                                                                                           <  1:293853/150‑1 (MQ=255)
            gATAACCGGTTTGTCTTTGTGCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGAt                                                                                                                                           >  1:1951247/1‑150 (MQ=255)
               aaCCGGTTTGTCTTTGTGCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTATTCTTTGTTGATCAg                                                                                                                                        >  1:2832914/1‑150 (MQ=255)
                    gTTTGTCTTTGTGCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGAcgcg                                                                                                                                   <  2:2616777/150‑1 (MQ=255)
                    gTTTGTCTTTGTGCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTATTCTTTGTTGATCAGAcgcg                                                                                                                                   <  1:1215278/150‑1 (MQ=255)
                      ttGTCTTTGTGCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTa                                                                                                                                 <  1:1869454/150‑1 (MQ=255)
                      ttGTCTTTGTGCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTATTCTTTGTTGATCAGACGCGTa                                                                                                                                 >  1:3507669/1‑150 (MQ=255)
                       tGTCTTTGTGCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAg                                                                                                                                <  2:2948940/150‑1 (MQ=255)
                       tGTCTTTGTGCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAg                                                                                                                                <  1:2505058/150‑1 (MQ=255)
                        gTCTTTGTGCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGc                                                                                                                               >  2:4122170/1‑150 (MQ=255)
                            ttGTGCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTca                                                                                                                           <  2:1154558/150‑1 (MQ=255)
                             tgtgCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTcac                                                                                                                          >  2:2120151/1‑150 (MQ=255)
                             tgtgCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTcac                                                                                                                          >  1:943478/1‑150 (MQ=255)
                              gtgCTTCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTcacc                                                                                                                         <  2:1255782/150‑1 (MQ=255)
                                  ttCTCAAGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGa                                                                                                                     >  1:1218976/1‑150 (MQ=255)
                                       aaGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAAt                                                                                                                >  1:449236/1‑150 (MQ=255)
                                       aaGCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAAt                                                                                                                <  2:54422/150‑1 (MQ=255)
                                         gCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCa                                                                                                              <  2:1696732/150‑1 (MQ=255)
                                         gCTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCa                                                                                                              >  1:4093430/1‑150 (MQ=255)
                                          cTCACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCAc                                                                                                             >  2:2186831/1‑150 (MQ=255)
                                            cACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCt                                                                                                           <  1:617466/150‑1 (MQ=255)
                                            cACCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCt                                                                                                           >  2:2921687/1‑150 (MQ=255)
                                             aCCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCtt                                                                                                          >  2:1890207/1‑150 (MQ=255)
                                             aCCAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCtt                                                                                                          >  1:4145699/1‑150 (MQ=255)
                                              ccAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTc                                                                                                         <  2:2412059/150‑1 (MQ=255)
                                              ccAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTc                                                                                                         <  2:103113/150‑1 (MQ=255)
                                               cAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCa                                                                                                        <  1:4361624/150‑1 (MQ=255)
                                               cAACATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCa                                                                                                        >  2:77662/1‑150 (MQ=255)
                                                  cATTGTTGGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACtt                                                                                                     <  1:2064739/150‑1 (MQ=255)
                                                     tgttgGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCg                                                                                                  >  1:1091539/1‑150 (MQ=255)
                                                     tgttgGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCg                                                                                                  <  2:1781877/150‑1 (MQ=255)
                                                       ttgGCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAg                                                                                                >  1:1433434/1‑150 (MQ=255)
                                                         ggCTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGt                                                                                              >  2:474180/1‑150 (MQ=255)
                                                           cTTTGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCa                                                                                            >  1:391554/1‑150 (MQ=255)
                                                            tttGATTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAg                                                                                           >  1:551520/1‑150 (MQ=255)
                                                                aTTTCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTATTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGctc                                                                                       <  1:3978727/150‑1 (MQ=255)
                                                                  ttCGCTCGGCAACAGATTAATAGAACCTTCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCtt                                                                                     <  1:2876682/150‑1 (MQ=255)
                                                                  ttCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCtt                                                                                     <  1:1333352/150‑1 (MQ=255)
                                                                  ttCGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCtt                                                                                     <  1:4633868/150‑1 (MQ=255)
                                                                    cGCTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTaa                                                                                   >  1:1716925/1‑150 (MQ=255)
                                                                      cTCGGCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAAc                                                                                 >  2:4224067/1‑150 (MQ=255)
                                                                         ggCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTa                                                                              <  2:1227239/150‑1 (MQ=255)
                                                                         ggCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTa                                                                              <  2:3003212/150‑1 (MQ=255)
                                                                          gCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAg                                                                             <  2:3229381/150‑1 (MQ=255)
                                                                          gCAACAGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAg                                                                             <  1:1154563/150‑1 (MQ=255)
                                                                               aGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTcaca                                                                        >  2:3898112/1‑150 (MQ=255)
                                                                               aGATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTATTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTcaca                                                                        >  2:230053/1‑150 (MQ=255)
                                                                                gATTAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACaa                                                                       >  1:3398644/1‑150 (MQ=255)
                                                                                   tAATAGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGaa                                                                    <  1:3223216/150‑1 (MQ=255)
                                                                                       aGAACCTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGc                                                                >  1:2011248/1‑150 (MQ=255)
                                                                                           ccTGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGccacc                                                            <  1:3454405/150‑1 (MQ=255)
                                                                                             tGCGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAg                                                          <  2:3204478/150‑1 (MQ=255)
                                                                                               cGATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTa                                                        <  2:1957967/150‑1 (MQ=255)
                                                                                                gATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTaa                                                       <  1:530531/150‑1 (MQ=255)
                                                                                                gATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTaa                                                       <  2:4547787/150‑1 (MQ=255)
                                                                                                gATATGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTaa                                                       >  2:1017234/1‑150 (MQ=255)
                                                                                                 atatGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAAc                                                      <  1:1800988/150‑1 (MQ=255)
                                                                                                  tatGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAAcg                                                     >  1:875513/1‑150 (MQ=255)
                                                                                                  tatGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAAcg                                                     <  1:1430628/150‑1 (MQ=255)
                                                                                                  tatGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAAcg                                                     >  1:2439340/1‑150 (MQ=255)
                                                                                                  tatGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAAcg                                                     >  2:3619622/1‑150 (MQ=255)
                                                                                                    tGGCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAAcgcg                                                   <  2:2849796/150‑1 (MQ=255)
                                                                                                     ggCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGa                                                  >  1:4223700/1‑150 (MQ=255)
                                                                                                      gCCTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGAt                                                 >  2:1657215/1‑150 (MQ=255)
                                                                                                       ccTTTACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATc                                                <  1:1744195/150‑1 (MQ=255)
                                                                                                         tttACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCa                                              <  2:2786954/150‑1 (MQ=255)
                                                                                                          ttACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAg                                             >  1:2646073/1‑150 (MQ=255)
                                                                                                           tACGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTATTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGATGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAgg                                            <  2:626953/150‑1 (MQ=255)
                                                                                                            aCGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGc                                           >  2:1588046/1‑150 (MQ=255)
                                                                                                             cGGAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTCACGCGATCCAGGCg                                          >  1:4560674/1‑150 (MQ=255)
                                                                                                              ggAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGa                                         >  1:3715247/1‑150 (MQ=255)
                                                                                                              ggAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTATTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGa                                         >  2:2340085/1‑150 (MQ=255)
                                                                                                               gAAGTCATCACGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGAt                                        <  2:4184384/150‑1 (MQ=255)
                                                                                                                  gTCATCACGCTGACGTAAATCCACAACCACAGCGTATTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACt                                     <  2:2506242/150‑1 (MQ=255)
                                                                                                                      tcaCGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGt                                 <  1:3187367/150‑1 (MQ=255)
                                                                                                                       caCGCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTa                                <  2:403316/150‑1 (MQ=255)
                                                                                                                        aCGCTGACGTAAATCCACAACCACAGCGTATTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTAt                               >  2:1222809/1‑150 (MQ=255)
                                                                                                                          gCTGACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATgg                             <  1:2136561/150‑1 (MQ=255)
                                                                                                                             gACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGAt                          <  2:4291733/150‑1 (MQ=255)
                                                                                                                             gACGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGAt                          >  2:612791/1‑150 (MQ=255)
                                                                                                                              aCGTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGATg                         <  2:3765971/150‑1 (MQ=255)
                                                                                                                                gTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGATGAc                       <  2:2977819/150‑1 (MQ=255)
                                                                                                                                gTAAATCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGATGAc                       <  1:2529013/150‑1 (MQ=255)
                                                                                                                                  aaaTCCACAACCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGATGACgg                     <  1:700778/150‑1 (MQ=255)
                                                                                                                                          aaCCACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGATGACGGCCAACaaa             >  1:2720202/1‑150 (MQ=255)
                                                                                                                                            ccACAGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGATGACGGCCAACAAAtt           <  1:3336841/150‑1 (MQ=255)
                                                                                                                                                aGCGTCTTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGATGACGGCCAACAAATTGCAt       >  1:844122/1‑150 (MQ=255)
                                                                                                                                                  cGTATTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGATGACGGCCAACAAATTGCATaa     <  1:1553850/150‑1 (MQ=255)
                                                                                                                                                  cGTATTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGATGACGGCCAACAAATTGCATaa     <  1:3611833/150‑1 (MQ=255)
                                                                                                                                                  cGTATTCTTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGATGACGGCCAACAAATTGCATaa     <  1:3611811/150‑1 (MQ=255)
                                                                                                                                                    tcttctTTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGATGACGGCCAACAAATTGCATAAtt   <  2:3336473/150‑1 (MQ=255)
                                                                                                                                                     cttcttTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGATGACGGCCAACAAATTGCATAAttt  <  2:3986572/150‑1 (MQ=255)
                                                                                                                                                     cttcttTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGATGACGGCCAACAAATTGCATAAttt  <  1:182542/150‑1 (MQ=255)
                                                                                                                                                     cttcttTGTTGATCAGACGCGTAGCTTCACCACGAGTAATCACCTTCACTTTCGAGGTCAGGCTCTTAAACGTAGTCACAAGAACCGCCACCAGTAACGCGATCCAGGCGATACTCAGTATGGGATGACGGCCAACAAATTGCATAAttt  >  1:2690810/1‑150 (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: