Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,920,013 G→T 8.4% D117E (GAC→GAA asnC ← transcriptional activator of asnA; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,920,0130GT 8.4% 1049.5 / 45.6 285D117E (GAC→GAAasnCtranscriptional activator of asnA; autorepressor
Reads supporting (aligned to +/- strand):  ref base G (116/145);  new base T (6/18);  total (122/163)
Fisher's exact test for biased strand distribution p-value = 8.42e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 285 total aligned reads displayed.

GGGATCTACCTGTGGAAAAGTATGGGATTAAAAAAGCCGATCAGGGCTTGATGGTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCT  >  CP009273/3919865‑3920160
                                                                                                                                                    |                                                                                                                                                   
gggATCTACCTGTGGAAAAGTATGGGATTAAAAAAGCCGATCAGGGCTTGATGGTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGt                                                                                                                                                    <  2:3576128/150‑1 (MQ=255)
 ggATCTACCTGTGGAAAAGTATGGGATTAAAAAAGCCGATCAGGGCTTGATGGTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGtc                                                                                                                                                   <  1:2665495/150‑1 (MQ=255)
      tACCTGTGGAAAAGTATGGGATTAAAAAAGCCGATCAGGGCTTGATGGTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGtcgatcg                                                                                                                                              <  1:960218/150‑1 (MQ=255)
      tACCTGTGGAAAAGTATGGGATTAAAAAAGCCGATCAGGGCTTGATGGTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGtcgatcg                                                                                                                                              >  1:1352687/1‑150 (MQ=255)
       cccggtgGAAAAGTATGGGATTAAAAAAGCCGATCAGGGCTTGATGGTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCTtcgatcgc                                                                                                                                             <  1:1057209/146‑2 (MQ=255)
         cTGTGGAAAAGTATGGGATTAAAAAAGCCGATCAGGGCTTGATGGTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCTTCGATCGAAc                                                                                                                                           <  2:2982005/150‑1 (MQ=255)
         cTGTGGAAAAGTATGGGATTAAAAAAGCCGATCAGGGCTTGATGGTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAAc                                                                                                                                           <  2:1876389/150‑1 (MQ=255)
          tgtgGAAAAGTATGGGATTAAAAAAGCCGATCAGGGCTTGATGGTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACg                                                                                                                                          <  1:2494862/150‑1 (MQ=255)
            tgGAAAAGTATGGGATTAAAAAAGCCGATCAGGGCTTGATGGTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGc                                                                                                                                        <  2:3536608/150‑1 (MQ=255)
               aaaaGTATGGGATTAAAAAAGCCGATCAGGGCTTGATGGTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGcaca                                                                                                                                     >  1:3603026/1‑150 (MQ=255)
                  aGTATGGGATTAAAAAAGCCGATCAGGGCTTGATGGTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCa                                                                                                                                  <  1:4183720/150‑1 (MQ=255)
                      tGGGATTAAAAAAGCCGATCAGGGCTTGATGGTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACttt                                                                                                                              <  2:4449669/150‑1 (MQ=255)
                       gggATTAAAAAAGCCGATCAGGGCTTGATGGTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACtttt                                                                                                                             <  1:867177/150‑1 (MQ=255)
                        ggATTATAAAAGCCGATCAGGGCTTGATGGTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTta                                                                                                                            <  2:3171809/150‑1 (MQ=255)
                        ggATTAAAAAAGCCGATCAGGGCTTGATGGTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCAATTTta                                                                                                                            >  1:3277347/1‑150 (MQ=255)
                          aTTAAAAAAGCCGATCAGGGCTTGATGGTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTtata                                                                                                                          <  2:401642/150‑1 (MQ=255)
                          aTTAAAAAAGCCGATCAGGGCTTGATGGTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCGTGTTGATCAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTGTttta                                                                                                                          >  2:2599179/1‑150 (MQ=255)
                           ttAAAAAAGCCGATCAGGGATTGATGGTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATaa                                                                                                                         <  2:2122371/150‑1 (MQ=255)
                            tAAAAAAGCCGATCAGGGCTTGATGGTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATaaa                                                                                                                        >  1:2276327/1‑150 (MQ=255)
                                 aaGCCGATCAGGGCTTGATGGTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGc                                                                                                                   >  2:420006/1‑150 (MQ=255)
                                  aGCCGATCAGGGCTTGATGGTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCt                                                                                                                  <  1:3349115/150‑1 (MQ=255)
                                   gCCGATCAGGGCTTGATGGTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACAGCACTTTTATAAAGATGCTg                                                                                                                 >  2:550269/1‑150 (MQ=255)
                                   gCCGATCAGGGCTTGATGGTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGACCATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTg                                                                                                                 <  1:3693307/150‑1 (MQ=255)
                                       aTCAGGGCTTGATGGTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGtcgg                                                                                                             >  2:130929/1‑147 (MQ=255)
                                       aTCAGGGCTTGATGGTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGt                                                                                                             >  1:4106946/1‑150 (MQ=255)
                                           gggCTTGATGGTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGcc                                                                                                         >  2:4138614/1‑150 (MQ=255)
                                               ttGATGGTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGtt                                                                                                     <  1:3916910/150‑1 (MQ=255)
                                                tGATGGTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTtg                                                                                                    <  2:3979601/150‑1 (MQ=255)
                                                tGATGGTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTtg                                                                                                    >  1:4022600/1‑150 (MQ=255)
                                                 gATGGTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTtgt                                                                                                   <  2:639353/150‑1 (MQ=255)
                                                    ggTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTa                                                                                                >  1:891330/1‑150 (MQ=255)
                                                    ggTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTa                                                                                                <  2:7391/150‑1 (MQ=255)
                                                    ggTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTGGTGTa                                                                                                >  2:1347890/1‑150 (MQ=255)
                                                         gCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAgg                                                                                           >  1:2158010/1‑150 (MQ=255)
                                                         gCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAgg                                                                                           >  1:3642441/1‑150 (MQ=255)
                                                           aTGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGGTTGGATCTTGTTGGTAAGTACATGCTGGAGAGCGGCGATCGACCGGCACATCACTTGTATAAAGATGCTGTAGTGGGCGGGGGTGTATTAGGcg                                                                                         >  2:1079421/1‑149 (MQ=255)
                                                             gATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTc                                                                                       <  1:2360824/150‑1 (MQ=255)
                                                                 gggTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCGTCAGTg                                                                                   >  2:3838275/1‑150 (MQ=255)
                                                                    ttttGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACt                                                                                <  2:1161761/150‑1 (MQ=255)
                                                                      ttGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTtc                                                                              <  2:3100305/150‑1 (MQ=255)
                                                                      ttGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGGTGTGTAATAGGCGTCAGTGACTtc                                                                              >  2:815187/1‑150 (MQ=255)
                                                                         cAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTtcatc                                                                           <  2:3350724/150‑1 (MQ=255)
                                                                         cAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTtcatc                                                                           >  1:3654907/1‑150 (MQ=255)
                                                                          aGGACGATAAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTtcatca                                                                          <  2:4392977/150‑1 (MQ=255)
                                                                            gACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAg                                                                        >  1:481978/1‑150 (MQ=255)
                                                                             aCGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAgg                                                                       <  1:1987799/150‑1 (MQ=255)
                                                                               gATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCTTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCt                                                                     <  1:2101326/150‑1 (MQ=255)
                                                                                 tctaTGTCTCGGTGGACTGAATTTCATCGATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCttt                                                                   <  1:3538351/147‑1 (MQ=255)
                                                                                   aaTGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTcc                                                                 <  2:2430298/150‑1 (MQ=255)
                                                                                    aTGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCa                                                                <  2:162705/150‑1 (MQ=255)
                                                                                     tGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAg                                                               <  1:4398615/150‑1 (MQ=255)
                                                                                        ctcGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCtt                                                            <  1:38475/150‑1 (MQ=255)
                                                                                           ggTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGc                                                         <  2:2464064/150‑1 (MQ=255)
                                                                                            gTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGcc                                                        <  2:1488607/150‑1 (MQ=255)
                                                                                             tGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCa                                                       >  1:1143250/1‑150 (MQ=255)
                                                                                              ggACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAg                                                      <  2:1765726/150‑1 (MQ=255)
                                                                                               gACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAgc                                                     <  1:4297939/150‑1 (MQ=255)
                                                                                                  tGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGcgcg                                                  >  1:213442/1‑150 (MQ=255)
                                                                                                  tGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGcgcg                                                  <  1:275956/150‑1 (MQ=255)
                                                                                                  tGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGcgcg                                                  <  2:1345014/150‑1 (MQ=255)
                                                                                                   gAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCgg                                                 >  1:1899671/1‑150 (MQ=255)
                                                                                                    aaTTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGa                                                >  2:3729286/1‑150 (MQ=255)
                                                                                                       ttCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTTTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACGTCATCAAGGCTTTCCCGCATTGCCAGAGCGGAAgg                                             >  2:1665347/1‑150 (MQ=255)
                                                                                                       ttCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAgg                                             <  1:3930199/150‑1 (MQ=255)
                                                                                                        tcttcaATTGTTAGGATATTGTTGATAAGTACATGCTGGAGAGCTTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAggg                                            <  2:2916225/147‑1 (MQ=255)
                                                                                                         catcaATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGt                                           >  2:280106/1‑150 (MQ=255)
                                                                                                          atcaATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCTTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTa                                          >  2:1244127/1‑150 (MQ=255)
                                                                                                           tcaATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGATGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAg                                         >  2:2221504/1‑150 (MQ=255)
                                                                                                            caAGTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCTTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGt                                        <  1:78729/150‑1 (MQ=255)
                                                                                                             aaTTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTc                                       >  1:2210680/1‑150 (MQ=255)
                                                                                                                tGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCttt                                    <  2:1352690/150‑1 (MQ=255)
                                                                                                                tGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCttt                                    <  2:367772/150‑1 (MQ=255)
                                                                                                                tGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCttt                                    <  1:3849429/150‑1 (MQ=255)
                                                                                                                  tttGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTgg                                  <  2:1130300/150‑1 (MQ=255)
                                                                                                                  tttGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTgg                                  >  2:1698873/1‑150 (MQ=255)
                                                                                                                   ttGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGgc                                 <  1:1206651/150‑1 (MQ=255)
                                                                                                                    tGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGgcg                                >  2:2618049/1‑150 (MQ=255)
                                                                                                                    tGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACAACACTTTTATAAAGATGCTGTAGTGGCCGGTTTTGTAATAGGCTTCAGTGACTTAATCAAGGCTTTACAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGgcg                                >  2:2017436/1‑150 (MQ=255)
                                                                                                                      gATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGct                              >  1:1365903/1‑150 (MQ=255)
                                                                                                                      gATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGct                              >  1:3759355/1‑150 (MQ=255)
                                                                                                                       aTCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGctc                             >  2:3940612/1‑150 (MQ=255)
                                                                                                                        tCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGctct                            <  1:3325553/150‑1 (MQ=255)
                                                                                                                          ttgttgATAAGTTCATGCTGGAGCGCGTCGATCTCACGGCACATCTCTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCttt                          <  2:277969/150‑1 (MQ=255)
                                                                                                                           tgttgATAAGTACATGCTGGAGAGCTTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTa                         >  2:1948037/1‑150 (MQ=255)
                                                                                                                             ttgATAAGTACATGCTGGAGAGCTTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAat                       <  1:3254102/150‑1 (MQ=255)
                                                                                                                               gATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAatat                     <  1:2860748/150‑1 (MQ=255)
                                                                                                                                aTAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAAtata                    <  1:4403806/150‑1 (MQ=255)
                                                                                                                                   aGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATg                 <  1:1719380/150‑1 (MQ=255)
                                                                                                                                     tACATGCTGGAGAGCGTCGATCGAACTGAACATCACTTTTATAATGATGGTGTAGTGGCCGGTTTTGTTATAGGCTTCAGTGACTTCATAAATGCTTTCAAGCTTTGCCAGCGCGGATGGGTAGTCTTTGGCGCTCTttaatttaattca               <  2:4284703/150‑4 (MQ=255)
                                                                                                                                     tACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGcc               >  2:138611/1‑150 (MQ=255)
                                                                                                                                     tACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGcc               <  1:4004749/150‑1 (MQ=255)
                                                                                                                                      aCATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCg              <  2:330000/150‑1 (MQ=255)
                                                                                                                                      aCATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCg              <  2:672212/150‑1 (MQ=255)
                                                                                                                                       cATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGa             >  2:1900900/1‑150 (MQ=255)
                                                                                                                                        aTGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGAt            >  2:1546670/1‑150 (MQ=255)
                                                                                                                                           cTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATaaa         >  1:3470280/1‑150 (MQ=255)
                                                                                                                                             ggAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAagc       <  2:3802751/150‑1 (MQ=255)
                                                                                                                                              gagaGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTTGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAagca      <  1:1083198/150‑1 (MQ=255)
                                                                                                                                              gagaGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAagca      <  1:1429327/150‑1 (MQ=255)
                                                                                                                                                  gCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCt  <  1:3114241/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGGATCTACCTGTGGAAAAGTATGGGATTAAAAAAGCCGATCAGGGCTTGATGGTACGCATGATCGGGTTTTGCAGGACGATCAATGTCTCGGTGGACTGAATTTCATCAATTGTTTGGATCTTGTTGATAAGTACATGCTGGAGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGTAATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCT  >  CP009273/3919865‑3920160

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