Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,590,930 T→A 7.1% intergenic (‑10/‑414) livK ← / → panM leucine transporter subunit/PanD autocleavage accelerator, panothenate synthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,590,9300TA 7.1% 1335.1 / 68.5 396intergenic (‑10/‑414)livK/panMleucine transporter subunit/PanD autocleavage accelerator, panothenate synthesis
Reads supporting (aligned to +/- strand):  ref base T (153/215);  new base A (15/13);  total (168/228)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

Only 100 of 397 total aligned reads displayed.

TTAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGATAATAGCGATTAAAATCCCTATTTTTCAGTCG  >  CP009273/3590781‑3591078
                                                                                                                                                     |                                                                                                                                                    
ttAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGAGAGTTTTCGCATTCCGTTTCATAGTCAAAAAt                                                                                                                                                      >  2:2736813/1‑150 (MQ=255)
 tAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATc                                                                                                                                                     <  1:3622183/150‑1 (MQ=255)
  aaaTTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATcc                                                                                                                                                    <  2:2366700/150‑1 (MQ=255)
  aaaTTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATcc                                                                                                                                                    <  2:240720/150‑1 (MQ=255)
     ttCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCa                                                                                                                                                 <  2:329185/150‑1 (MQ=255)
       ccATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCAtt                                                                                                                                               <  1:4450498/150‑1 (MQ=255)
       ccATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCAtt                                                                                                                                               >  1:2684361/1‑150 (MQ=255)
         atatCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCg                                                                                                                                             <  1:3029192/150‑1 (MQ=255)
         atatCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCg                                                                                                                                             <  1:2592205/150‑1 (MQ=255)
              gCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGAtg                                                                                                                                        >  2:1902131/1‑150 (MQ=255)
                cccACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGAtgtt                                                                                                                                      <  1:3487295/150‑1 (MQ=255)
                  cACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGAtgttgt                                                                                                                                    >  1:2136373/1‑150 (MQ=255)
                  cACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGAtgttgt                                                                                                                                    >  2:3426731/1‑150 (MQ=255)
                  cACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGAtgttgt                                                                                                                                    <  2:2808718/150‑1 (MQ=255)
                   aCTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGAtgttgtg                                                                                                                                   >  1:2716126/1‑150 (MQ=255)
                    cTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGAtgttgtgt                                                                                                                                  >  1:2400263/1‑150 (MQ=255)
                     tGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGAtgttgtgtt                                                                                                                                 <  1:922519/150‑1 (MQ=255)
                      gggCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGAtgttgtgttg                                                                                                                                <  2:4714360/150‑1 (MQ=255)
                       ggCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGc                                                                                                                               >  1:4939976/1‑150 (MQ=255)
                             cGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGtt                                                                                                                         >  1:2234546/1‑150 (MQ=255)
                             cGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGtt                                                                                                                         <  2:435241/150‑1 (MQ=255)
                             cGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGtt                                                                                                                         >  2:3909826/1‑150 (MQ=255)
                               ggCCGGACATCGCGCCGTCAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCAATTCGTGATGTTGTGTTGCTTTGtttt                                                                                                                       <  2:2595499/150‑1 (MQ=255)
                               ggCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGtttt                                                                                                                       >  1:3100258/1‑150 (MQ=255)
                                      cATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATgtgt                                                                                                                >  2:3046466/1‑150 (MQ=255)
                                       aTCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGtt                                                                                                               >  1:652213/1‑150 (MQ=255)
                                         cgcgCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTACGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTaa                                                                                                             <  1:2279717/150‑1 (MQ=255)
                                           cgcCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACa                                                                                                           <  1:1120818/150‑1 (MQ=255)
                                            gcCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACaa                                                                                                          <  1:4628128/150‑1 (MQ=255)
                                             ccGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACaaa                                                                                                         <  2:4806536/150‑1 (MQ=255)
                                              cGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAAt                                                                                                        <  2:1056052/150‑1 (MQ=255)
                                               gACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATc                                                                                                       >  1:15062/1‑150 (MQ=255)
                                                aCAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCa                                                                                                      >  2:1617866/1‑150 (MQ=255)
                                                    cGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACt                                                                                                  <  2:1249423/150‑1 (MQ=255)
                                                       cGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAAACCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGtt                                                                                               <  1:3522146/150‑1 (MQ=255)
                                                        gACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTc                                                                                              <  2:3100381/150‑1 (MQ=255)
                                                         aCTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCt                                                                                             >  2:1280574/1‑150 (MQ=255)
                                                          cTTTAAGATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCtt                                                                                            <  2:4420823/150‑1 (MQ=255)
                                                           tttAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCttt                                                                                           >  1:4522881/1‑150 (MQ=255)
                                                           tttAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCttt                                                                                           <  2:1754066/150‑1 (MQ=255)
                                                           tttAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCttt                                                                                           >  2:888851/1‑150 (MQ=255)
                                                                tatCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTta                                                                                      <  1:3827110/150‑1 (MQ=255)
                                                                       aGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCa                                                                               >  2:3236770/1‑150 (MQ=255)
                                                                           aTAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGt                                                                           <  1:1235727/150‑1 (MQ=255)
                                                                               cGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGttttt                                                                       <  2:3616797/150‑1 (MQ=255)
                                                                                ggTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTg                                                                      >  1:4264039/1‑150 (MQ=255)
                                                                                ggTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTg                                                                      >  1:1788367/1‑150 (MQ=255)
                                                                                 gtgtgtGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGc                                                                     <  1:4991421/150‑1 (MQ=255)
                                                                                  tgtgtgAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGcc                                                                    <  1:3420614/150‑1 (MQ=255)
                                                                                    tgtgAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTg                                                                  <  1:1008746/150‑1 (MQ=255)
                                                                                     gtgAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGt                                                                 >  1:1872041/1‑150 (MQ=255)
                                                                                     gtgAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGt                                                                 >  2:4964316/1‑150 (MQ=255)
                                                                                       gAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCt                                                               >  1:332134/1‑150 (MQ=255)
                                                                                       gAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTGTTATAATGCACTGTTTTTGCCTGTCt                                                               >  2:3661219/1‑150 (MQ=255)
                                                                                             gCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTaa                                                         <  1:1318414/150‑1 (MQ=255)
                                                                                               cAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAgg                                                       >  1:2830177/1‑150 (MQ=255)
                                                                                                aGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAggg                                                      <  2:2981063/150‑1 (MQ=255)
                                                                                                aGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAggg                                                      >  1:2526243/1‑150 (MQ=255)
                                                                                                 gTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAgggg                                                     <  2:986697/150‑1 (MQ=255)
                                                                                                    cAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTa                                                  >  2:5016399/1‑150 (MQ=255)
                                                                                                      atcatcCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAgc                                                >  2:1385520/1‑150 (MQ=255)
                                                                                                      atcatcCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAgc                                                <  2:1126303/150‑1 (MQ=255)
                                                                                                       tcatcCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAgcg                                               <  2:15236/150‑1 (MQ=255)
                                                                                                        catcCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAgcgc                                              >  1:2273109/1‑150 (MQ=255)
                                                                                                         atcCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCa                                             >  2:3702063/1‑150 (MQ=255)
                                                                                                         atcCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCa                                             >  1:4828289/1‑150 (MQ=255)
                                                                                                          tcCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAg                                            >  2:546964/1‑150 (MQ=255)
                                                                                                            ccTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTa                                          <  1:1617866/150‑1 (MQ=255)
                                                                                                            ccTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTa                                          >  2:1821119/1‑150 (MQ=255)
                                                                                                            ccTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTa                                          <  1:993880/150‑1 (MQ=255)
                                                                                                             cTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTAt                                         >  2:714307/1‑150 (MQ=255)
                                                                                                             cTGCGATGATAGTCTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTAt                                         <  2:898471/150‑1 (MQ=255)
                                                                                                                cGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTAtttt                                      <  1:5088047/150‑1 (MQ=255)
                                                                                                                  atgatAGATTTCGCAGTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTtgt                                    <  2:2873615/150‑1 (MQ=255)
                                                                                                                   tgatAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTtgtg                                   >  2:2523566/1‑150 (MQ=255)
                                                                                                                        gTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGataat                              <  2:1196211/150‑1 (MQ=255)
                                                                                                                         ttttCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGataata                             <  2:335882/150‑1 (MQ=255)
                                                                                                                         ttttCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGataata                             <  1:3928979/150‑1 (MQ=255)
                                                                                                                         ttttCGCATTCCGTTTCATAGTCAAAAAACCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGataata                             <  2:2240818/150‑1 (MQ=255)
                                                                                                                          tttCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGATAATAg                            <  2:4123588/150‑1 (MQ=255)
                                                                                                                          tttCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGATAATAg                            <  2:4607609/150‑1 (MQ=255)
                                                                                                                           ttCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGATAATAGc                           >  2:4551400/1‑150 (MQ=255)
                                                                                                                           ttCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGATAATAGc                           <  1:1737075/150‑1 (MQ=255)
                                                                                                                           ttCGCATTCCGTTTCATAGTCAAAAAACCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGATAATAGc                           >  1:4766699/1‑150 (MQ=255)
                                                                                                                             cGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGATAATAGCGa                         >  2:3515782/1‑150 (MQ=255)
                                                                                                                               cATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGATAATAGCGAtt                       <  2:3853891/150‑1 (MQ=255)
                                                                                                                                aTTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGATAATAGCGATTa                      <  1:1394550/150‑1 (MQ=255)
                                                                                                                                 ttCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGATAATAGCGATTaa                     <  1:545510/150‑1 (MQ=255)
                                                                                                                                 ttCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGATAATAGCGATTaa                     <  1:1989124/150‑1 (MQ=255)
                                                                                                                                  tCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGATAATAGCGATTaaa                    >  1:4199327/1‑150 (MQ=255)
                                                                                                                                    cGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGATAATAGCGATTAAAAt                  <  2:2170063/150‑1 (MQ=255)
                                                                                                                                    cGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGATAATAGCGATTAAAAt                  >  2:833780/1‑150 (MQ=255)
                                                                                                                                      tttCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGATAATAGCGATTAAAATcc                <  2:2655901/150‑1 (MQ=255)
                                                                                                                                       ttCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGATAATAGCGATTAAAATccc               <  2:2038515/150‑1 (MQ=255)
                                                                                                                                          aTAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGATAATAGCGATTAAAATCCCTAt            <  2:3970132/150‑1 (MQ=255)
                                                                                                                                             gTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGATAATAGCGATTAAAATCCCTAtttt         <  1:219698/150‑1 (MQ=255)
                                                                                                                                              tCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGATAATAGCGATTAAAATCCCTAttttt        <  2:4268208/150‑1 (MQ=255)
                                                                                                                                                 aaaaTCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGATAATAGCGATTAAAATCCCTATTTTTCAg     <  1:1313897/150‑1 (MQ=255)
                                                                                                                                                  aaaTCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGATAATAGCGATTAAAATCCCTATTTTTCAGt    <  1:3737788/150‑1 (MQ=255)
                                                                                                                                                    aTCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGATAATAGCGATTAAAATCCCTATTTTTCAGTCg  <  1:605371/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTAAATTCCATATCGCCCCACTGGGCAATCGGGCCGGACATCGCGCCGACAACGGCGACTTTAATATCGTCAGCCATAGCGGTGTGTGAAATTGCCAGTGCAATCATCCCTGCGATGATAGTTTTCGCATTCCGTTTCATAGTCAAAAATCCCCATTCGTGATGTTGTGTTGCTTTGTTTTTATGTGTTAACAAATCAGACTGTTCTTTTTTTATACTGCACTGTTTTTGCCTGTCTGATTAAGGGGTTAGCGCAGTATTTTGTGATAATAGCGATTAAAATCCCTATTTTTCAGTCG  >  CP009273/3590781‑3591078

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: