Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,100,615 A→C 13.6% 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,6150AC13.6% 715.8 / 125.6 314N67H (AAT→CAT) csgCcurli assembly protein
Reads supporting (aligned to +/- strand):  ref base A (132/139);  new base C (23/20);  total (155/159)
Fisher's exact test for biased strand distribution p-value = 6.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

Only 100 of 319 total aligned reads displayed.

ACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATAT  >  CP009273/1100468‑1100762
                                                                                                                                                   |                                                                                                                                                   
aCGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTa                                                                                                                                                     >  2:1853717/1‑148 (MQ=255)
aCGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTa                                                                                                                                                     >  1:22102/1‑148 (MQ=255)
aCGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTa                                                                                                                                                     <  2:3019110/148‑1 (MQ=255)
      cAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAAc                                                                                                                                               >  1:2504405/1‑148 (MQ=255)
          aaGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCAt                                                                                                                                           <  1:3750002/148‑1 (MQ=255)
             gggATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGc                                                                                                                                        >  2:2270504/1‑148 (MQ=255)
                    gtATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACg                                                                                                                                 >  1:862562/1‑148 (MQ=255)
                     tataCCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGa                                                                                                                                <  1:2481462/148‑1 (MQ=255)
                        ggcgcaattCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGt                                                                                                                             <  2:671698/142‑1 (MQ=255)
                          cATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTg                                                                                                                           >  1:1058121/1‑148 (MQ=255)
                          cATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTg                                                                                                                           >  1:2272400/1‑148 (MQ=255)
                           attattCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGa                                                                                                                          <  2:1454420/148‑1 (MQ=255)
                             tattCCTGAAGTCGCTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGt                                                                                                                        >  1:1943308/1‑148 (MQ=255)
                             tattCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGt                                                                                                                        <  1:467489/148‑1 (MQ=255)
                             tattCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGt                                                                                                                        >  2:1261056/1‑148 (MQ=255)
                              attCCTGAAGTCTCTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGtt                                                                                                                       <  2:598050/148‑1 (MQ=255)
                              attCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGtt                                                                                                                       <  2:1670772/148‑1 (MQ=255)
                                   tGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAat                                                                                                                  <  2:183015/148‑1 (MQ=255)
                                   tGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAat                                                                                                                  <  1:958493/148‑1 (MQ=255)
                                   tGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAat                                                                                                                  >  1:2425945/1‑148 (MQ=255)
                                     aaGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAatat                                                                                                                >  2:1174097/1‑148 (MQ=255)
                                     aaGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATAGCTTTGACGAAGTTGAGTTTAAatat                                                                                                                >  2:3079907/1‑148 (MQ=255)
                                      aGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATAtt                                                                                                               >  2:1021712/1‑148 (MQ=255)
                                       gTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATAttt                                                                                                              >  2:834322/1‑148 (MQ=255)
                                        tCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTc                                                                                                             >  1:1415727/1‑148 (MQ=255)
                                           ctctTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTcccc                                                                                                          >  2:3065703/1‑148 (MQ=255)
                                            tctTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCg                                                                                                         >  1:2005090/1‑148 (MQ=255)
                                             ctTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCCAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCgg                                                                                                        >  2:3233213/1‑148 (MQ=255)
                                             ctTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCgg                                                                                                        >  2:1929965/1‑148 (MQ=255)
                                                aCTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACg                                                                                                     <  1:144385/148‑1 (MQ=255)
                                                  tCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGAt                                                                                                   >  1:2354823/1‑148 (MQ=255)
                                                     aTCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCgg                                                                                                >  2:919669/1‑148 (MQ=255)
                                                       cTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGt                                                                                              >  1:15617/1‑148 (MQ=255)
                                                       cTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGt                                                                                              >  2:3745689/1‑148 (MQ=255)
                                                            cTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGaaaa                                                                                         >  1:1189037/1‑148 (MQ=255)
                                                             tgtgCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAAt                                                                                        <  1:1510136/148‑1 (MQ=255)
                                                              gtgCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATa                                                                                       <  1:1021548/148‑1 (MQ=255)
                                                               tgCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAg                                                                                      <  1:702610/148‑1 (MQ=255)
                                                               tgCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAg                                                                                      <  1:1260524/148‑1 (MQ=255)
                                                                  agagTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAgttg                                                                                   <  2:1804006/148‑1 (MQ=255)
                                                                  agagTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAgttg                                                                                   >  2:1568819/1‑148 (MQ=255)
                                                                    agTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAgttgtt                                                                                 >  2:2165868/1‑148 (MQ=255)
                                                                         aaaTATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGt                                                                            >  1:578541/1‑148 (MQ=255)
                                                                         aaaTATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGt                                                                            <  2:2883374/148‑1 (MQ=255)
                                                                            tatTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTc                                                                         >  2:2697614/1‑148 (MQ=255)
                                                                               tGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGa                                                                      <  2:1943308/148‑1 (MQ=255)
                                                                               tGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGa                                                                      <  2:949434/148‑1 (MQ=255)
                                                                                gTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGAt                                                                     <  2:1774078/148‑1 (MQ=255)
                                                                                 tCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATg                                                                    >  1:1301652/1‑148 (MQ=255)
                                                                                 tCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATg                                                                    >  1:3251939/1‑148 (MQ=255)
                                                                                  cccTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATgg                                                                   >  1:3011693/1‑148 (MQ=255)
                                                                                  cccTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATgg                                                                   >  2:2418299/1‑148 (MQ=255)
                                                                                   ccTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGa                                                                  >  1:3363432/1‑148 (MQ=255)
                                                                                    cTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGAc                                                                 >  2:3400445/1‑148 (MQ=255)
                                                                                    cTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGAc                                                                 >  2:1356637/1‑148 (MQ=255)
                                                                                      gcgcGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAg                                                               <  2:1816500/148‑1 (MQ=255)
                                                                                      gcgcGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAg                                                               <  2:776223/148‑1 (MQ=255)
                                                                                       cgcgAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGt                                                              <  1:2080323/148‑1 (MQ=255)
                                                                                          gAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCAc                                                           >  1:1381139/1‑148 (MQ=255)
                                                                                               cAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCAt                                                      <  2:1419499/148‑1 (MQ=255)
                                                                                               cAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCAt                                                      >  1:3280746/1‑148 (MQ=255)
                                                                                                 gTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCAttt                                                    <  2:2869382/148‑1 (MQ=255)
                                                                                                  ttCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTa                                                   >  1:737307/1‑148 (MQ=255)
                                                                                                     aGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTCTCTGATGGACAGTCACTTCATTTATca                                                <  1:1838442/148‑1 (MQ=255)
                                                                                                       ggCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATcaca                                              >  1:490369/1‑148 (MQ=255)
                                                                                                       ggCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATcaca                                              <  1:3136653/148‑1 (MQ=255)
                                                                                                         cAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCacaac                                            >  1:1673845/1‑148 (MQ=255)
                                                                                                            aGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAAt                                         <  2:2228967/148‑1 (MQ=255)
                                                                                                             gTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATg                                        >  1:642537/1‑148 (MQ=255)
                                                                                                             gTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATg                                        >  1:1333771/1‑148 (MQ=255)
                                                                                                             gTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATg                                        <  2:3430428/148‑1 (MQ=255)
                                                                                                                aGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGgcc                                     >  2:2967256/1‑148 (MQ=255)
                                                                                                                  aCGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGgccgc                                   >  1:3514138/1‑148 (MQ=255)
                                                                                                                   cGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGgccgcc                                  <  2:2739553/148‑1 (MQ=255)
                                                                                                                     aaGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCct                                >  1:1781955/1‑148 (MQ=255)
                                                                                                                       gCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCctct                              >  1:3773918/1‑148 (MQ=255)
                                                                                                                        cAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCtt                             >  1:767948/1‑148 (MQ=255)
                                                                                                                            aaaaGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGa                         >  1:2943898/1‑148 (MQ=255)
                                                                                                                              aaGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGaaa                       <  1:817963/148‑1 (MQ=255)
                                                                                                                              aaGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGaaa                       <  1:2900461/148‑1 (MQ=255)
                                                                                                                               aGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAATTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGaaaa                      <  1:2831945/148‑1 (MQ=255)
                                                                                                                               aGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGaaaa                      >  1:2115395/1‑148 (MQ=255)
                                                                                                                               aGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGaaaa                      >  1:2602995/1‑148 (MQ=255)
                                                                                                                               aGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGaaaa                      <  1:1238460/148‑1 (MQ=255)
                                                                                                                               aGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGaaaa                      <  2:473578/148‑1 (MQ=255)
                                                                                                                               aGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGaaaa                      <  2:505529/148‑1 (MQ=255)
                                                                                                                                 aCCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGt                    >  2:3655064/1‑148 (MQ=255)
                                                                                                                                  cccTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCATAAAAGTc                   >  1:2742779/1‑148 (MQ=255)
                                                                                                                                   ccTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCt                  <  2:2300916/148‑1 (MQ=255)
                                                                                                                                   ccTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCt                  >  2:3497846/1‑148 (MQ=255)
                                                                                                                                      ttCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTaa               >  2:3121869/1‑148 (MQ=255)
                                                                                                                                      ttCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTaa               <  2:3553674/148‑1 (MQ=255)
                                                                                                                                       tCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAAt              <  2:2370925/148‑1 (MQ=255)
                                                                                                                                           tGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATttgt          <  1:1135790/148‑1 (MQ=255)
                                                                                                                                           tGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATttgt          <  2:3566669/148‑1 (MQ=255)
                                                                                                                                            gCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATttgtt         >  1:889658/1‑148 (MQ=255)
                                                                                                                                               ttcTAATCAACCCAGTGCTTTGAAGAAGTTGAGTTTAAATATTTCCCCGGAAGATCGGGTGAAATTAGTTTTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGaa      <  1:480513/146‑1 (MQ=255)
                                                                                                                                                gCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGaaa     <  1:2923275/148‑1 (MQ=255)
                                                                                                                                                gCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGaaa     >  2:2577498/1‑148 (MQ=255)
                                                                                                                                                   aaTCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAatat  >  2:1300195/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: