Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,079,951 G→T 31.2% D374Y (GAT→TAT)  efeB → deferrrochelatase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,079,9510GT31.2% 500.9 / 205.0 252D374Y (GAT→TAT) efeBdeferrrochelatase, periplasmic
Reads supporting (aligned to +/- strand):  ref base G (81/92);  new base T (42/37);  total (123/129)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 256 total aligned reads displayed.

TCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTATTGCGGGTT  >  CP009273/1079802‑1080100
                                                                                                                                                     |                                                                                                                                                     
   cgcTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATc                                                                                                                                                    >  2:1645136/1‑150 (MQ=255)
     cTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACTATCTg                                                                                                                                                  <  1:788626/150‑1 (MQ=255)
     cTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACTATCTg                                                                                                                                                  <  1:788747/150‑1 (MQ=255)
      tGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTgg                                                                                                                                                 <  1:2035992/150‑1 (MQ=255)
            gCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAg                                                                                                                                           <  1:578212/150‑1 (MQ=255)
            gCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAg                                                                                                                                           <  1:1970991/150‑1 (MQ=255)
            gCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAg                                                                                                                                           <  1:1971073/150‑1 (MQ=255)
                 atCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTc                                                                                                                                      >  2:4613538/1‑150 (MQ=255)
                  tCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTcc                                                                                                                                     <  1:2068997/150‑1 (MQ=255)
                    cGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTg                                                                                                                                   >  1:4443327/1‑150 (MQ=255)
                    cGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTg                                                                                                                                   >  1:3012912/1‑150 (MQ=255)
                    cGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTg                                                                                                                                   <  2:2695020/150‑1 (MQ=255)
                     ggCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGa                                                                                                                                  <  1:4123702/150‑1 (MQ=255)
                       cTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACa                                                                                                                                >  2:3200663/1‑150 (MQ=255)
                          gCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTa                                                                                                                             >  1:1915564/1‑150 (MQ=255)
                          gCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTa                                                                                                                             >  2:3111306/1‑150 (MQ=255)
                           cGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTAc                                                                                                                            >  2:1490631/1‑150 (MQ=255)
                           cGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTAc                                                                                                                            <  2:3004297/150‑1 (MQ=255)
                           cGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTAc                                                                                                                            <  2:2319395/150‑1 (MQ=255)
                            gAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACa                                                                                                                           <  2:463239/150‑1 (MQ=255)
                                ccccGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACaaaaa                                                                                                                       >  1:1368062/1‑150 (MQ=255)
                                  ccGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAg                                                                                                                     >  1:414917/1‑150 (MQ=255)
                                     cACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCt                                                                                                                  >  1:733827/1‑150 (MQ=255)
                                      aCGGCGGAGTGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTc                                                                                                                 >  1:1800761/1‑150 (MQ=255)
                                       cggcggAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCa                                                                                                                <  2:3936392/150‑1 (MQ=255)
                                          cggAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATg                                                                                                             <  2:2334805/150‑1 (MQ=255)
                                            gagaGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGc                                                                                                           <  2:3457477/150‑1 (MQ=255)
                                             agagTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCg                                                                                                          >  1:640053/1‑150 (MQ=255)
                                             agagTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCg                                                                                                          >  1:4567535/1‑150 (MQ=255)
                                                gtgagtCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAg                                                                                                       <  2:3755763/150‑1 (MQ=255)
                                                  gagtCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAgcg                                                                                                     <  1:1467886/150‑1 (MQ=255)
                                                   agtCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAgcgc                                                                                                    <  1:39715/150‑1 (MQ=255)
                                                     tCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTg                                                                                                  <  2:2445614/150‑1 (MQ=255)
                                                      ccAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTgg                                                                                                 <  2:3058547/150‑1 (MQ=255)
                                                          ccTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTggagga                                                                                             >  1:1138176/1‑150 (MQ=255)
                                                            tGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAAt                                                                                           >  1:259719/1‑150 (MQ=255)
                                                               tgctgcGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACg                                                                                        >  2:3948172/1‑150 (MQ=255)
                                                                      tcgtGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAAcc                                                                                 <  2:3366162/150‑1 (MQ=255)
                                                                       cgtGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCt                                                                                <  2:3872458/150‑1 (MQ=255)
                                                                        gtGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCTCTGGAGGAATACGTTAAACCTa                                                                               <  1:2096191/150‑1 (MQ=255)
                                                                        gtGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTa                                                                               >  1:3705953/1‑150 (MQ=255)
                                                                        gtGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTa                                                                               >  1:1435828/1‑150 (MQ=255)
                                                                         tGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTAt                                                                              >  1:1709365/1‑150 (MQ=255)
                                                                          ggCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATc                                                                             <  1:2197514/150‑1 (MQ=255)
                                                                             tACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATcggc                                                                          <  2:4607062/150‑1 (MQ=255)
                                                                                aGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATcggcggc                                                                       <  2:2816883/150‑1 (MQ=255)
                                                                                  ttATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATcggcggcgg                                                                     <  2:297772/150‑1 (MQ=255)
                                                                                   tATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGt                                                                    <  1:3169569/150‑1 (MQ=255)
                                                                                     ttCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTa                                                                  <  2:1266790/150‑1 (MQ=255)
                                                                                      tCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTAt                                                                 >  2:1853105/1‑150 (MQ=255)
                                                                                      tCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTAt                                                                 >  2:1852925/1‑150 (MQ=255)
                                                                                      tCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTAt                                                                 >  2:1967759/1‑150 (MQ=255)
                                                                                        aCTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTAttt                                                               <  2:4135433/150‑1 (MQ=255)
                                                                                        aCTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTAttt                                                               >  2:2678738/1‑150 (MQ=255)
                                                                                         cTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTAtttt                                                              <  1:1625480/150‑1 (MQ=255)
                                                                                          tGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTAttttt                                                             >  1:3096931/1‑150 (MQ=255)
                                                                                           gggCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTAtttttt                                                            <  1:2939163/150‑1 (MQ=255)
                                                                                           gggCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTAtttttt                                                            >  2:4420742/1‑150 (MQ=255)
                                                                                            ggCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTAttttttt                                                           <  2:1737307/150‑1 (MQ=255)
                                                                                             gCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTg                                                          <  1:4613538/150‑1 (MQ=255)
                                                                                             gCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATATGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTg                                                          <  1:2126872/150‑1 (MQ=255)
                                                                                              cGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTgc                                                         <  2:4348423/150‑1 (MQ=255)
                                                                                               gTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTgcg                                                        <  2:1554612/150‑1 (MQ=255)
                                                                                                tCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTgcgc                                                       >  2:2192755/1‑150 (MQ=255)
                                                                                                tCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTgcgc                                                       >  2:3919270/1‑150 (MQ=255)
                                                                                                tCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTgcgc                                                       >  2:2492345/1‑150 (MQ=255)
                                                                                                  aCCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTg                                                     <  2:1015148/150‑1 (MQ=255)
                                                                                                  aCCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTg                                                     <  2:3301839/150‑1 (MQ=255)
                                                                                                    cAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGcc                                                   <  1:4619620/150‑1 (MQ=255)
                                                                                                     aaCTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCg                                                  >  2:1865425/1‑150 (MQ=255)
                                                                                                        tCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCgggg                                               >  1:2863433/1‑150 (MQ=255)
                                                                                                        tCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCgggg                                               >  1:4381829/1‑150 (MQ=255)
                                                                                                            ggCAACTGGATATGGGGTTGCTGTTAGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGa                                           >  2:1086952/1‑150 (MQ=255)
                                                                                                               aaCTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAgg                                        <  2:2933361/150‑1 (MQ=255)
                                                                                                                aCTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAAGa                                       <  1:903575/150‑1 (MQ=255)
                                                                                                                    gATATGGGGTTGCTGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGAcgcg                                   >  2:1977885/1‑150 (MQ=255)
                                                                                                                      tatGGGGTTGCTGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGaa                                 <  1:4213454/150‑1 (MQ=255)
                                                                                                                       atGGGGTTGCTGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAAc                                <  1:911930/150‑1 (MQ=255)
                                                                                                                       atGGGGTTGCTGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAAc                                <  1:3951645/150‑1 (MQ=255)
                                                                                                                         ggggTTGCTGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGa                              >  1:2192678/1‑150 (MQ=255)
                                                                                                                            gTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTTCGCTGCCGGGGGTGAAGGACGCGAACGatta                           >  1:1290894/1‑150 (MQ=255)
                                                                                                                             ttGCTGTTTGTCTGCTACCAACAATATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGattat                          <  2:1243504/150‑1 (MQ=255)
                                                                                                                              tGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGattatt                         >  2:2408831/1‑150 (MQ=255)
                                                                                                                              tGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGattatt                         >  1:3833438/1‑150 (MQ=255)
                                                                                                                               gCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCCGCCGGGGGTGAAGGACGCGAACGATTAttt                        <  1:3892570/150‑1 (MQ=255)
                                                                                                                                cTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAGCCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTc                       >  2:3174416/1‑150 (MQ=255)
                                                                                                                                 tGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCg                      <  2:2266499/150‑1 (MQ=255)
                                                                                                                                  gTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCgg                     >  1:241520/1‑150 (MQ=255)
                                                                                                                                     tGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAg                  >  2:677104/1‑150 (MQ=255)
                                                                                                                                      gTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAgc                 <  1:2202080/150‑1 (MQ=255)
                                                                                                                                         tGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGcgcg              <  2:4395751/150‑1 (MQ=255)
                                                                                                                                           cTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGtt            <  2:3619180/150‑1 (MQ=255)
                                                                                                                                              ccAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTAtt         <  1:4124986/150‑1 (MQ=255)
                                                                                                                                              ccAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTAtt         <  2:3818888/150‑1 (MQ=255)
                                                                                                                                               cAACACTATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTATTg        <  1:3467746/150‑1 (MQ=255)
                                                                                                                                               cAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTATTg        <  2:490345/150‑1 (MQ=255)
                                                                                                                                                aaCACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTATTGc       >  2:4154294/1‑150 (MQ=255)
                                                                                                                                                aaCACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTATTGc       >  2:3338136/1‑150 (MQ=255)
                                                                                                                                                 acacGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCGGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTATTGCg      <  2:4246255/150‑1 (MQ=255)
                                                                                                                                                  cacGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTATTGCgg     <  2:145833/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTATTGCGGGTT  >  CP009273/1079802‑1080100

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: