Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,642,278 T→C 13.0% intergenic (+37/‑17) arsR → / → arsB arsenical resistance operon transcriptional repressor; autorepressor/arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,642,2780TC13.0% 687.7 / 80.0 241intergenic (+37/‑17)arsR/arsBarsenical resistance operon transcriptional repressor; autorepressor/arsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base T (112/97);  new base C (20/11);  total (133/108)
Fisher's exact test for biased strand distribution p-value = 3.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 251 total aligned reads displayed.

CGGCGAAAATTATTGATGAGGCCTGGCGATGTGAACAGGAAAAGGTTCAGGCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGTACTGGCGTTAGTTACGGGCGTGGTCCATCCGGGTG  >  CP009273/3642129‑3642427
                                                                                                                                                     |                                                                                                                                                     
 ggcgAAAATTATTGATGAGGCCTGGCGATGTGAACAGGAAAAGGTTCAGGCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGttt                                                                                                                                                      <  2:803447/150‑1 (MQ=255)
   cgAAAATTATTGATGAGGCCTGGCGATGTGAACAGGAAAAGGTTCAGGCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTa                                                                                                                                                    <  1:3464709/150‑1 (MQ=255)
     aaaaTTATTGATGAGGCCTGGCGATGTGAACAGGAAAAGGTTCAGGCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATc                                                                                                                                                  <  2:2076502/150‑1 (MQ=255)
         ttattGATGAGGCCTGGCGATGTGAACAGGAAAAGGTTCAGGCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCggg                                                                                                                                              <  1:1261158/150‑1 (MQ=255)
           attGATGAGGCCTGGCGATGTGAACAGGAAAAGGTTCAGGCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAg                                                                                                                                            >  2:2357324/1‑150 (MQ=255)
              gatgaGGCCTGGCGATGTGAACAGGAAAAGGTTCAGGCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCa                                                                                                                                         <  2:3254832/150‑1 (MQ=255)
                 gaGGCCTGGCGATGTGAACAGGAAAAGGTTCAGGCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTa                                                                                                                                      <  1:4035785/150‑1 (MQ=255)
                 gaGGCCTGGCGATGTGAACAGGAAAAGGTTCAGGCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTa                                                                                                                                      <  2:1552417/150‑1 (MQ=255)
                 gaGGCCTGGCGATGTGAACAGGAAAAGGTTCAGGCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTa                                                                                                                                      >  2:4112032/1‑150 (MQ=255)
                  aGGCCTGGCGATGTGAACAGGAAAAGGTTCAGGCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTAt                                                                                                                                     <  1:1613730/150‑1 (MQ=255)
                  aGGCCTGGCGATGTGAACAGGAAAAGGTTCAGGCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTAt                                                                                                                                     >  2:4553388/1‑150 (MQ=255)
                  aGGCCTGGCGATGTGAACAGGAAAAGGTTCAGGCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTCTTATCCGGGAGGCATTAt                                                                                                                                     >  2:2674475/1‑150 (MQ=255)
                   ggCCTGGCGATGTGAACAGGAAAAGGTTCAGGCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATg                                                                                                                                    >  2:2631959/1‑150 (MQ=255)
                          cGATGTGAACAGGAAAAGGTTCAGGCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTCTTATCCGGGAGGCATTATGTTACTgg                                                                                                                             >  1:1987412/1‑150 (MQ=255)
                            aTGTGAACAGGAAAAGGTTCAGGCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCa                                                                                                                           <  2:3164195/150‑1 (MQ=255)
                            aTGTGAACAGGAAAAGGTTCAGGCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCa                                                                                                                           <  2:4465553/150‑1 (MQ=255)
                             tgtgAACAGGAAAAGGTTCAGGCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAg                                                                                                                          <  2:733045/150‑1 (MQ=255)
                              gtgAACAGGAAAAGGTTCAGGCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAgg                                                                                                                         >  1:2580510/1‑150 (MQ=255)
                              gtgAACAGGAAAAGGTTCAGGCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAgg                                                                                                                         >  1:4699570/1‑150 (MQ=255)
                              gtgAACAGGAAAAGGTTCAGGCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAgg                                                                                                                         >  1:127419/1‑150 (MQ=255)
                                       aaaaGGTTCAGGCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCtt                                                                                                                >  1:3563506/1‑150 (MQ=255)
                                       aaaaGGTTCAGGCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCtt                                                                                                                <  1:2788954/150‑1 (MQ=255)
                                        aaaGGTTCAGGCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCttt                                                                                                               >  1:3341782/1‑150 (MQ=255)
                                             ttCAGGCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTCTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCt                                                                                                          <  1:4105512/150‑1 (MQ=255)
                                              tCAGGCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTg                                                                                                         <  2:2530446/150‑1 (MQ=255)
                                               cAGGCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTCTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCTTGa                                                                                                        <  1:93945/150‑1 (MQ=255)
                                                aGGCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGAc                                                                                                       <  2:188787/150‑1 (MQ=255)
                                                aGGCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGAc                                                                                                       >  2:3368235/1‑150 (MQ=255)
                                                aGGCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGAc                                                                                                       >  1:3547481/1‑150 (MQ=255)
                                                aGGCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGAc                                                                                                       >  2:2800610/1‑150 (MQ=255)
                                                  gCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCa                                                                                                     <  1:845621/150‑1 (MQ=255)
                                                  gCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCa                                                                                                     <  1:3772967/150‑1 (MQ=255)
                                                    gATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATc                                                                                                   >  2:1763887/1‑150 (MQ=255)
                                                        gTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTAt                                                                                               <  2:4087049/150‑1 (MQ=255)
                                                                ccTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATc                                                                                       <  1:548227/150‑1 (MQ=255)
                                                                ccTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATc                                                                                       <  1:2256615/150‑1 (MQ=255)
                                                                   ggCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTCTTGGTTATCTgg                                                                                    >  1:573848/1‑150 (MQ=255)
                                                                     cTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCa                                                                                  <  2:4314130/150‑1 (MQ=255)
                                                                         aCAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCg                                                                              <  2:4376014/150‑1 (MQ=255)
                                                                            aaaCTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACCTATGAATTTTCAGATGTGTCTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGaaa                                                                           <  1:196274/150‑1 (MQ=255)
                                                                              aCTGTTCCGGGGACAGTAAGAACATTTGGAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTCTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAgg                                                                         >  1:1350902/1‑150 (MQ=255)
                                                                              aCTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTCTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAgg                                                                         >  1:1351200/1‑150 (MQ=255)
                                                                                tGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTCTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGtt                                                                       <  1:2985005/150‑1 (MQ=255)
                                                                                ggTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGCATTGTCAGAGGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGtt                                                                       <  2:1110341/149‑1 (MQ=255)
                                                                                 gTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGttt                                                                      <  2:1060889/150‑1 (MQ=255)
                                                                                 gTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGttt                                                                      <  2:299601/150‑1 (MQ=255)
                                                                                  ttCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTa                                                                     <  2:2960921/150‑1 (MQ=255)
                                                                                    ccGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAgg                                                                   >  1:430563/1‑150 (MQ=255)
                                                                                    ccGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAgg                                                                   <  1:389955/150‑1 (MQ=255)
                                                                                     cGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGc                                                                  <  2:4202454/150‑1 (MQ=255)
                                                                                      ggggACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCa                                                                 >  2:1748901/1‑150 (MQ=255)
                                                                                      ggggACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCa                                                                 <  2:3579338/150‑1 (MQ=255)
                                                                                      ggggACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCa                                                                 >  2:3623773/1‑150 (MQ=255)
                                                                                      ggggACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCa                                                                 >  2:3808057/1‑150 (MQ=255)
                                                                                      ggggACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTCTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCa                                                                 >  1:3780453/1‑150 (MQ=255)
                                                                                        ggACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATc                                                               <  1:1898939/150‑1 (MQ=255)
                                                                                           cAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTCTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGc                                                            >  2:493329/1‑150 (MQ=255)
                                                                                             gTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTg                                                          >  1:1045886/1‑150 (MQ=255)
                                                                                                aGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAg                                                       >  2:2264226/1‑150 (MQ=255)
                                                                                                aGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAg                                                       >  2:36835/1‑150 (MQ=255)
                                                                                                aGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAg                                                       >  2:2466274/1‑150 (MQ=255)
                                                                                                    cATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCa                                                   >  1:2565465/1‑150 (MQ=255)
                                                                                                     aTTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTAGTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCaa                                                  >  2:379325/1‑150 (MQ=255)
                                                                                                      tttGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAAc                                                 <  1:669620/150‑1 (MQ=255)
                                                                                                            gTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTAGTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTctg                                           >  1:3743679/1‑148 (MQ=255)
                                                                                                             ttAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGgc                                          >  2:2137944/1‑150 (MQ=255)
                                                                                                                 aaaTTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAg                                      >  2:1842853/1‑150 (MQ=255)
                                                                                                                  aaTTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGt                                     >  2:4715877/1‑150 (MQ=255)
                                                                                                                  aaTTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTAACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGt                                     >  2:3529182/1‑150 (MQ=255)
                                                                                                                   aTTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGTa                                    <  1:4416715/150‑1 (MQ=255)
                                                                                                                   aTTTAGCTAAACACATATGAAGTTTCAGATGTGTCTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGTa                                    <  1:1623872/150‑1 (MQ=255)
                                                                                                                      tAGCTAAACACATATGAATTTTCAGATGTGTCTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGTACTg                                 >  1:3175515/1‑150 (MQ=255)
                                                                                                                       aGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGTACTgg                                >  1:32677/1‑150 (MQ=255)
                                                                                                                        gCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGTACTGGc                               >  2:4172002/1‑150 (MQ=255)
                                                                                                                         cTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGTAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGTACTGGCg                              >  2:1426806/1‑150 (MQ=255)
                                                                                                                         cTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGTACTGGCg                              <  2:2215506/150‑1 (MQ=255)
                                                                                                                         cTAAACACATATGAATTTTCAGATGTGTCTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGTACTGGCg                              >  2:4702523/1‑150 (MQ=255)
                                                                                                                         cTAAACACAGATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGTACTGGCg                              <  1:3368235/150‑1 (MQ=255)
                                                                                                                          tAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGTACTGGCgt                             >  1:198939/1‑150 (MQ=255)
                                                                                                                            aaCACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGTACTGGCgtta                           >  1:1195146/1‑150 (MQ=255)
                                                                                                                            aaCACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGTACTGGCgtta                           >  2:1330151/1‑150 (MQ=255)
                                                                                                                            aaCACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGTACTGGCgtta                           >  1:3239645/1‑150 (MQ=255)
                                                                                                                             acacATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGTACTGGCgttag                          <  1:2354197/150‑1 (MQ=255)
                                                                                                                              cacaTATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGTACTGGCgttagt                         >  1:544789/1‑150 (MQ=255)
                                                                                                                              cacaTATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGTACTGGCgttagt                         <  1:2757475/150‑1 (MQ=255)
                                                                                                                               acaTATGAATTTTCAGATGTGTCTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGTACTGGCgttagtt                        >  1:1014562/1‑150 (MQ=255)
                                                                                                                                  tatGAATTTTCAGATGTGTTTTATCCGGGAGGCAATATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGTACAGGCGTTAGTTACg                     >  2:563057/1‑150 (MQ=255)
                                                                                                                                     gAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGTACTGGCGTTAGTTACGGGc                  <  1:120348/150‑1 (MQ=255)
                                                                                                                                      aaTTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGTACTGGCGTTAGTTACGGGCg                 >  2:1011680/1‑150 (MQ=255)
                                                                                                                                      aaTTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGTACTGGCGTTAGTTACGGGCg                 >  2:287762/1‑150 (MQ=255)
                                                                                                                                         tttCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGTACTGGCGTTAGTTACGGGCGTgg              >  2:4827809/1‑150 (MQ=255)
                                                                                                                                          ttCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGTACTGGCGTTAGTTACGGGCGTGGt             >  1:3455373/1‑150 (MQ=255)
                                                                                                                                          tacAGATGTGTATAATCCGGGAGGCATTATGTAACTGGCAGGCGCTAACTATGTCCTGACCATCGTATTGGTTATCAGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCACGGCGCAGTACTGGCGTTAGTTACGGGCGTGGt             >  2:1370815/3‑150 (MQ=255)
                                                                                                                                            cAGATGTGTCTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGTACTGGCGTTAGTTACGGGCGTGGTcc           <  1:4833718/150‑1 (MQ=255)
                                                                                                                                             aGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGTACTGGCGTTAGTTACGGGCGTGGTCCa          >  2:2537963/1‑150 (MQ=255)
                                                                                                                                               aTGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGTACTGGCGTTAGTTACGGGCGTGGTCCATc        >  2:307472/1‑150 (MQ=255)
                                                                                                                                               aTGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGTACTGGCGTTAGTTACGGGCGTGGTCCATc        >  1:208271/1‑150 (MQ=255)
                                                                                                                                               aTGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGTACTGGCGTTAGTTACGGGCGTGGTCCATc        <  1:4624421/150‑1 (MQ=255)
                                                                                                                                                  tgtTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGTACTGGCGTTAGTTACGGGCGTGGTCCATCCgg     >  1:2145018/1‑150 (MQ=255)
                                                                                                                                                    ttttATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGTACTGGCGTTAGTTACGGGCGTGGTCCATCCGGGt   >  1:811085/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGGCGAAAATTATTGATGAGGCCTGGCGATGTGAACAGGAAAAGGTTCAGGCGATTGTCCGCAACCTGGCTCGACAAAACTGTTCCGGGGACAGTAAGAACATTTGCAGTTAAAAATTTAGCTAAACACATATGAATTTTCAGATGTGTTTTATCCGGGAGGCATTATGTTACTGGCAGGCGCTATCTTTGTCCTGACCATCGTATTGGTTATCTGGCAGCCGAAAGGTTTAGGCATCGGCTGGAGTGCAACGCTCGGCGCAGTACTGGCGTTAGTTACGGGCGTGGTCCATCCGGGTG  >  CP009273/3642129‑3642427

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: