Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,100,615 A→C 35.2% N67H (AAT→CAT)  csgC → curli assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,100,6150AC35.2% 138.5 / 275.3 231N67H (AAT→CAT) csgCcurli assembly protein
Reads supporting (aligned to +/- strand):  ref base A (74/76);  new base C (36/45);  total (110/121)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

Only 100 of 238 total aligned reads displayed.

ACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATAT  >  CP009273/1100468‑1100762
                                                                                                                                                   |                                                                                                                                                   
aCGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTa                                                                                                                                                     >  2:1878800/1‑148 (MQ=255)
aCGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTa                                                                                                                                                     <  1:1523703/148‑1 (MQ=255)
  gACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAAt                                                                                                                                                   >  1:2435052/1‑148 (MQ=255)
         cAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCa                                                                                                                                            >  1:167668/1‑148 (MQ=255)
         cAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCa                                                                                                                                            >  1:3026651/1‑148 (MQ=255)
         cAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCa                                                                                                                                            >  1:1248707/1‑148 (MQ=255)
         cAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCa                                                                                                                                            >  1:1927819/1‑148 (MQ=255)
           aGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCAtt                                                                                                                                          >  1:2053054/1‑148 (MQ=255)
                 tgtgTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTg                                                                                                                                    <  2:431660/148‑1 (MQ=255)
                      ataCCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAGGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGaa                                                                                                                               >  2:3563051/1‑148 (MQ=255)
                         ccATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGtt                                                                                                                            >  2:1379946/1‑148 (MQ=255)
                          cATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTg                                                                                                                           <  1:253982/148‑1 (MQ=255)
                            ttattCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAg                                                                                                                         <  1:825527/148‑1 (MQ=255)
                               ttCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGttt                                                                                                                      >  2:2474971/1‑148 (MQ=255)
                                tCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTa                                                                                                                     >  1:1302678/1‑148 (MQ=255)
                                 ccTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTaa                                                                                                                    <  1:477990/148‑1 (MQ=255)
                                 ccTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTaa                                                                                                                    <  2:1220686/148‑1 (MQ=255)
                                   tGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAat                                                                                                                  <  1:3353482/148‑1 (MQ=255)
                                   tGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAat                                                                                                                  <  2:3024610/148‑1 (MQ=255)
                                    gAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAata                                                                                                                 <  1:2706300/148‑1 (MQ=255)
                                      aGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATAtt                                                                                                               <  1:3267296/148‑1 (MQ=255)
                                         cACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTcc                                                                                                            <  1:1746727/148‑1 (MQ=255)
                                                aCTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACg                                                                                                     >  2:3291798/1‑148 (MQ=255)
                                                 cTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGa                                                                                                    <  1:3432425/148‑1 (MQ=255)
                                                      tCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCggg                                                                                               <  1:1878800/148‑1 (MQ=255)
                                                      tCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCggg                                                                                               >  2:779542/1‑148 (MQ=255)
                                                         tGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGa                                                                                            <  1:1125302/148‑1 (MQ=255)
                                                             tgtgCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAAt                                                                                        <  1:1570883/148‑1 (MQ=255)
                                                             tgtgCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAAt                                                                                        <  1:1092623/148‑1 (MQ=255)
                                                             tgtgCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAAt                                                                                        <  1:1092509/148‑1 (MQ=255)
                                                             tgtgCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAAt                                                                                        <  2:829751/148‑1 (MQ=255)
                                                              gtgCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATa                                                                                       >  1:470815/1‑148 (MQ=255)
                                                                gCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAgt                                                                                     <  2:3186689/148‑1 (MQ=255)
                                                                gCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAgt                                                                                     <  1:1434251/148‑1 (MQ=255)
                                                                   gagTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAgttgt                                                                                  <  1:227288/148‑1 (MQ=255)
                                                                   gagTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAgttgt                                                                                  >  1:1889555/1‑148 (MQ=255)
                                                                    agTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAgttgtt                                                                                 <  1:895477/148‑1 (MQ=255)
                                                                        cAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTg                                                                             <  1:2903996/148‑1 (MQ=255)
                                                                           atatTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGttt                                                                          <  1:592949/148‑1 (MQ=255)
                                                                            tatTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTc                                                                         <  1:730953/148‑1 (MQ=255)
                                                                               tGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGa                                                                      >  2:1735786/1‑148 (MQ=255)
                                                                               tGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGa                                                                      >  1:2824288/1‑148 (MQ=255)
                                                                                gTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGAt                                                                     <  1:2004925/148‑1 (MQ=255)
                                                                                    cTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGAc                                                                 >  1:1788778/1‑148 (MQ=255)
                                                                                    cTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGAc                                                                 <  1:3291798/148‑1 (MQ=255)
                                                                                    cTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGAc                                                                 <  1:608772/148‑1 (MQ=255)
                                                                                      gcgcGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAg                                                               >  1:1595025/1‑148 (MQ=255)
                                                                                       cgcgAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGt                                                              >  1:3449649/1‑148 (MQ=255)
                                                                                       cgcgAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGt                                                              <  2:2503877/148‑1 (MQ=255)
                                                                                        gcgAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTc                                                             <  1:22200/148‑1 (MQ=255)
                                                                                              gCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCa                                                       >  2:1740787/1‑148 (MQ=255)
                                                                                                 gTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCAttt                                                    <  2:676899/148‑1 (MQ=255)
                                                                                                   tCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTAt                                                  >  2:918564/1‑148 (MQ=255)
                                                                                                   tCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTAt                                                  <  2:743231/148‑1 (MQ=255)
                                                                                                   tCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTAt                                                  <  2:2464439/148‑1 (MQ=255)
                                                                                                      gggCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATcac                                               <  2:2789137/148‑1 (MQ=255)
                                                                                                       ggCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATcaca                                              <  1:798374/148‑1 (MQ=255)
                                                                                                           aaGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCAcaacaa                                          <  2:3260154/148‑1 (MQ=255)
                                                                                                             gTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATg                                        >  2:2853985/1‑148 (MQ=255)
                                                                                                               cAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGgc                                      >  1:2505721/1‑148 (MQ=255)
                                                                                                                  aCGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGgccgc                                   >  1:2600848/1‑148 (MQ=255)
                                                                                                                   cGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGgccgcc                                  <  2:2053054/148‑1 (MQ=255)
                                                                                                                     aaGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGCCGATCGGGTGAAAATAGTTGTTACGGTTTCTGATGGACAGTCACTTCATTTATCACCACAATGGCCGCCct                                >  2:602433/1‑148 (MQ=255)
                                                                                                                      aGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCctc                               >  2:2362525/1‑148 (MQ=255)
                                                                                                                      aGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCctc                               <  2:167668/148‑1 (MQ=255)
                                                                                                                      aGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCctc                               <  2:729461/148‑1 (MQ=255)
                                                                                                                      aGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCctc                               >  2:586529/1‑148 (MQ=255)
                                                                                                                      aGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACACTGGCCGCCctc                               >  2:3207656/1‑148 (MQ=255)
                                                                                                                       gCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCctct                              <  1:390359/148‑1 (MQ=255)
                                                                                                                       gCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCctct                              <  1:685817/148‑1 (MQ=255)
                                                                                                                        cAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCtt                             >  1:240831/1‑148 (MQ=255)
                                                                                                                        cAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCtt                             <  1:2473068/148‑1 (MQ=255)
                                                                                                                          aGAAAAGACCCTTTCAGTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCa                           <  1:2906752/148‑1 (MQ=255)
                                                                                                                           gAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAg                          >  2:2113410/1‑148 (MQ=255)
                                                                                                                           gAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAg                          <  1:992167/148‑1 (MQ=255)
                                                                                                                           gAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAg                          >  1:2675565/1‑148 (MQ=255)
                                                                                                                            aaaaGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGa                         >  1:2761605/1‑148 (MQ=255)
                                                                                                                            aaaaGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGa                         >  1:1012116/1‑148 (MQ=255)
                                                                                                                              aaGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGaaa                       <  1:3122005/148‑1 (MQ=255)
                                                                                                                              aaGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGaaa                       >  1:939564/1‑148 (MQ=255)
                                                                                                                               aGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGaaaa                      <  2:1369394/148‑1 (MQ=255)
                                                                                                                               aGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGaaaa                      <  2:2521802/148‑1 (MQ=255)
                                                                                                                                gACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAg                     >  2:1174876/1‑148 (MQ=255)
                                                                                                                                gACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAg                     >  2:114484/1‑148 (MQ=255)
                                                                                                                                gACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAg                     >  2:397478/1‑148 (MQ=255)
                                                                                                                                gACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAg                     >  2:3100193/1‑148 (MQ=255)
                                                                                                                                  cccTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTc                   <  2:3103574/148‑1 (MQ=255)
                                                                                                                                   ccTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTCCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCt                  <  2:3378697/148‑1 (MQ=255)
                                                                                                                                    cTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCtt                 >  1:2351935/1‑148 (MQ=255)
                                                                                                                                    cTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCtt                 >  1:2230405/1‑148 (MQ=255)
                                                                                                                                     tttCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTa                <  2:554830/148‑1 (MQ=255)
                                                                                                                                     tttCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTa                <  1:2200224/148‑1 (MQ=255)
                                                                                                                                        cATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAAtt             >  2:3201477/1‑148 (MQ=255)
                                                                                                                                          ttGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATttg           <  1:1390550/148‑1 (MQ=255)
                                                                                                                                          ttGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATttg           <  1:9277/148‑1 (MQ=255)
                                                                                                                                            gCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATttgtt         <  1:402706/148‑1 (MQ=255)
                                                                                                                                            gCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATttgtt         <  1:402713/148‑1 (MQ=255)
                                                                                                                                            gCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATttgtt         >  2:1044175/1‑148 (MQ=255)
                                                                                                                                                 cTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAat    >  1:3315011/1‑148 (MQ=255)
                                                                                                                                                   aaTCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAatat  >  1:3478415/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATAT  >  CP009273/1100468‑1100762

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: