Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,566,392 A→T 30.0% 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,3920AT30.0% 435.7 / 254.6 260intergenic (‑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 (34/44);  total (129/131)
Fisher's exact test for biased strand distribution p-value = 2.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 264 total aligned reads displayed.

GCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAA  >  CP009273/1566244‑1566540
                                                                                                                                                    |                                                                                                                                                    
 cACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACaaa                                                                                                                                                    <  2:4642316/150‑1 (MQ=255)
   ccAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACTAAAt                                                                                                                                                  <  2:2225969/150‑1 (MQ=255)
   ccAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAAt                                                                                                                                                  <  1:843631/150‑1 (MQ=255)
   ccAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAAt                                                                                                                                                  <  2:4745838/150‑1 (MQ=255)
    cAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATc                                                                                                                                                 <  1:3666856/150‑1 (MQ=255)
       aaCGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTa                                                                                                                                              <  2:4313110/150‑1 (MQ=255)
         cGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAAt                                                                                                                                            >  1:5055066/1‑150 (MQ=255)
           tGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGt                                                                                                                                          >  1:2329072/1‑150 (MQ=255)
            gAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGtt                                                                                                                                         >  2:52045/1‑150 (MQ=255)
              aTCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACTAAATCCTAATGttat                                                                                                                                       <  1:862355/150‑1 (MQ=255)
              aTCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGttat                                                                                                                                       >  1:3764732/1‑150 (MQ=255)
              aTCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGttat                                                                                                                                       <  2:4802198/150‑1 (MQ=255)
                cGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACTAAATCCTAATGttattt                                                                                                                                     <  1:3817063/150‑1 (MQ=255)
                  agtagtTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGttatttat                                                                                                                                   <  2:2064137/150‑1 (MQ=255)
                  agtagtTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGttatttat                                                                                                                                   <  2:4989228/150‑1 (MQ=255)
                   gtagtTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATc                                                                                                                                  >  1:3253108/1‑150 (MQ=255)
                      gtTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTg                                                                                                                               <  1:699187/150‑1 (MQ=255)
                       ttCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTga                                                                                                                              >  2:2356291/1‑150 (MQ=255)
                        tCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACTAAATCCTAATGTTATTTATCGTgag                                                                                                                             <  2:467432/150‑1 (MQ=255)
                        tCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTgag                                                                                                                             >  1:3345192/1‑150 (MQ=255)
                          cGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGat                                                                                                                           <  1:3074822/150‑1 (MQ=255)
                           gACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGata                                                                                                                          >  2:3138354/1‑150 (MQ=255)
                            aCCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGatat                                                                                                                         <  2:361136/150‑1 (MQ=255)
                             ccTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACTAAATCCTAATGTTATTTATCGTGAGATAtt                                                                                                                        >  1:4558060/1‑150 (MQ=255)
                             ccTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACTAAATCCTAATGTTATTTATCGTGAGATAtt                                                                                                                        >  2:2252251/1‑150 (MQ=255)
                             ccTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATAtt                                                                                                                        >  1:1352225/1‑150 (MQ=255)
                             ccTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATAtt                                                                                                                        >  2:2385756/1‑150 (MQ=255)
                             ccTTAAATCCGTTACTTGCTTCTTATCCAATTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATAtt                                                                                                                        <  2:3091887/150‑1 (MQ=255)
                              cTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTa                                                                                                                       >  1:1020554/1‑150 (MQ=255)
                                tAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTAcg                                                                                                                     <  1:4991874/150‑1 (MQ=255)
                                tAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTAcg                                                                                                                     >  1:4450409/1‑150 (MQ=255)
                                    tCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGaa                                                                                                                 >  1:486935/1‑150 (MQ=255)
                                     ccGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACTAAATCCTAATGTTATTTATCGTGAGATATTACGCGaat                                                                                                                >  1:2192639/1‑150 (MQ=255)
                                       gTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGaataa                                                                                                              >  1:2907959/1‑150 (MQ=255)
                                         tACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAataata                                                                                                            >  2:1707131/1‑150 (MQ=255)
                                           cTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACTAACTCCTAATGTTATTTATCGTGAGATATTACGCGAATAATAtt                                                                                                          <  1:3202802/150‑1 (MQ=255)
                                             tGCTTCTTATCCATTTTAAACTCCTTAAAATGATTTGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAATGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCTAATAATAtttt                                                                                                        <  1:1139134/150‑1 (MQ=255)
                                              gCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATAttttt                                                                                                       <  2:566381/150‑1 (MQ=255)
                                                 tnttATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCa                                                                                                    >  1:2000995/1‑150 (MQ=255)
                                                  cttATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCAt                                                                                                   >  1:1434447/1‑150 (MQ=255)
                                                  cttATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCAt                                                                                                   >  2:3054685/1‑150 (MQ=255)
                                                  cttATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCAt                                                                                                   >  2:3054959/1‑150 (MQ=255)
                                                  cttATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCAt                                                                                                   >  1:2420550/1‑150 (MQ=255)
                                                  cttATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCAt                                                                                                   <  2:1157247/150‑1 (MQ=255)
                                                    tATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACTAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTCATTGa                                                                                                <  1:1889273/150‑1 (MQ=255)
                                                    tATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTg                                                                                                 >  1:353584/1‑150 (MQ=255)
                                                      tCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGaa                                                                                               <  1:52045/150‑1 (MQ=255)
                                                        cATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACTAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTCATTGaaaaa                                                                                            >  2:3541796/1‑150 (MQ=255)
                                                         aTTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGaaaaa                                                                                            >  1:4111775/1‑150 (MQ=255)
                                                         aTTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGaaaaa                                                                                            >  2:379178/1‑150 (MQ=255)
                                                          ttttAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGAGAAAATGACACGTGTTTAAATAACAAAATCATAATGTTATTTATCGTGAGATATTACGCTAATAATATTTTTTCATTGAAAAac                                                                                           >  1:2611431/1‑150 (MQ=255)
                                                           tttAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACTAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAaca                                                                                          <  2:3075518/150‑1 (MQ=255)
                                                           tttAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACTAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAaca                                                                                          <  2:3262156/150‑1 (MQ=255)
                                                            ttAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAacaa                                                                                         >  1:3598939/1‑150 (MQ=255)
                                                             tAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAacaat                                                                                        <  2:1309714/150‑1 (MQ=255)
                                                               aaCTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACTAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAacaatac                                                                                      <  1:4148311/150‑1 (MQ=255)
                                                                    cTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAatat                                                                                 >  2:2856371/1‑150 (MQ=255)
                                                                      tAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACTAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGa                                                                               <  1:3294770/150‑1 (MQ=255)
                                                                       aaaaTGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGaa                                                                              <  2:1002026/150‑1 (MQ=255)
                                                                       aaaaTGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGaa                                                                              >  1:585462/1‑150 (MQ=255)
                                                                       aaaaTGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGaa                                                                              >  1:27166/1‑150 (MQ=255)
                                                                        aaaTGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACTAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGaaa                                                                             >  2:4612286/1‑150 (MQ=255)
                                                                        aaaTGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGaaa                                                                             <  2:1477843/150‑1 (MQ=255)
                                                                        aaaTGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGaaa                                                                             >  2:2952977/1‑150 (MQ=255)
                                                                           tGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACTAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTc                                                                          >  2:4077883/1‑150 (MQ=255)
                                                                             aTTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACTAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCtt                                                                        <  1:3747639/150‑1 (MQ=255)
                                                                              ttGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTg                                                                       >  2:2746623/1‑150 (MQ=255)
                                                                                ggATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTggg                                                                     <  2:193150/150‑1 (MQ=255)
                                                                                ggATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTggg                                                                     <  1:2133031/150‑1 (MQ=255)
                                                                                 gATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGggt                                                                    <  2:2187760/150‑1 (MQ=255)
                                                                                  aTCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACTAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGggtg                                                                   <  1:1501858/150‑1 (MQ=255)
                                                                                     gCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGggtggtg                                                                <  2:4407498/150‑1 (MQ=255)
                                                                                       aTTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACTAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGggtggtggt                                                              <  2:4573731/150‑1 (MQ=255)
                                                                                        ttAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTa                                                             >  1:4399319/1‑150 (MQ=255)
                                                                                          aaaaaaGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAg                                                           <  1:2261556/150‑1 (MQ=255)
                                                                                            aaaaGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGgt                                                         <  2:1432795/150‑1 (MQ=255)
                                                                                              aaGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGgtgt                                                       >  2:3950506/1‑150 (MQ=255)
                                                                                               aGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGtt                                                      >  2:3887052/1‑150 (MQ=255)
                                                                                                   ggATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTAt                                                  >  1:4010978/1‑150 (MQ=255)
                                                                                                    gATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACTAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATg                                                 <  2:1263487/150‑1 (MQ=255)
                                                                                                     aTTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACTAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCt                                               >  1:3938412/1‑150 (MQ=255)
                                                                                                        tATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGt                                             >  1:4806976/1‑150 (MQ=255)
                                                                                                              tAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACTAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTAtttt                                       >  1:3186167/1‑150 (MQ=255)
                                                                                                               aaaGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTAttttt                                      >  2:2277588/1‑150 (MQ=255)
                                                                                                                   tAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATgc                                  <  1:1859737/150‑1 (MQ=255)
                                                                                                                      gCAAGTGCAAAGGACTCGTGTTTAAATAACTAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCAt                              <  1:418805/150‑1 (MQ=255)
                                                                                                                      gCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCa                               <  1:1809750/150‑1 (MQ=255)
                                                                                                                       cAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCAt                              <  2:1191462/150‑1 (MQ=255)
                                                                                                                       cAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCAt                              <  1:1396912/150‑1 (MQ=255)
                                                                                                                         aGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATTc                            <  1:4268825/150‑1 (MQ=255)
                                                                                                                                aGGACTCGTGTTTAAATAACTAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGtctc                    >  1:4637282/1‑150 (MQ=255)
                                                                                                                                aGGACTCGTGTTTAAATAACTAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGtctc                    >  1:4637216/1‑150 (MQ=255)
                                                                                                                                   aCTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCt                  <  1:3073002/150‑1 (MQ=255)
                                                                                                                                    cTCGTGTTTAAATAACTAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTg                 >  1:808513/1‑150 (MQ=255)
                                                                                                                                      cGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGaa               >  2:766578/1‑150 (MQ=255)
                                                                                                                                       gtgtTTAAATAACTAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAAtt             >  2:4578571/1‑150 (MQ=255)
                                                                                                                                            tAAATAACTAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCa         <  2:1818964/150‑1 (MQ=255)
                                                                                                                                            tAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCa         >  2:1415279/1‑150 (MQ=255)
                                                                                                                                             aaaTAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCAc        >  2:3420679/1‑150 (MQ=255)
                                                                                                                                              aaTAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACg       >  1:1666473/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATTCTGTGTCTCCTGAATTATCACGTAAAA  >  CP009273/1566244‑1566540

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: