Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,577,849 G→A 11.3% intergenic (+180/‑270) yrhA → / → yrhB pseudogene, interrupted by IS1E/stable heat shock chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,577,8490GA11.3% 852.9 / 72.8 266intergenic (+180/‑270)yrhA/yrhBpseudogene, interrupted by IS1E/stable heat shock chaperone
Reads supporting (aligned to +/- strand):  ref base G (114/122);  new base A (4/26);  total (118/148)
Fisher's exact test for biased strand distribution p-value = 3.01e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 271 total aligned reads displayed.

CCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATC  >  CP009273/3577700‑3577998
                                                                                                                                                     |                                                                                                                                                     
 cATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGc                                                                                                                                                      >  2:2269001/1‑150 (MQ=255)
 cATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGc                                                                                                                                                      >  1:139007/1‑150 (MQ=255)
 cATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGc                                                                                                                                                      <  1:3108462/150‑1 (MQ=255)
   tGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTa                                                                                                                                                    <  1:2859050/150‑1 (MQ=255)
   tGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTa                                                                                                                                                    >  2:1037183/1‑150 (MQ=255)
       aTCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGACTatgat                                                                                                                                                <  2:760032/150‑1 (MQ=255)
          tttttAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATgattat                                                                                                                                             <  1:3835443/150‑1 (MQ=255)
           ttttAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGAtt                                                                                                                                            <  2:2877294/150‑1 (MQ=255)
           ttttAGAGATTACCAGCAGAAAATTCCTTGNTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGAtt                                                                                                                                            >  1:3856322/1‑150 (MQ=255)
            tttAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGACTATGATGAttt                                                                                                                                           <  2:860588/150‑1 (MQ=255)
             ttAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGAtttt                                                                                                                                          <  2:2416502/150‑1 (MQ=255)
               agagATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTaa                                                                                                                                        >  2:2256650/1‑150 (MQ=255)
               agagATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTaa                                                                                                                                        >  1:1486218/1‑150 (MQ=255)
               agagATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTaa                                                                                                                                        >  1:163071/1‑150 (MQ=255)
               agagATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTaa                                                                                                                                        >  1:3765117/1‑150 (MQ=255)
                   aTTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGAc                                                                                                                                    >  1:2901806/1‑150 (MQ=255)
                   aTTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGACTATGATGATTTTAATGAc                                                                                                                                    <  1:883838/150‑1 (MQ=255)
                    ttACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACa                                                                                                                                   <  1:2735967/150‑1 (MQ=255)
                         agcagAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTg                                                                                                                              <  1:2211289/150‑1 (MQ=255)
                         agcagAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTg                                                                                                                              <  2:1722607/150‑1 (MQ=255)
                           cagAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGat                                                                                                                            <  1:3143344/150‑1 (MQ=255)
                             gAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGattt                                                                                                                          <  2:1597900/150‑1 (MQ=255)
                             gAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGattt                                                                                                                          <  2:1597936/150‑1 (MQ=255)
                               aaaTTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGatttat                                                                                                                        >  2:1403056/1‑150 (MQ=255)
                               aaaTTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGatttat                                                                                                                        >  2:4057745/1‑150 (MQ=255)
                                 aTTCCTTGTTCGAAATTTAGCTGTGCAATCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATTGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGatttattt                                                                                                                      <  1:2186078/150‑1 (MQ=255)
                                  ttCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGAtttattta                                                                                                                     >  1:485673/1‑150 (MQ=255)
                                   tCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAg                                                                                                                    <  1:4134130/150‑1 (MQ=255)
                                    ccTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGa                                                                                                                   <  2:2151839/150‑1 (MQ=255)
                                           cGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGACTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTAc                                                                                                            <  1:977543/150‑1 (MQ=255)
                                            gAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTAcc                                                                                                           >  1:1227081/1‑150 (MQ=255)
                                              aaTTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCg                                                                                                         >  2:435947/1‑150 (MQ=255)
                                              aaTTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCg                                                                                                         <  2:2005106/150‑1 (MQ=255)
                                                 ttAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAAt                                                                                                      <  2:3995737/150‑1 (MQ=255)
                                                  tAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATc                                                                                                     <  1:3206688/150‑1 (MQ=255)
                                                   aGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCa                                                                                                    >  1:2299534/1‑150 (MQ=255)
                                                   aGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGACTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCa                                                                                                    <  2:655174/150‑1 (MQ=255)
                                                     cAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATa                                                                                                  <  1:1513054/150‑1 (MQ=255)
                                                         gCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGACTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGAt                                                                                              <  2:2396663/150‑1 (MQ=255)
                                                            aGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCt                                                                                           >  1:1054903/1‑150 (MQ=255)
                                                            aGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGACTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCt                                                                                           >  1:609464/1‑150 (MQ=255)
                                                                 tGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGatat                                                                                      >  1:2437623/1‑150 (MQ=255)
                                                                  ggAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATc                                                                                     >  2:2987340/1‑150 (MQ=255)
                                                                  ggAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATc                                                                                     <  1:247191/150‑1 (MQ=255)
                                                                   gAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCt                                                                                    <  1:2356639/150‑1 (MQ=255)
                                                                   gAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCt                                                                                    <  2:2171646/150‑1 (MQ=255)
                                                                   gAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCt                                                                                    >  2:3890149/1‑150 (MQ=255)
                                                                    aaaaaCTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGACTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCtt                                                                                   <  1:3910809/150‑1 (MQ=255)
                                                                     aaaaCTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTa                                                                                  <  2:3501124/150‑1 (MQ=255)
                                                                         cTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATa                                                                              >  2:809251/1‑150 (MQ=255)
                                                                              aaCGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGaaa                                                                         <  2:1320763/150‑1 (MQ=255)
                                                                              aaCGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGaaa                                                                         >  2:3080640/1‑150 (MQ=255)
                                                                              aaCGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGaaa                                                                         <  2:1815388/150‑1 (MQ=255)
                                                                               aCGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGaaaa                                                                        <  2:1766948/150‑1 (MQ=255)
                                                                                 ggAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGACTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATc                                                                      <  2:4094025/150‑1 (MQ=255)
                                                                                   aaCTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATccc                                                                    <  2:2759165/150‑1 (MQ=255)
                                                                                      tACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGtt                                                                 <  1:1111055/150‑1 (MQ=255)
                                                                                      tACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGtt                                                                 <  2:1232985/150‑1 (MQ=255)
                                                                                       acacTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTc                                                                <  1:2931750/150‑1 (MQ=255)
                                                                                        cacTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCg                                                               <  2:1693337/150‑1 (MQ=255)
                                                                                         acTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGa                                                              >  1:1285911/1‑150 (MQ=255)
                                                                                         acTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGa                                                              <  2:1541000/150‑1 (MQ=255)
                                                                                         acTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGa                                                              >  2:614690/1‑150 (MQ=255)
                                                                                           tCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACg                                                            >  1:4003020/1‑150 (MQ=255)
                                                                                           tCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACg                                                            <  1:961861/150‑1 (MQ=255)
                                                                                            ccTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGt                                                           >  1:3035694/1‑150 (MQ=255)
                                                                                              ttAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTc                                                         <  1:3288214/150‑1 (MQ=255)
                                                                                               tAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCt                                                        >  1:2898335/1‑150 (MQ=255)
                                                                                                 aGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGACTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGa                                                      <  2:1183385/150‑1 (MQ=255)
                                                                                                    gagaTGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGaaaa                                                   >  1:3539733/1‑150 (MQ=255)
                                                                                                     agaTGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGaaaaa                                                  >  2:3127155/1‑150 (MQ=255)
                                                                                                       aTGGACCGAAGACATGTAGTCTCTATATTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGa                                                <  1:1930993/150‑1 (MQ=255)
                                                                                                               aaGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATTGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATaa                                        <  1:760112/150‑1 (MQ=255)
                                                                                                               aaGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATaa                                        <  1:3360962/150‑1 (MQ=255)
                                                                                                                  aCATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAACATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAATTCATatatt                                     <  2:2680977/150‑5 (MQ=255)
                                                                                                                   cATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAata                                    >  2:1787994/1‑150 (MQ=255)
                                                                                                                   cATGTACTCTCTATATTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTAATTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAata                                    <  2:3604120/150‑1 (MQ=255)
                                                                                                                    aTGTAGTCTCTATTTTAGCCGAAATAGCTAATGACTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAatat                                   >  2:1870645/1‑150 (MQ=255)
                                                                                                                     tGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATAtt                                  >  2:487047/1‑150 (MQ=255)
                                                                                                                     tGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATAtt                                  >  1:1853819/1‑150 (MQ=255)
                                                                                                                     tGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATAtt                                  >  1:3896103/1‑150 (MQ=255)
                                                                                                                      gTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTa                                 >  1:798333/1‑150 (MQ=255)
                                                                                                                      gTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTa                                 >  1:2413247/1‑150 (MQ=255)
                                                                                                                        aGTCTCTATTTTAGCCGAAATAGCTAATGACTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGc                               >  1:1678619/1‑150 (MQ=255)
                                                                                                                            tctATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGAt                           <  2:2936183/150‑1 (MQ=255)
                                                                                                                            tctATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGAt                           <  2:1972746/150‑1 (MQ=255)
                                                                                                                              tATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATAc                         <  2:1524849/150‑1 (MQ=255)
                                                                                                                               aTTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATACt                        >  1:684241/1‑150 (MQ=255)
                                                                                                                                ttttAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATACtt                       <  2:3772567/150‑1 (MQ=255)
                                                                                                                                    aGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTAttt                   <  2:865345/150‑1 (MQ=255)
                                                                                                                                     gCCGAAATAGCTAATGACTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTc                  <  2:2462174/150‑1 (MQ=255)
                                                                                                                                      ccGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCg                 <  2:2652444/150‑1 (MQ=255)
                                                                                                                                      ccGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCg                 <  2:2559609/150‑1 (MQ=255)
                                                                                                                                          aaTAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGAc             >  2:561342/1‑150 (MQ=255)
                                                                                                                                              gCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACtata         >  2:1472273/1‑150 (MQ=255)
                                                                                                                                                tAATGACTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATc       <  1:1032179/150‑1 (MQ=255)
                                                                                                                                                tAATGACTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCACGTTCGACGTTCTGATAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATc       <  1:1032235/150‑1 (MQ=255)
                                                                                                                                                   tGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTa    <  1:556954/150‑1 (MQ=255)
                                                                                                                                                   tGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTa    <  1:2128321/150‑1 (MQ=255)
                                                                                                                                                   tGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTa    >  1:1315400/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTGAGATGGACCGAAGACATGTAGTCTCTATTTTAGCCGAAATAGCTAATGGCTATGATGATTTTAATGACATGTTGATTTATTTAGAATTTTACCCGAATCATAAGATTCTGATATCTTAAATAAGAAAATCCCGTTCGACGTTCTGAAAAAGAAATCATAAAATATTAGCAGATACTTATTTCGTGACTATATCTTATC  >  CP009273/3577700‑3577998

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: