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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,079,9510GT 7.4% 1021.2 / 36.6 270D374Y (GAT→TAT) efeBdeferrrochelatase, periplasmic
Reads supporting (aligned to +/- strand):  ref base G (122/128);  new base T (8/12);  total (130/140)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 275 total aligned reads displayed.

TCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTATTGCGGGTT  >  CP009273/1079802‑1080100
                                                                                                                                                     |                                                                                                                                                     
tCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACg                                                                                                                                                       <  1:3932128/150‑1 (MQ=255)
        gACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGaa                                                                                                                                               >  2:2342607/1‑150 (MQ=255)
          cAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGaaaa                                                                                                                                             <  2:2795125/150‑1 (MQ=255)
             ccATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAgg                                                                                                                                          >  1:19678/1‑150 (MQ=255)
             ccATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAgg                                                                                                                                          >  1:2844277/1‑150 (MQ=255)
               atatCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCt                                                                                                                                        >  1:4527850/1‑150 (MQ=255)
               atatCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCt                                                                                                                                        >  1:3998140/1‑150 (MQ=255)
               atatCCGGCTGGAGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCt                                                                                                                                        <  2:516534/150‑1 (MQ=255)
                 atCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTc                                                                                                                                      >  2:247840/1‑150 (MQ=255)
                 atCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTc                                                                                                                                      <  2:276579/150‑1 (MQ=255)
                  tCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTcc                                                                                                                                     <  1:1263504/150‑1 (MQ=255)
                  tCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTcc                                                                                                                                     <  2:2660839/150‑1 (MQ=255)
                   ccGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCt                                                                                                                                    <  1:3206627/150‑1 (MQ=255)
                   ccGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCt                                                                                                                                    <  1:539158/150‑1 (MQ=255)
                    cgtcTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTg                                                                                                                                   >  2:1773176/4‑150 (MQ=255)
                    cGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTg                                                                                                                                   >  1:2413204/1‑150 (MQ=255)
                      gCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACGGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGAc                                                                                                                                 <  2:2216496/150‑1 (MQ=255)
                              aTCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACaaa                                                                                                                         <  1:3765666/150‑1 (MQ=255)
                                  ccGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAg                                                                                                                     <  2:3731041/150‑1 (MQ=255)
                                  ccGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAg                                                                                                                     <  2:1441652/150‑1 (MQ=255)
                                     cACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCt                                                                                                                  <  2:3637120/150‑1 (MQ=255)
                                          cggAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATg                                                                                                             <  1:2173401/150‑1 (MQ=255)
                                           ggAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATgg                                                                                                            >  1:2281/1‑150 (MQ=255)
                                           ggAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATgg                                                                                                            >  1:2199/1‑150 (MQ=255)
                                              gagtgagtCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGa                                                                                                         <  2:2165354/150‑1 (MQ=255)
                                               agtgagtCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGaa                                                                                                        >  2:4072482/1‑150 (MQ=255)
                                               agtgagtCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGaa                                                                                                        <  2:3771722/150‑1 (MQ=255)
                                                   agtCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAgcgc                                                                                                    <  1:1745895/150‑1 (MQ=255)
                                                    gtCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCt                                                                                                   <  2:1268650/150‑1 (MQ=255)
                                                    gtCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCt                                                                                                   <  2:4097805/150‑1 (MQ=255)
                                                       cAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCATCTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTgga                                                                                                <  1:508764/150‑1 (MQ=255)
                                                        aGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTggag                                                                                               >  2:729948/1‑150 (MQ=255)
                                                           cTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTAACAACACGATCTGTAAAAAGGCTTCCTGACAGTACAAAAAAGGATCAATGGCGAAGCGCTGGAGGaa                                                                                            >  1:2043425/1‑150 (MQ=255)
                                                            tGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAAt                                                                                           >  2:1992875/1‑150 (MQ=255)
                                                            tGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCAACAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAAt                                                                                           <  2:4511063/150‑1 (MQ=255)
                                                                 ctgcGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGtt                                                                                      <  1:2666875/150‑1 (MQ=255)
                                                                  tgcGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTa                                                                                     >  1:3455543/1‑150 (MQ=255)
                                                                  tgcGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTa                                                                                     <  2:786319/150‑1 (MQ=255)
                                                                       cgtGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCt                                                                                <  2:1521586/150‑1 (MQ=255)
                                                                        gtGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTa                                                                               <  2:2058047/150‑1 (MQ=255)
                                                                          ggCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATc                                                                             <  1:598324/150‑1 (MQ=255)
                                                                            cTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATcgg                                                                           <  1:901238/150‑1 (MQ=255)
                                                                              aCAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATcggcg                                                                         <  1:1822198/150‑1 (MQ=255)
                                                                               cAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATcggcgg                                                                        >  2:3570187/1‑150 (MQ=255)
                                                                                aGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATcggcggc                                                                       >  1:2127919/1‑150 (MQ=255)
                                                                                  ttATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATcggcggcgg                                                                     >  1:1369790/1‑150 (MQ=255)
                                                                                  ttATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATcggcggcgg                                                                     >  1:2732661/1‑150 (MQ=255)
                                                                                   tATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGt                                                                    >  2:2045047/1‑150 (MQ=255)
                                                                                     ttCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTa                                                                  >  1:811411/1‑150 (MQ=255)
                                                                                        aCTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTAttt                                                               >  1:4451914/1‑150 (MQ=255)
                                                                                         cTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTAtttt                                                              <  1:2321878/150‑1 (MQ=255)
                                                                                          tGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTAttttt                                                             >  1:2019774/1‑150 (MQ=255)
                                                                                          tGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTAttttt                                                             >  1:4346439/1‑150 (MQ=255)
                                                                                           gggCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTAtttttt                                                            <  2:4623627/150‑1 (MQ=255)
                                                                                            ggCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTAttttttt                                                           <  2:1121883/150‑1 (MQ=255)
                                                                                            ggCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTAttttttt                                                           <  2:1885054/150‑1 (MQ=255)
                                                                                                tCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTgcgc                                                       >  2:3380354/1‑150 (MQ=255)
                                                                                                 cACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCt                                                      <  2:1260214/150‑1 (MQ=255)
                                                                                                  aCCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTg                                                     <  2:1368664/150‑1 (MQ=255)
                                                                                                  aCCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTg                                                     <  2:3183503/150‑1 (MQ=255)
                                                                                                   ccAACTCCGGGCAACTGGATATGGGGTTGATGTTTGTCTGCTACCAACACGATATGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATTGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTTCGCTGc                                                    >  2:2847819/1‑150 (MQ=255)
                                                                                                    cAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGTCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGAGCTGca                                                   >  1:3488471/1‑149 (MQ=255)
                                                                                                    cAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGcc                                                   >  2:2391650/1‑150 (MQ=255)
                                                                                                    cAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGcc                                                   <  2:3202233/150‑1 (MQ=255)
                                                                                                    cAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGcc                                                   >  1:4628352/1‑150 (MQ=255)
                                                                                                     aaCTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCg                                                  >  2:816958/1‑150 (MQ=255)
                                                                                                           ggggAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTCTTTTTTTGCGCTGCCGGGGGTg                                            >  2:308378/1‑150 (MQ=255)
                                                                                                           gggCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTg                                            <  1:3734000/150‑1 (MQ=255)
                                                                                                             gCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGaa                                          >  2:3239260/1‑150 (MQ=255)
                                                                                                                 cTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGAc                                      <  1:480384/150‑1 (MQ=255)
                                                                                                                  tGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGAcg                                     <  1:2899383/150‑1 (MQ=255)
                                                                                                                   ggATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGAcgc                                    >  2:864098/1‑150 (MQ=255)
                                                                                                                    gATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGAcgcg                                   >  2:674751/1‑150 (MQ=255)
                                                                                                                    gATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGAcgcg                                   >  2:161155/1‑150 (MQ=255)
                                                                                                                     atatGTGGTTGCTGTTTGTCTGCTACCAACACGATCGGGAAAAAGGCTTCCTGACAGTACAAAACAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGTTTATTTTTTTGCGCTGCCGGGGGTGAAGGAAGCGa                                  <  1:3240358/150‑1 (MQ=255)
                                                                                                                     atatGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGa                                  <  1:1210661/150‑1 (MQ=255)
                                                                                                                         ggggTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGa                              >  1:1560750/1‑150 (MQ=255)
                                                                                                                         ggggTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGa                              <  2:627922/150‑1 (MQ=255)
                                                                                                                                cTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTc                       <  2:3493655/150‑1 (MQ=255)
                                                                                                                                cTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTc                       >  2:4264370/1‑150 (MQ=255)
                                                                                                                                cTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTc                       <  1:3660895/150‑1 (MQ=255)
                                                                                                                                  gTTTGTCTGCTACCAACACTATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCgg                     >  1:692938/1‑150 (MQ=255)
                                                                                                                                     tGTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAg                  <  1:3979977/150‑1 (MQ=255)
                                                                                                                                      gTCTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAgc                 >  2:3321488/1‑150 (MQ=255)
                                                                                                                                        cTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAgcgc               >  1:3609040/1‑150 (MQ=255)
                                                                                                                                        cTGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAgcgc               >  1:561954/1‑150 (MQ=255)
                                                                                                                                         tGCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATGATTTCGGAAGcgcg              >  1:4463883/1‑150 (MQ=255)
                                                                                                                                          gCTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGt             <  2:539365/150‑1 (MQ=255)
                                                                                                                                           cTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGtt            <  2:4107693/150‑1 (MQ=255)
                                                                                                                                           cTACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGtt            <  2:1727543/150‑1 (MQ=255)
                                                                                                                                            tACCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTa           >  1:411274/1‑150 (MQ=255)
                                                                                                                                             aCCAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTAt          <  1:4301120/150‑1 (MQ=255)
                                                                                                                                              ccAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTAtt         <  1:435791/150‑1 (MQ=255)
                                                                                                                                              ccAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTAtt         >  2:1667759/1‑150 (MQ=255)
                                                                                                                                              ccAACACGATCTGGAAAAAGGATTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCAGGAGGAATACGTTAAACATATCGGCGGCGGATATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGGTAtt         >  2:392314/1‑150 (MQ=255)
                                                                                                                                               cAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGTGCGTTATTg        >  1:1691862/1‑150 (MQ=255)
                                                                                                                                               cAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTATTg        <  2:4616404/150‑1 (MQ=255)
                                                                                                                                               cAACACGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTATTg        >  2:428641/1‑150 (MQ=255)
                                                                                                                                                 acacGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTATTGCg      >  2:1518339/1‑150 (MQ=255)
                                                                                                                                                   acGATCTGGAAAAAGGCTTCCTGACAGTACAAAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTATTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTATTGCggg    <  2:4420681/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: