Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,099,378 A→C 10.6% 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,3780AC10.6% 726.3 / 78.0 274intergenic (‑726/‑29)csgD/csgBcsgBAC operon transcriptional regulator/curlin nucleator protein, minor subunit in curli complex
Reads supporting (aligned to +/- strand):  ref base A (115/130);  new base C (11/18);  total (126/148)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 277 total aligned reads displayed.

TTTAATACTTTGGTATGAACTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAG  >  CP009273/1099230‑1099526
                                                                                                                                                    |                                                                                                                                                    
   aaTACTTTGGTATGAACTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAAt                                                                                                                                                   >  2:4321331/1‑149 (MQ=255)
    aTACTTTGGTATGAACTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATg                                                                                                                                                  >  2:4158109/1‑149 (MQ=255)
    aTACTTTGGTATGAACTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATg                                                                                                                                                  >  2:4236594/1‑149 (MQ=255)
     tACTTTGGTATGAACTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGt                                                                                                                                                 >  1:2029942/1‑149 (MQ=255)
      aCTTTGGTATGAACTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTa                                                                                                                                                >  1:2325954/1‑149 (MQ=255)
        tttGGTATGAACTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACg                                                                                                                                              >  1:4250494/1‑149 (MQ=255)
          tGGTATGAACTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGAc                                                                                                                                            <  1:881218/149‑1 (MQ=255)
            gTATGAACTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCa                                                                                                                                          <  2:2486314/149‑1 (MQ=255)
               tGAACTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGt                                                                                                                                       >  2:4451840/1‑149 (MQ=255)
               tGAACTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGt                                                                                                                                       <  1:1900299/149‑1 (MQ=255)
                   cTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTCAATGTACGACCAGGTCCAg                                                                                                                                   <  2:4529800/149‑1 (MQ=255)
                   cTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAg                                                                                                                                   <  1:4233819/149‑1 (MQ=255)
                    tAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAgg                                                                                                                                  >  1:498610/1‑149 (MQ=255)
                      aaaaaGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGt                                                                                                                                >  2:3557457/1‑149 (MQ=255)
                       aaaaGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTCAATGTACGACCAGGTCCAGGGTg                                                                                                                               >  1:831527/1‑149 (MQ=255)
                       aaaaGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTg                                                                                                                               >  1:4608429/1‑149 (MQ=255)
                       aaaaGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTg                                                                                                                               >  1:294909/1‑149 (MQ=255)
                       aaaaGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTg                                                                                                                               >  1:4383054/1‑149 (MQ=255)
                        aaaGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGa                                                                                                                              >  1:911840/1‑149 (MQ=255)
                        aaaGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGa                                                                                                                              <  1:2310236/149‑1 (MQ=255)
                         aaGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGac                                                                                                                             <  2:3510301/149‑1 (MQ=255)
                                  aCAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGaaa                                                                                                                    <  1:3458201/149‑1 (MQ=255)
                                   cAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTCAATGTACGACCAGGTCCAGGGTGACAACATGaaaa                                                                                                                   >  2:1267292/1‑149 (MQ=255)
                                    aaCGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGaaaaa                                                                                                                  <  1:4520163/149‑1 (MQ=255)
                                     aCGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAAc                                                                                                                 >  2:1020879/1‑149 (MQ=255)
                                         gcgGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAAt                                                                                                             <  1:357463/149‑1 (MQ=255)
                                            ggTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGt                                                                                                          <  1:1385349/149‑1 (MQ=255)
                                            ggTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGt                                                                                                          >  1:3901499/1‑149 (MQ=255)
                                               gAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGttat                                                                                                       <  2:1180842/149‑1 (MQ=255)
                                               gAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGttat                                                                                                       >  2:2519051/1‑149 (MQ=255)
                                                 gTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGttattt                                                                                                     <  2:383520/149‑1 (MQ=255)
                                                 gTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGttattt                                                                                                     >  2:2423211/1‑149 (MQ=255)
                                                 gTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGttattt                                                                                                     <  1:221784/149‑1 (MQ=255)
                                                  ttattaAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGGTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGttattta                                                                                                    <  2:2882551/149‑1 (MQ=255)
                                                       aaaaaTATTTCCGCAGACATACTTTCCATCGTAACGCAGCTTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTatgatg                                                                                               <  2:2521616/149‑1 (MQ=255)
                                                       aaaaaTATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTatgatg                                                                                               >  2:3753121/1‑149 (MQ=255)
                                                            tatTTCCGCAGACATACTGTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAAc                                                                                          <  1:2998616/149‑1 (MQ=255)
                                                              tttCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACaa                                                                                        <  2:3060849/149‑1 (MQ=255)
                                                               ttCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAAt                                                                                       <  1:3985753/149‑1 (MQ=255)
                                                               ttCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAAt                                                                                       <  1:3252064/149‑1 (MQ=255)
                                                                     aGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTggg                                                                                 >  2:2698471/1‑149 (MQ=255)
                                                                      taaaTAATTTCCATCTTAACGAAGAGTTAACAAAATACAGGTTGCGTTAACAACAAATTTGAAATGATTTAATTTATTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGt                                                                                <  1:2481884/146‑1 (MQ=255)
                                                                        cATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTgc                                                                              <  1:486787/149‑1 (MQ=255)
                                                                               tCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTggg                                                                       >  2:2635071/1‑149 (MQ=255)
                                                                                 cATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTCAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGAt                                                                     <  1:1136967/149‑1 (MQ=255)
                                                                                  aTCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGAtt                                                                    >  2:315150/1‑149 (MQ=255)
                                                                                  aTCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCATGGGAtt                                                                    >  2:1383844/1‑149 (MQ=255)
                                                                                     gTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCa                                                                 >  1:674157/1‑149 (MQ=255)
                                                                                      tAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAg                                                                >  2:3542301/1‑149 (MQ=255)
                                                                                        aCGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGcc                                                              >  2:1110292/1‑149 (MQ=255)
                                                                                        aCGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGcc                                                              >  2:1110498/1‑149 (MQ=255)
                                                                                             gctttAACAAAATACAGGTTGCTTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGAACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCgcagc                                                         <  1:4059296/146‑1 (MQ=255)
                                                                                             gCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCgcagc                                                         >  2:4059224/1‑149 (MQ=255)
                                                                                             gCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCgcagc                                                         >  1:4657694/1‑149 (MQ=255)
                                                                                               gTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGcagcag                                                       >  2:4633183/1‑149 (MQ=255)
                                                                                                 tAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGt                                                     >  2:526846/1‑149 (MQ=255)
                                                                                                  aaCAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGtt                                                    <  1:3620475/149‑1 (MQ=255)
                                                                                                   aCAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTa                                                   <  1:1394957/149‑1 (MQ=255)
                                                                                                    cAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTAt                                                  <  1:4562105/149‑1 (MQ=255)
                                                                                                      aaaTACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGa                                                <  2:3690/149‑1 (MQ=255)
                                                                                                      aaaTACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGa                                                >  2:1789234/1‑149 (MQ=255)
                                                                                                      aaaTACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGa                                                >  2:4696715/1‑149 (MQ=255)
                                                                                                      aaaTACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGa                                                <  1:1588051/149‑1 (MQ=255)
                                                                                                         tACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGGTATGAttt                                             >  1:2877932/1‑149 (MQ=255)
                                                                                                             ggTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCt                                         <  2:1621033/149‑1 (MQ=255)
                                                                                                                tGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAAt                                      <  2:850036/149‑1 (MQ=255)
                                                                                                                tGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAAt                                      <  2:806104/149‑1 (MQ=255)
                                                                                                                 gCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAAtt                                     <  1:2166560/149‑1 (MQ=255)
                                                                                                                   gTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCa                                   <  1:2712270/149‑1 (MQ=255)
                                                                                                                   gTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCa                                   >  2:3860172/1‑149 (MQ=255)
                                                                                                                   gTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCa                                   >  2:1154682/1‑149 (MQ=255)
                                                                                                                   gTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCa                                   <  2:14012/149‑1 (MQ=255)
                                                                                                                   gTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCa                                   <  1:1852454/149‑1 (MQ=255)
                                                                                                                    ttAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAg                                  >  1:107041/1‑149 (MQ=255)
                                                                                                                        caacCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAata                              <  2:1942009/149‑1 (MQ=255)
                                                                                                                        caacCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTAAACAATACTGGGTGCGCCTGGAATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAata                              <  1:4421395/149‑1 (MQ=255)
                                                                                                                          acCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAAtata                            <  2:3595439/149‑1 (MQ=255)
                                                                                                                              aGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACtt                        <  1:2145923/149‑1 (MQ=255)
                                                                                                                               gTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTc                       <  2:127678/149‑1 (MQ=255)
                                                                                                                               gTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTc                       <  1:3608554/149‑1 (MQ=255)
                                                                                                                                ttGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTcg                      <  1:4376372/149‑1 (MQ=255)
                                                                                                                                 tGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTcgc                     <  1:380852/149‑1 (MQ=255)
                                                                                                                                  gAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTcgcg                    >  2:3661751/1‑149 (MQ=255)
                                                                                                                                   aaaTGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCgg                   >  1:1739472/1‑149 (MQ=255)
                                                                                                                                   aaaTGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCgg                   >  2:3031184/1‑149 (MQ=255)
                                                                                                                                    aaTGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGt                  >  1:3743897/1‑149 (MQ=255)
                                                                                                                                     aTGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTa                 >  1:2216174/1‑149 (MQ=255)
                                                                                                                                     aTGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTa                 <  2:3442582/149‑1 (MQ=255)
                                                                                                                                      tGAGTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACAGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTaa                >  1:4106345/1‑149 (MQ=255)
                                                                                                                                        atatAATTTCTTCAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAAt              <  2:4038697/146‑1 (MQ=255)
                                                                                                                                         tttAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATg             <  1:1242724/149‑1 (MQ=255)
                                                                                                                                         tttAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATg             >  1:3122188/1‑149 (MQ=255)
                                                                                                                                            aaTTTCTTCAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAAt          >  1:1063554/1‑149 (MQ=255)
                                                                                                                                            aaTTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAAt          <  2:4349202/149‑1 (MQ=255)
                                                                                                                                            aaTTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAAt          >  2:3221857/1‑149 (MQ=255)
                                                                                                                                            aaTTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAAt          <  2:1576913/149‑1 (MQ=255)
                                                                                                                                              tttCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTg        <  2:4395293/149‑1 (MQ=255)
                                                                                                                                              tttCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTg        >  2:2319689/1‑149 (MQ=255)
                                                                                                                                                tCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAg      <  2:2211544/149‑1 (MQ=255)
                                                                                                                                                    aaaTGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAg  >  2:325112/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTTAATACTTTGGTATGAACTAAAAAAGAAAAATACAACGCGCGGGTGAGTTATTAAAAATATTTCCGCAGACATACTTTCCATCGTAACGCAGCGTTAACAAAATACAGGTTGCGTTAACAACCAAGTTGAAATGATTTAATTTCTTAAATGTACGACCAGGTCCAGGGTGACAACATGAAAAACAAATTGTTATTTATGATGTTAACAATACTGGGTGCGCCTGGGATTGCAGCCGCAGCAGGTTATGATTTAGCTAATTCAGAATATAACTTCGCGGTAAATGAATTGAGTAAG  >  CP009273/1099230‑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: