Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,566,392 A→T 6.3% intergenic (‑90/+272) gadB ← / ← pqqL glutamate decarboxylase B, PLP‑dependent/putative periplasmic M16 family zinc metalloendopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,566,3920AT 6.3% 682.7 / 26.5 195intergenic (‑90/+272)gadB/pqqLglutamate decarboxylase B, PLP‑dependent/putative periplasmic M16 family zinc metalloendopeptidase
Reads supporting (aligned to +/- strand):  ref base A (95/87);  new base T (2/10);  total (98/97)
Fisher's exact test for biased strand distribution p-value = 3.29e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 209 total aligned reads displayed.

GCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAA  >  CP009273/1566244‑1566541
                                                                                                                                                    |                                                                                                                                                     
gCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACaa                                                                                                                                                      >  2:863852/1‑150 (MQ=255)
 cACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACaaa                                                                                                                                                     <  1:1163342/150‑1 (MQ=255)
 cACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACaaa                                                                                                                                                     >  2:3934669/1‑150 (MQ=255)
 cACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACaaa                                                                                                                                                     >  2:2573965/1‑150 (MQ=255)
 cACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACaaa                                                                                                                                                     >  2:1658416/1‑150 (MQ=255)
    cAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATc                                                                                                                                                  <  2:3378525/150‑1 (MQ=255)
    aaaaaCTTGAATAGAGTAGTTCCGAACTTAAATCCGTTCCTTTCTTCTTATCCATTTTAAAATCCTTAAAATGATTGGATCGCATTAAAAAAGTAGTATTTATCGATAAAGTAAGCAATTGCAAAGGACTCGTGTTTAAATAACAAAATc                                                                                                                                                  <  1:2668778/149‑1 (MQ=255)
     aaaaCGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTACGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATcc                                                                                                                                                 >  2:3215885/1‑150 (MQ=255)
     aaaaCGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATcc                                                                                                                                                 >  1:1723220/1‑150 (MQ=255)
      aaaCGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCt                                                                                                                                                <  1:723481/150‑1 (MQ=255)
             aaTCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGgta                                                                                                                                         >  2:4071146/1‑150 (MQ=255)
                 gAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGttattta                                                                                                                                     >  2:3555387/1‑150 (MQ=255)
                  agtagtTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGttatttat                                                                                                                                    <  1:93901/150‑1 (MQ=255)
                  agtagtTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGttatttat                                                                                                                                    >  1:1977816/1‑150 (MQ=255)
                    tagtTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCg                                                                                                                                  >  2:1736097/1‑150 (MQ=255)
                     agtTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGt                                                                                                                                 <  2:2043846/150‑1 (MQ=255)
                      gtTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTg                                                                                                                                >  2:2824730/1‑150 (MQ=255)
                       ttCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTga                                                                                                                               >  2:2189503/1‑150 (MQ=255)
                          cGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGat                                                                                                                            <  1:3720642/150‑1 (MQ=255)
                             ccTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATAtt                                                                                                                         >  1:3208394/1‑150 (MQ=255)
                               ttAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTAc                                                                                                                       >  1:4242282/1‑150 (MQ=255)
                                tAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTAcg                                                                                                                      >  1:2902568/1‑150 (MQ=255)
                                  aaTCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTAcgcg                                                                                                                    >  1:4487407/1‑150 (MQ=255)
                                   aTCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGa                                                                                                                   <  1:239296/150‑1 (MQ=255)
                                      cGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGaata                                                                                                                >  1:840944/1‑150 (MQ=255)
                                         tACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAataata                                                                                                             >  2:1195189/1‑150 (MQ=255)
                                          aCTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAatat                                                                                                            >  1:1763116/1‑150 (MQ=255)
                                          aCTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAatat                                                                                                            <  1:1121269/150‑1 (MQ=255)
                                          aCTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAatat                                                                                                            <  1:4419431/150‑1 (MQ=255)
                                          aCTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAatat                                                                                                            <  2:4521787/150‑1 (MQ=255)
                                           cTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATAtt                                                                                                           <  2:3986343/150‑1 (MQ=255)
                                             tGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATAtttt                                                                                                         <  1:2444995/150‑1 (MQ=255)
                                             tGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATAtttt                                                                                                         >  2:23114/1‑150 (MQ=255)
                                               cttcttATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATAtttttt                                                                                                       >  2:3421927/1‑150 (MQ=255)
                                                  cttATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACTAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCAt                                                                                                    <  2:3047340/150‑1 (MQ=255)
                                                  cttATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCAt                                                                                                    >  1:3065012/1‑150 (MQ=255)
                                                    tATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTg                                                                                                  <  2:1451349/150‑1 (MQ=255)
                                                     aTCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGa                                                                                                 <  2:3862204/150‑1 (MQ=255)
                                                      tCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGaa                                                                                                >  1:1639726/1‑150 (MQ=255)
                                                       ccATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACGAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGaaa                                                                                               >  1:1470418/1‑150 (MQ=255)
                                                         aTTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGaaaaa                                                                                             <  1:2736702/150‑1 (MQ=255)
                                                          ttttAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAac                                                                                            <  2:375236/150‑1 (MQ=255)
                                                          ttttAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAac                                                                                            <  1:2685160/150‑1 (MQ=255)
                                                           tttAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGaaaaaaa                                                                                           >  1:4195915/1‑148 (MQ=255)
                                                             tAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAacaat                                                                                         <  1:2610861/150‑1 (MQ=255)
                                                             tAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAacaat                                                                                         >  1:464757/1‑150 (MQ=255)
                                                             tAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAacaat                                                                                         <  1:1207446/150‑1 (MQ=255)
                                                             tAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAacaat                                                                                         <  1:3933514/150‑1 (MQ=255)
                                                              aaaCTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAacaata                                                                                        <  1:432433/150‑1 (MQ=255)
                                                                aCTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATCTTACGCGAATAATATTTTTTCATTGAAAAacaataca                                                                                      >  1:568663/1‑150 (MQ=255)
                                                                     ttAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATg                                                                                 <  2:456548/150‑1 (MQ=255)
                                                                     ttAAAATGATTGGATCGAATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATg                                                                                 <  1:2291706/150‑1 (MQ=255)
                                                                       aaaaTGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACTAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGaa                                                                               >  2:2475906/1‑150 (MQ=255)
                                                                       aaaaTGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGaa                                                                               <  1:2000379/150‑1 (MQ=255)
                                                                       aaaaTGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGaa                                                                               >  2:1966199/1‑150 (MQ=255)
                                                                         aaTGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAAt                                                                             >  2:4451914/1‑150 (MQ=255)
                                                                          aTGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAAtt                                                                            <  1:3421777/150‑1 (MQ=255)
                                                                               tGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTgg                                                                       >  1:80340/1‑150 (MQ=255)
                                                                                        ttAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTa                                                              <  2:1314566/150‑1 (MQ=255)
                                                                                          aaaaaaGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAg                                                            <  2:2648305/150‑1 (MQ=255)
                                                                                           aaaaaGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAgg                                                           >  1:1276539/1‑150 (MQ=255)
                                                                                           aaaaaGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAgg                                                           >  1:1276481/1‑150 (MQ=255)
                                                                                           aaaaaGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAgg                                                           <  1:3809383/150‑1 (MQ=255)
                                                                                             aaaGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGgtg                                                         >  2:2453853/1‑150 (MQ=255)
                                                                                              aaGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGgtgt                                                        <  2:4393532/150‑1 (MQ=255)
                                                                                              aaGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGgtgt                                                        <  2:695355/150‑1 (MQ=255)
                                                                                              aaGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGgtgt                                                        >  2:3238784/1‑150 (MQ=255)
                                                                                                 tAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGtttt                                                     >  2:2944153/1‑150 (MQ=255)
                                                                                                  aGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTa                                                    >  2:3188505/1‑150 (MQ=255)
                                                                                                  aGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTa                                                    >  2:3138192/1‑150 (MQ=255)
                                                                                                   ggATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTAt                                                   <  2:32629/150‑1 (MQ=255)
                                                                                                   ggATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTAt                                                   >  2:3391358/1‑150 (MQ=255)
                                                                                                       ttATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTg                                               <  2:3111889/150‑1 (MQ=255)
                                                                                                       ttATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTg                                               <  2:3112072/150‑1 (MQ=255)
                                                                                                            gATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAATATCCTAATGTTATTTATCGTGAGATATTAAGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGAGGTGGTAAGGTGTTTTATGCTGTTAtt                                          >  2:4004036/1‑150 (MQ=255)
                                                                                                             aTAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTAttt                                         <  2:3354811/150‑1 (MQ=255)
                                                                                                               aaaGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTAttttt                                       <  1:1691882/150‑1 (MQ=255)
                                                                                                               aaaGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTAttttt                                       >  1:2983456/1‑150 (MQ=255)
                                                                                                               aaaGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTAttttt                                       >  1:975107/1‑150 (MQ=255)
                                                                                                                  gTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACTAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATg                                    >  2:4330790/1‑150 (MQ=255)
                                                                                                                   tAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATgc                                   >  2:1695511/1‑150 (MQ=255)
                                                                                                                   tAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATgc                                   >  1:165922/1‑150 (MQ=255)
                                                                                                                       cAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCAt                               >  2:1985039/1‑150 (MQ=255)
                                                                                                                        aaGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCAtt                              <  2:1549360/150‑1 (MQ=255)
                                                                                                                             cAAAGGACTCGTGTTTAAATAACTAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCtgtg                         <  2:1182008/150‑1 (MQ=255)
                                                                                                                             cAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGGAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCtgtg                         <  2:3489555/150‑1 (MQ=255)
                                                                                                                              aaaGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTgtgt                        <  2:2946656/150‑1 (MQ=255)
                                                                                                                              aaaGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTgtgt                        <  1:1380891/150‑1 (MQ=255)
                                                                                                                               aaGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGtc                       <  2:3852918/150‑1 (MQ=255)
                                                                                                                               aaGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGtc                       >  1:1856000/1‑150 (MQ=255)
                                                                                                                                 ggACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGtctc                     >  2:1390349/1‑150 (MQ=255)
                                                                                                                                 ggACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGtctc                     >  2:2426826/1‑150 (MQ=255)
                                                                                                                                   aCTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCt                   <  1:1921532/150‑1 (MQ=255)
                                                                                                                                     tCGTGTTTAAATAACAAGATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGa                 <  1:3300321/150‑1 (MQ=255)
                                                                                                                                           ttAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATc           >  2:949050/1‑150 (MQ=255)
                                                                                                                                            tAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCa          >  2:4099476/1‑150 (MQ=255)
                                                                                                                                             aaaTAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCAc         >  1:399257/1‑150 (MQ=255)
                                                                                                                                                   cAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTaaaa   <  1:4140496/150‑1 (MQ=255)
                                                                                                                                                    aaaaTCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTaaaaa  >  1:109016/1‑150 (MQ=255)
                                                                                                                                                    aaaaTCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCAAGTaaaaa  >  2:2665115/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAAA  >  CP009273/1566244‑1566541

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: