Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,100,615 A→C 40.9% 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,6150AC40.9% 48.6 / inf 225N67H (AAT→CAT) csgCcurli assembly protein
Reads supporting (aligned to +/- strand):  ref base A (55/78);  new base C (45/47);  total (100/125)
Fisher's exact test for biased strand distribution p-value = 2.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 231 total aligned reads displayed.

TACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATC  >  CP009273/1100467‑1100763
                                                                                                                                                    |                                                                                                                                                    
tACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTa                                                                                                                                                      >  2:370677/1‑149 (MQ=255)
tACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTa                                                                                                                                                      >  2:2674329/1‑149 (MQ=255)
   gACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATc                                                                                                                                                   >  1:1898228/1‑149 (MQ=255)
     cccAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCaa                                                                                                                                                 <  2:7081/149‑1 (MQ=255)
         gCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCa                                                                                                                                             >  2:1792972/1‑149 (MQ=255)
         gCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCa                                                                                                                                             <  1:3717541/149‑1 (MQ=255)
           aaGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCAtt                                                                                                                                           >  1:1357003/1‑149 (MQ=255)
              gggATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCt                                                                                                                                        >  1:1079943/1‑149 (MQ=255)
               ggATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCtt                                                                                                                                       >  1:631183/1‑149 (MQ=255)
                gATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACACATTGCttt                                                                                                                                      >  2:3969487/1‑149 (MQ=255)
                gATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCttt                                                                                                                                      >  1:3015598/1‑149 (MQ=255)
                 aTGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTg                                                                                                                                     <  1:1925573/149‑1 (MQ=255)
                  tgtgTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGa                                                                                                                                    >  1:1988751/1‑149 (MQ=255)
                    tgtATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACg                                                                                                                                  <  2:4061356/149‑1 (MQ=255)
                        taCCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGt                                                                                                                              <  2:885222/149‑1 (MQ=255)
                           cATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGa                                                                                                                           <  1:1965079/149‑1 (MQ=255)
                            attattCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAg                                                                                                                          >  1:385677/1‑149 (MQ=255)
                            attattCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAg                                                                                                                          >  2:1282235/1‑149 (MQ=255)
                            attattCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAg                                                                                                                          >  2:3303891/1‑149 (MQ=255)
                              tattCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGtt                                                                                                                        >  2:161579/1‑149 (MQ=255)
                                 tCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTaa                                                                                                                     <  2:1334287/149‑1 (MQ=255)
                                 tCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTaa                                                                                                                     <  2:3384154/149‑1 (MQ=255)
                                    tGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAata                                                                                                                  >  2:2117641/1‑149 (MQ=255)
                                     gAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAatat                                                                                                                 >  2:3091244/1‑149 (MQ=255)
                                      aaGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCCAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATAtt                                                                                                                >  2:2370668/1‑149 (MQ=255)
                                        gTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTc                                                                                                              <  2:1349928/149‑1 (MQ=255)
                                         tCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTcc                                                                                                             <  2:4014546/149‑1 (MQ=255)
                                         tCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTcc                                                                                                             <  1:1591589/149‑1 (MQ=255)
                                             tctTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCgg                                                                                                         >  2:2025958/1‑149 (MQ=255)
                                               ttACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGAc                                                                                                       >  2:2229150/1‑149 (MQ=255)
                                                  cTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGAt                                                                                                    <  1:2456172/149‑1 (MQ=255)
                                                  cTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGAt                                                                                                    >  2:605443/1‑149 (MQ=255)
                                                   tCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATc                                                                                                   >  1:3194914/1‑149 (MQ=255)
                                                     aaTCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCgg                                                                                                 >  1:3878009/1‑149 (MQ=255)
                                                     aaTCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCgg                                                                                                 >  1:1840403/1‑149 (MQ=255)
                                                      aTCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCggg                                                                                                <  2:303242/149‑1 (MQ=255)
                                                         ttGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGa                                                                                             >  2:2894821/1‑149 (MQ=255)
                                                          tGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGaa                                                                                            <  1:1823753/149‑1 (MQ=255)
                                                            tCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGaaaa                                                                                          <  1:192273/149‑1 (MQ=255)
                                                                tgCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAgt                                                                                      <  2:1389724/149‑1 (MQ=255)
                                                                     agTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTa                                                                                 >  2:3442288/1‑149 (MQ=255)
                                                                      gTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTAc                                                                                <  1:1296048/149‑1 (MQ=255)
                                                                      gTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTACAGTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTAc                                                                                <  2:2666397/149‑1 (MQ=255)
                                                                       tACAACTATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACt                                                                               <  2:3897982/149‑1 (MQ=255)
                                                                          aaaTATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGtt                                                                            <  1:3868823/149‑1 (MQ=255)
                                                                           aaTATTGTCCCTGCGCTAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGttt                                                                           <  1:3791305/149‑1 (MQ=255)
                                                                            atatTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTc                                                                          <  1:2278051/149‑1 (MQ=255)
                                                                               ttGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGa                                                                       >  1:1749142/1‑149 (MQ=255)
                                                                                 gTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATg                                                                     >  1:1474686/1‑149 (MQ=255)
                                                                                  tCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATgg                                                                    <  1:810176/149‑1 (MQ=255)
                                                                                   cccTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGa                                                                   <  2:1003266/149‑1 (MQ=255)
                                                                                   cccTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGa                                                                   <  1:2602234/149‑1 (MQ=255)
                                                                                    ccTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATGTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGAc                                                                  <  1:3779094/149‑1 (MQ=255)
                                                                                          cgAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCAc                                                            <  1:3838755/149‑1 (MQ=255)
                                                                                           gAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACt                                                           >  2:1961497/1‑149 (MQ=255)
                                                                                            aaGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACtt                                                          <  1:3656125/149‑1 (MQ=255)
                                                                                             aGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTc                                                         >  1:2838299/1‑149 (MQ=255)
                                                                                              ggCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCa                                                        <  1:896085/149‑1 (MQ=255)
                                                                                                cAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCAtt                                                      <  1:3546301/149‑1 (MQ=255)
                                                                                                 aGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCAttt                                                     <  1:3395474/149‑1 (MQ=255)
                                                                                                 aGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCAttt                                                     <  2:462256/149‑1 (MQ=255)
                                                                                                    tCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATc                                                  <  2:2388810/149‑1 (MQ=255)
                                                                                                    tCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATc                                                  <  1:306876/149‑1 (MQ=255)
                                                                                                     cAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCAGTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATca                                                 <  1:790780/149‑1 (MQ=255)
                                                                                                     cAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATca                                                 <  2:1631760/149‑1 (MQ=255)
                                                                                                      aGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATcac                                                <  2:2817606/149‑1 (MQ=255)
                                                                                                       gggCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATcaca                                               >  1:2146743/1‑149 (MQ=255)
                                                                                                         gCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCacaac                                             <  1:274157/149‑1 (MQ=255)
                                                                                                         gCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCacaac                                             >  1:3445825/1‑149 (MQ=255)
                                                                                                         gCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCCGATGGACAGTCACTTCATTTATCacaac                                             <  1:719137/149‑1 (MQ=255)
                                                                                                           aaaGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCAcaacaa                                           <  2:3724415/149‑1 (MQ=255)
                                                                                                           aaaGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCAcaacaa                                           <  1:1484620/149‑1 (MQ=255)
                                                                                                              gTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATgg                                        <  1:1750669/149‑1 (MQ=255)
                                                                                                               tCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGgc                                       >  2:3166091/1‑149 (MQ=255)
                                                                                                                  gaagAAGCTAGAAAAGACCCTTTCATTGCGGGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCCGGTGAAAATAGTTGTAACTGTATCTGATGGACAGTCACTTAATTTATCACAACAATGgacga                                    >  2:4066735/4‑148 (MQ=255)
                                                                                                                        gCAATAAAAGACCCTTTCATTGCCTGCTCATCAACCCTTTGCTTTGAAGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCtt                              <  1:3969487/149‑1 (MQ=255)
                                                                                                                        gCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCtt                              <  1:561217/149‑1 (MQ=255)
                                                                                                                              aaaGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGaaa                        >  1:151834/1‑149 (MQ=255)
                                                                                                                               aaGACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGaaaa                       <  1:1955385/149‑1 (MQ=255)
                                                                                                                               aaGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGaaaa                       >  2:2265020/1‑149 (MQ=255)
                                                                                                                               aaGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGaaaa                       >  2:3839118/1‑149 (MQ=255)
                                                                                                                                aGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAg                      >  1:3925869/1‑149 (MQ=255)
                                                                                                                                 gACCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGt                     <  2:2518143/149‑1 (MQ=255)
                                                                                                                                  aCCCTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTc                    >  1:584032/1‑149 (MQ=255)
                                                                                                                                  aCCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTc                    <  1:370677/149‑1 (MQ=255)
                                                                                                                                  aCCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTc                    <  1:2674329/149‑1 (MQ=255)
                                                                                                                                  aCCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTc                    <  1:2476206/149‑1 (MQ=255)
                                                                                                                                   cccTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCt                   >  2:1914558/1‑149 (MQ=255)
                                                                                                                                    ccTTTCATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCtt                  <  1:3567379/149‑1 (MQ=255)
                                                                                                                                      tttCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTaa                <  1:3433847/149‑1 (MQ=255)
                                                                                                                                         cATTGCCTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAAttt             >  2:908995/1‑149 (MQ=255)
                                                                                                                                           ttGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATttgt           >  2:3679828/1‑149 (MQ=255)
                                                                                                                                             gCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATttgttg         >  1:1132776/1‑149 (MQ=255)
                                                                                                                                               cTGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGaa       >  2:2481522/1‑149 (MQ=255)
                                                                                                                                                tGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGaaa      <  2:148443/149‑1 (MQ=255)
                                                                                                                                                tGCTCATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGaaa      <  2:880077/149‑1 (MQ=255)
                                                                                                                                                tGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGaaa      <  1:2077039/149‑1 (MQ=255)
                                                                                                                                                  cTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAata    >  2:1542501/1‑149 (MQ=255)
                                                                                                                                                  cTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAata    >  2:1542495/1‑149 (MQ=255)
                                                                                                                                                    aaTCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATc  >  2:1395271/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TACGACCCAGCAAGGGGATGTGTATACCATTATTCCTGAAGTCACTCTTACTCAATCTTGTCTGTGCAGAGTACAAATATTGTCCCTGCGCGAAGGCAGTTCAGGGCAAAGTCAGACGAAGCAAGAAAAGACCCTTTCATTGCCTGCTAATCAACCCATTGCTTTGACGAAGTTGAGTTTAAATATTTCCCCGGACGATCGGGTGAAAATAGTTGTTACTGTTTCTGATGGACAGTCACTTCATTTATCACAACAATGGCCGCCCTCTTCAGAAAAGTCTTAATTTGTTGAAATATC  >  CP009273/1100467‑1100763

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: