Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,209,074 T→A 8.3% *244K (TAA→AAA)  mlrA → transcriptional activator of csgD and csgBA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,209,0740TA 8.3% 1020.2 / 60.6 313*244K (TAA→AAA) mlrAtranscriptional activator of csgD and csgBA
Reads supporting (aligned to +/- strand):  ref base T (138/149);  new base A (12/14);  total (150/163)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 320 total aligned reads displayed.

CCTCGCCCACTCGCTCAATCAACTACGCCCTGAACTATTCGAAGGCCGTACATTGCTGGTGTGGTGCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCT  >  CP009273/2208926‑2209223
                                                                                                                                                    |                                                                                                                                                     
ccTCGCCCACTCGCTCAATCAACTACGCCCTGAACTATTCGAAGGCCGTACATTGCTGGTGTGGTGCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTa                                                                                                                                                      <  1:492942/150‑1 (MQ=255)
ccTCGCCCACTCGCTCAATCAACTACGCCCTGAACTATTCGAAGGCCGTACATTGCTGGTGTGGTGCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTa                                                                                                                                                      >  1:2299772/1‑150 (MQ=255)
     cccACTCGCTCAATCAACTACGCCCTGAACTATTCGAAGGCCGTACATTGCTGGTGTGGTGCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGAt                                                                                                                                                 <  2:2884517/150‑1 (MQ=255)
       cACTCGCTCAATCAACTACGCCCTGAACTATTCGAAGGCCGTACATTGCTGGTGTGGTGCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTc                                                                                                                                               <  2:1073937/150‑1 (MQ=255)
       cACTCGCTCAATCAACTACGCCCTGAACTATTCGAAGGCCGTACATTGCTGGTGTGGTGCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTc                                                                                                                                               <  2:4262983/150‑1 (MQ=255)
        aCTCGCTCAATCAACTACGCCCTGAACTATTCGAAGGCCGTACATTGCTGGTGTGGTGCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCg                                                                                                                                              <  2:3775234/150‑1 (MQ=255)
         cTCGCTCAATCAACTACGCCCTGAACTATTCGAAGGCCGTACATTGCTGGTGTGGTGCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGt                                                                                                                                             >  1:270633/1‑150 (MQ=255)
          tCGCTCAATCAACTACGCCCTGAACTATTCGAAGGCCGTACATTGCTGGTGTGGTGCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGtt                                                                                                                                            >  2:3258073/1‑150 (MQ=255)
          tCGCTCAATCAACTACGCCCTGAACTATTCGAAGGCCGTACATTGCTGGTGTGGTGCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGtt                                                                                                                                            <  2:1412727/150‑1 (MQ=255)
           cGCTAAATAAACTACGCCCTGAACTATTCGAAGTCCGTAGATGGCTGTTGTGGTGCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAATGCCATGTTATTTTCCCACTCGGCATTTAATGATTCGTTa                                                                                                                                           <  2:3244899/150‑1 (MQ=255)
            gCTCAATCAACTACGCCCTGAACTATTCGAAGGCCGTACATTGCTGGTGTGGTGCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTaa                                                                                                                                          >  2:3587827/1‑150 (MQ=255)
              tcaatcaaCTACGCCCTGAACTATTCGAAGGCCGTACATTGCTGGTGTGGTGCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACa                                                                                                                                        <  1:4162150/150‑1 (MQ=255)
              tcaatcaaCTACGCCCTGAACTATTCGAAGGCCGTACATTGCTGGTGTGGTGCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACa                                                                                                                                        <  2:2128234/150‑1 (MQ=255)
                   caaCTACGCCCTGAACTATTCGAAGGCCGTACATTGCTGGTGTGGTGCGGTGAAAATCGAACCTCCGTCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTAAATGATTCGTTAACAAATgc                                                                                                                                   >  1:456373/1‑150 (MQ=255)
                   caaCTACGCCCTGAACTATTCGAAGGCCGTACATTGCTGGTGTGGTGCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATgc                                                                                                                                   >  1:4387357/1‑150 (MQ=255)
                        aCGCCCTGAACTATTCGAAGGCCGTACATTGCTGGTGTGGTGCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCttt                                                                                                                              <  2:133160/150‑1 (MQ=255)
                        aCGCCCTGAACTATTCGAAGGCCGTACATTGCTGGTGTGGTGCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCttt                                                                                                                              >  1:199992/1‑150 (MQ=255)
                             cTGAACTATTCGAAGGCCGTACATTGCTGGTGTGGTGCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGt                                                                                                                         >  1:680064/1‑150 (MQ=255)
                             cTGAACTATTCGAAGGCCGTACATTGCTGGTGTGGTGCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGt                                                                                                                         >  1:1815318/1‑150 (MQ=255)
                               gAACTATTCGAAGGCCGTACATTGCTGGTGTGGTGCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTAc                                                                                                                       <  1:899468/150‑1 (MQ=255)
                                aaCTATTCGAAGGCCGTACATTGCTGGTGTGGTGCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACa                                                                                                                      <  1:2248909/150‑1 (MQ=255)
                                   tATTCGAAGGCCGTACATTGCTGGTGTGGTGCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATc                                                                                                                   <  1:2566351/150‑1 (MQ=255)
                                      tCGAAGGCCGTACATTGCTGGTGTGGTGCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCtt                                                                                                                <  1:783652/150‑1 (MQ=255)
                                       cGAAGGCCGTACATTGCTGGTGTGGTGCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCttt                                                                                                               >  2:146009/1‑150 (MQ=255)
                                           ggCCGTACATTGCTGGTGTGGTGCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGtt                                                                                                           <  2:2985801/150‑1 (MQ=255)
                                            gCCGTACATTGCTGGTGTGGTGCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTAAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTa                                                                                                          <  1:3925459/150‑1 (MQ=255)
                                            gCCGTACATTGCTGGTGTGGTGCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTAAATGATTCGTCAACAAATGCGCTTTACTGTACAATCCTTTCGTTa                                                                                                          <  1:980800/150‑1 (MQ=255)
                                              cGTACATTGCTGGTGTGGTGCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAAc                                                                                                        >  2:489337/1‑150 (MQ=255)
                                               gTACATTGCTGGTGTGGTGCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACa                                                                                                       >  1:110596/1‑150 (MQ=255)
                                                   aTTGCTGGTGTGGTGCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAg                                                                                                   >  1:4014210/1‑150 (MQ=255)
                                                   aTTGCTGGTGTGGTGCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTAAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAg                                                                                                   >  1:55602/1‑150 (MQ=255)
                                                       cTGGTGTGGTGCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGt                                                                                               <  1:2978613/150‑1 (MQ=255)
                                                        tggtgtggtgCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTAAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTg                                                                                              <  1:787227/150‑1 (MQ=255)
                                                          gtgtggtgCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCa                                                                                            <  2:3897333/150‑1 (MQ=255)
                                                            gtggtgCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCAtt                                                                                          <  2:3628671/150‑1 (MQ=255)
                                                            gtggtgCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCAtt                                                                                          <  1:25768/150‑1 (MQ=255)
                                                             tggtgCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGGTAACATAAGGAGTGCATTa                                                                                         >  1:1672121/1‑150 (MQ=255)
                                                               gtgCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATg                                                                                       >  2:2018377/1‑150 (MQ=255)
                                                               gtgCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATg                                                                                       <  1:960612/150‑1 (MQ=255)
                                                                tgCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGAGATTTTCCCACTCGGCATTTAATGATTAGTTAACAAATGCGCTTTACTGTACAATCATTTCGTTAACATAAGGAGTGCATTATgc                                                                                      <  1:3542471/150‑1 (MQ=255)
                                                                  cGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATgcgc                                                                                    <  2:1815318/150‑1 (MQ=255)
                                                                   ggTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAATGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCa                                                                                   >  2:2719907/1‑150 (MQ=255)
                                                                   ggTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCa                                                                                   >  2:4078689/1‑150 (MQ=255)
                                                                    gTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCAt                                                                                  <  2:3991826/150‑1 (MQ=255)
                                                                    gTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCAt                                                                                  <  2:1546083/150‑1 (MQ=255)
                                                                    gTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCAt                                                                                  >  1:1293068/1‑150 (MQ=255)
                                                                    gTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCAt                                                                                  <  2:1246747/150‑1 (MQ=255)
                                                                        aaaTCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCt                                                                              >  2:2811382/1‑150 (MQ=255)
                                                                        aaaTCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGAGCTTTACTGTACAAACCTTTCGATAAAATAAGGAGTGCATTATGCGCATAGCt                                                                              >  2:3044290/1‑150 (MQ=255)
                                                                         aaTCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTa                                                                             >  1:3316449/1‑150 (MQ=255)
                                                                         aaTCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTa                                                                             <  1:1007681/150‑1 (MQ=255)
                                                                            cGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTAAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCCTAGCTaaaa                                                                          >  1:3202307/1‑150 (MQ=255)
                                                                             gAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAAt                                                                         >  1:3440091/1‑150 (MQ=255)
                                                                              aaCCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAAtt                                                                        <  1:4499900/150‑1 (MQ=255)
                                                                              aaCCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGAATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAAtt                                                                        <  1:4443085/150‑1 (MQ=255)
                                                                               aCCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTg                                                                       <  2:4145507/150‑1 (MQ=255)
                                                                                 cTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTggg                                                                     <  2:2054288/150‑1 (MQ=255)
                                                                                  tCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTgggg                                                                    <  2:2596079/150‑1 (MQ=255)
                                                                                     gCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCAATCGGCATTAAATGATTCGTTAACAAATGCGCTTTACTGTACAAACCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCa                                                                 >  2:3937290/1‑150 (MQ=255)
                                                                                       ccAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCAta                                                               >  2:3654217/1‑149 (MQ=255)
                                                                                        cAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTAAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCAtag                                                              >  2:2257371/1‑148 (MQ=255)
                                                                                           cagcagCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGccc                                                           >  1:4282515/1‑150 (MQ=255)
                                                                                              cagcaACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCcctgt                                                        >  2:3426012/1‑150 (MQ=255)
                                                                                               agcaACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCcctgtt                                                       >  1:116572/1‑150 (MQ=255)
                                                                                                  aaCTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCcctgttcct                                                    >  2:4026661/1‑150 (MQ=255)
                                                                                                  aaCTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCcctgttcct                                                    <  1:774391/150‑1 (MQ=255)
                                                                                                     tCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCcctgttcctgtt                                                 >  2:681717/1‑150 (MQ=255)
                                                                                                      cccccggtGGCAAGAACACGGCCAGGACATCTACCCACTCGGCATTTACTGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGttt                                                <  1:4295291/143‑1 (MQ=255)
                                                                                                      cACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGttt                                                <  1:3462148/150‑1 (MQ=255)
                                                                                                       aCCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTAAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTa                                               >  2:1608061/1‑150 (MQ=255)
                                                                                                        ccAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTAt                                              <  2:1744661/150‑1 (MQ=255)
                                                                                                          aGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATgg                                            <  1:59584/150‑1 (MQ=255)
                                                                                                              ggCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGtt                                        >  1:3905125/1‑150 (MQ=255)
                                                                                                              ggCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGtt                                        <  1:2508431/150‑1 (MQ=255)
                                                                                                                 aaGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAgg                                     >  1:1131964/1‑150 (MQ=255)
                                                                                                                    aaCAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCgg                                  >  2:2513043/1‑150 (MQ=255)
                                                                                                                     aCAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGTCCTGTTCCTGTTTATGGCGTTAGGCGGa                                 <  1:3258073/150‑1 (MQ=255)
                                                                                                                     aCAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGa                                 <  1:3399265/150‑1 (MQ=255)
                                                                                                                             aTGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATtggtgg                         <  1:760419/150‑1 (MQ=255)
                                                                                                                              tGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGc                        <  2:407831/150‑1 (MQ=255)
                                                                                                                              tGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGc                        <  2:1294701/150‑1 (MQ=255)
                                                                                                                              tGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGc                        <  2:906226/150‑1 (MQ=255)
                                                                                                                               gATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCg                       <  2:1959493/150‑1 (MQ=255)
                                                                                                                                atatTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGt                      >  2:1148975/1‑150 (MQ=255)
                                                                                                                                 tatTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTc                     >  1:1396101/1‑150 (MQ=255)
                                                                                                                                 tatTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGACATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTc                     >  1:3664403/1‑150 (MQ=255)
                                                                                                                                  atTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCa                    <  2:31207/150‑1 (MQ=255)
                                                                                                                                  atTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCa                    <  2:31026/150‑1 (MQ=255)
                                                                                                                                     ttCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGACTTGGTGGCGTCATGc                 >  2:1535116/1‑150 (MQ=255)
                                                                                                                                      tCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCt                <  2:431830/150‑1 (MQ=255)
                                                                                                                                       cccACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTc               >  2:1057773/1‑150 (MQ=255)
                                                                                                                                        ccACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCg              <  2:3732297/150‑1 (MQ=255)
                                                                                                                                        ccACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCg              <  2:555287/150‑1 (MQ=255)
                                                                                                                                           cTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAg           <  2:3275204/150‑1 (MQ=255)
                                                                                                                                           cTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCAGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAg           >  1:3720910/1‑150 (MQ=255)
                                                                                                                                           cTCGGCATTAAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAg           >  2:817505/1‑150 (MQ=255)
                                                                                                                                             cGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGt         >  1:3590459/1‑150 (MQ=255)
                                                                                                                                                 aTTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTtata     >  1:2163508/1‑150 (MQ=255)
                                                                                                                                                 aTTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGAATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTtata     >  1:4404847/1‑150 (MQ=255)
                                                                                                                                                    tAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATAGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCt  <  1:106766/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCTCGCCCACTCGCTCAATCAACTACGCCCTGAACTATTCGAAGGCCGTACATTGCTGGTGTGGTGCGGTGAAAATCGAACCTCCGCCCAACAGCAGCAACTCACCAGTTGGCAAGAACAAGGCCATGATATTTTCCCACTCGGCATTTAATGATTCGTTAACAAATGCGCTTTACTGTACAATCCTTTCGTTAACATAAGGAGTGCATTATGCGCATAGCTAAAATTGGGGTCATCGCCCTGTTCCTGTTTATGGCGTTAGGCGGAATTGGTGGCGTCATGCTCGCAGGTTATACCT  >  CP009273/2208926‑2209223

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: