Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,099,378 A→C 14.0% intergenic (‑726/‑29) csgD ← / → csgB csgBAC operon transcriptional regulator/curlin nucleator protein, minor subunit in curli complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,099,3780AC14.0% 529.8 / 93.7 220intergenic (‑726/‑29)csgD/csgBcsgBAC operon transcriptional regulator/curlin nucleator protein, minor subunit in curli complex
Reads supporting (aligned to +/- strand):  ref base A (82/107);  new base C (20/11);  total (102/118)
Fisher's exact test for biased strand distribution p-value = 3.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 222 total aligned reads displayed.

TTAATACTTTGGTATGAACTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAG  >  CP009273/1099231‑1099526
                                                                                                                                                   |                                                                                                                                                    
    tACTTTGGTATGAACTAAAAAAGAAAAAGACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTa                                                                                                                                                <  2:808663/150‑1 (MQ=255)
     aCTTTGGTATGAACTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTAc                                                                                                                                               >  2:1607999/1‑150 (MQ=255)
     aCTTTGGTATGAACTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTAc                                                                                                                                               >  2:2526316/1‑150 (MQ=255)
     aCTTTGGTATGAACTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCACAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTAc                                                                                                                                               >  1:4338695/1‑150 (MQ=255)
       tttGGTATGAACTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGa                                                                                                                                             >  2:3645494/1‑150 (MQ=255)
          ggTATGAACTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCa                                                                                                                                          <  1:3004639/150‑1 (MQ=255)
            tATGAACTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTCAATGTACGACCAgg                                                                                                                                        >  2:1758639/1‑150 (MQ=255)
            tATGAACTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAgg                                                                                                                                        >  1:1974203/1‑150 (MQ=255)
            tATGAACTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAgg                                                                                                                                        >  1:2470398/1‑150 (MQ=255)
               gAACTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTcc                                                                                                                                     <  1:1013160/150‑1 (MQ=255)
               gAACTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTcc                                                                                                                                     >  1:410942/1‑150 (MQ=255)
                  cTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAATATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAgg                                                                                                                                  <  1:661845/150‑1 (MQ=255)
                  cTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAgg                                                                                                                                  <  1:1587113/150‑1 (MQ=255)
                  cTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAgg                                                                                                                                  <  2:2270091/150‑1 (MQ=255)
                     aaaaaGAAAAATACACCGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTg                                                                                                                               >  2:872050/1‑150 (MQ=255)
                     aaaaaGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTg                                                                                                                               >  2:4095610/1‑150 (MQ=255)
                       aaaGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGac                                                                                                                             <  2:2741412/150‑1 (MQ=255)
                       aaaGAAAAATACAACGCGCGGGGGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGac                                                                                                                             <  1:210781/150‑1 (MQ=255)
                        aaGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGaca                                                                                                                            >  2:1280004/1‑150 (MQ=255)
                           aaaaaTACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGacaaca                                                                                                                         <  1:858649/150‑1 (MQ=255)
                           aaaaaTACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGacaaca                                                                                                                         >  2:2880179/1‑150 (MQ=255)
                             aaaTACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATg                                                                                                                       <  2:3531478/150‑1 (MQ=255)
                              aaTACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGa                                                                                                                      >  1:1111635/1‑150 (MQ=255)
                               aTACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGaa                                                                                                                     <  2:1173345/150‑1 (MQ=255)
                                  cAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAACACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTCAATGTACGACCAGGTCCAGGGTGACAACATGaaaaa                                                                                                                  >  2:2926838/1‑150 (MQ=255)
                                   aaCGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAAc                                                                                                                 >  2:2092852/1‑150 (MQ=255)
                                   aaCGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAAc                                                                                                                 <  1:2047078/150‑1 (MQ=255)
                                      gcgcgGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACaaa                                                                                                              <  2:3236758/150‑1 (MQ=255)
                                        gcgGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAAtt                                                                                                            <  1:3568356/150‑1 (MQ=255)
                                        gcgGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAAtt                                                                                                            <  2:2336154/150‑1 (MQ=255)
                                         cgGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTg                                                                                                           <  2:3347609/150‑1 (MQ=255)
                                         cgGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTg                                                                                                           <  2:3347440/150‑1 (MQ=255)
                                          gggTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGt                                                                                                          >  1:3678669/1‑150 (MQ=255)
                                               aGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGttattt                                                                                                     >  2:3967150/1‑150 (MQ=255)
                                               aGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGttattt                                                                                                     <  1:2124252/150‑1 (MQ=255)
                                                   attaAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTCAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTatga                                                                                                 >  1:457370/1‑150 (MQ=255)
                                                   attaAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTatga                                                                                                 >  2:1027736/1‑150 (MQ=255)
                                                    ttaAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTatgat                                                                                                >  2:3327967/1‑150 (MQ=255)
                                                    ttaAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTatgat                                                                                                >  1:2294171/1‑150 (MQ=255)
                                                      aaaaaTATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGt                                                                                              >  1:4520444/1‑150 (MQ=255)
                                                       aaaaTATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTCAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGtt                                                                                             >  2:3613100/1‑150 (MQ=255)
                                                           tatTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTCAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACa                                                                                         <  2:328594/150‑1 (MQ=255)
                                                            atTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTCAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACaa                                                                                        >  1:3069468/1‑150 (MQ=255)
                                                               tCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATAc                                                                                     <  2:3103327/150‑1 (MQ=255)
                                                                ccGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACt                                                                                    >  1:22449/1‑150 (MQ=255)
                                                                ccGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACt                                                                                    <  1:2619827/150‑1 (MQ=255)
                                                                 cGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTg                                                                                   <  2:3388464/150‑1 (MQ=255)
                                                                   cAGACATACTTTCCATCGTAACGCAGCGTTAACAACATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTggg                                                                                 <  2:3458517/150‑1 (MQ=255)
                                                                     gACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTg                                                                               >  1:1247550/1‑150 (MQ=255)
                                                                      aCATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTgc                                                                              >  2:142899/1‑150 (MQ=255)
                                                                        aTACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTgcgc                                                                            >  2:1444094/1‑150 (MQ=255)
                                                                          aCTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCt                                                                          <  2:2023748/150‑1 (MQ=255)
                                                                             ttCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTCAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTggg                                                                       <  1:3056535/150‑1 (MQ=255)
                                                                             ttCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTggg                                                                       <  2:2896451/150‑1 (MQ=255)
                                                                               ccATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGAt                                                                     <  1:3214376/150‑1 (MQ=255)
                                                                                 aTCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTg                                                                   >  1:2373053/1‑150 (MQ=255)
                                                                                 aTCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTg                                                                   <  1:2711497/150‑1 (MQ=255)
                                                                                  tCGTAACGCAGCGTTAACAACATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGc                                                                  <  2:3437864/150‑1 (MQ=255)
                                                                                     tAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGc                                                               <  1:3418371/150‑1 (MQ=255)
                                                                                            gCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCgcagca                                                        <  2:3168348/150‑1 (MQ=255)
                                                                                                tAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGtt                                                    >  1:2898490/1‑150 (MQ=255)
                                                                                                 aaCAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTCAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTa                                                   >  1:906604/1‑150 (MQ=255)
                                                                                                 aaCAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTa                                                   >  1:1567244/1‑150 (MQ=255)
                                                                                                   cAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATg                                                 >  2:432547/1‑150 (MQ=255)
                                                                                                    aaaaTACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTCAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGa                                                >  2:3713628/1‑150 (MQ=255)
                                                                                                    aaaaTACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGa                                                >  1:4365279/1‑150 (MQ=255)
                                                                                                    aaaaTACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGa                                                >  1:2039266/1‑150 (MQ=255)
                                                                                                          cAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGc                                          >  2:4574204/1‑150 (MQ=255)
                                                                                                           tggTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCt                                         >  2:815731/2‑150 (MQ=255)
                                                                                                            ggTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTa                                        <  1:3622271/150‑1 (MQ=255)
                                                                                                            ggTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTa                                        <  2:1187605/150‑1 (MQ=255)
                                                                                                               tGCTTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAAtt                                     <  2:3197628/150‑1 (MQ=255)
                                                                                                               tGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAAtt                                     <  2:3674886/150‑1 (MQ=255)
                                                                                                                gCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTc                                    <  1:1465157/150‑1 (MQ=255)
                                                                                                                gCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTc                                    <  2:2766909/150‑1 (MQ=255)
                                                                                                                gCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTc                                    <  1:3498027/150‑1 (MQ=255)
                                                                                                                 cGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCa                                   <  1:2501236/150‑1 (MQ=255)
                                                                                                                  gTTAACAACCAAGTTGAAATGATTTAATTTCTTCAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAg                                  >  2:485155/1‑150 (MQ=255)
                                                                                                                   ttAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGa                                 >  2:2695554/1‑150 (MQ=255)
                                                                                                                    tAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGaa                                >  2:4228314/1‑150 (MQ=255)
                                                                                                                     aacaacCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAat                               <  2:943394/150‑1 (MQ=255)
                                                                                                                     aacaacCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAat                               <  1:667664/150‑1 (MQ=255)
                                                                                                                        aacCAAGTTGAAATGATTTAATTTCTTCAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAAtata                            >  1:1323304/1‑150 (MQ=255)
                                                                                                                        aacCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAAtata                            <  2:1327201/150‑1 (MQ=255)
                                                                                                                        aacCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAAtata                            <  2:4472282/150‑1 (MQ=255)
                                                                                                                          ccAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAAc                          <  2:2594503/150‑1 (MQ=255)
                                                                                                                            aaGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACtt                        >  2:2452647/1‑150 (MQ=255)
                                                                                                                            aaGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACtt                        <  1:394098/150‑1 (MQ=255)
                                                                                                                             aGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTc                       <  1:1049938/150‑1 (MQ=255)
                                                                                                                                  aaaTGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGt                  >  2:462597/1‑150 (MQ=255)
                                                                                                                                   aaTGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTa                 <  2:571985/150‑1 (MQ=255)
                                                                                                                                          tAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAAt          <  1:1911844/150‑1 (MQ=255)
                                                                                                                                           aaTTTCTTCAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATAATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAAtt         >  2:4358884/1‑150 (MQ=255)
                                                                                                                                             tttCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGa       >  1:670397/1‑150 (MQ=255)
                                                                                                                                             tttCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGa       >  2:3821616/1‑150 (MQ=255)
                                                                                                                                             tttCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGa       >  2:1920155/1‑150 (MQ=255)
                                                                                                                                              ttCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAg      >  1:2814281/1‑150 (MQ=255)
                                                                                                                                                cTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTa    >  2:1544648/1‑150 (MQ=255)
                                                                                                                                                 ttAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTaa   <  1:2603448/150‑1 (MQ=255)
                                                                                                                                                  tAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAg  <  1:801836/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTAATACTTTGGTATGAACTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAG  >  CP009273/1099231‑1099526

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: