Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,097,943 T→A 41.5% N19Y (AAT→TAT)  csgE ← curlin secretion specificity factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,097,9430TA41.5% 174.4 / inf 275N19Y (AAT→TAT) csgEcurlin secretion specificity factor
Reads supporting (aligned to +/- strand):  ref base T (71/90);  new base A (66/48);  total (137/138)
Fisher's exact test for biased strand distribution p-value = 2.78e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 278 total aligned reads displayed.

TGCACTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATGCGTT  >  CP009273/1097797‑1098089
                                                                                                                                                  |                                                                                                                                                  
tGCACTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTc                                                                                                                                                   <  2:2013720/148‑1 (MQ=255)
 gCACTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATAcc                                                                                                                                                  <  1:2903227/148‑1 (MQ=255)
 gCACTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATAcc                                                                                                                                                  <  2:1169153/148‑1 (MQ=255)
 gCACTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATAcc                                                                                                                                                  <  1:124807/148‑1 (MQ=255)
  cACTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTccc                                                                                                                                                 <  1:606927/148‑1 (MQ=255)
   aCTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCt                                                                                                                                                >  2:2719994/1‑148 (MQ=255)
    cTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTg                                                                                                                                               >  1:2932/1‑148 (MQ=255)
    cTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTg                                                                                                                                               >  1:1646022/1‑148 (MQ=255)
     tGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTgc                                                                                                                                              <  2:2613340/148‑1 (MQ=255)
      gggCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTgcg                                                                                                                                             >  1:2078108/1‑148 (MQ=255)
       ggCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTgcgg                                                                                                                                            >  2:599877/1‑148 (MQ=255)
            ttCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAAc                                                                                                                                       >  1:2344590/1‑148 (MQ=255)
            ttCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAAc                                                                                                                                       <  1:707619/148‑1 (MQ=255)
               aTTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGa                                                                                                                                    <  2:3435900/148‑1 (MQ=255)
               aTTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGa                                                                                                                                    >  2:3485849/1‑148 (MQ=255)
                 tAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGaaa                                                                                                                                  <  1:3146409/148‑1 (MQ=255)
                   aTCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAAtt                                                                                                                                <  1:807924/148‑1 (MQ=255)
                   aTCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAAtt                                                                                                                                <  1:739914/148‑1 (MQ=255)
                    tCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTc                                                                                                                               >  1:3363259/1‑148 (MQ=255)
                    tCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTc                                                                                                                               <  1:2619939/148‑1 (MQ=255)
                      gTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTg                                                                                                                             <  1:2531395/148‑1 (MQ=255)
                       ttAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTgc                                                                                                                            >  1:2092947/1‑148 (MQ=255)
                       ttAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTgc                                                                                                                            >  1:1365596/1‑148 (MQ=255)
                          aGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTgccgc                                                                                                                         >  2:3166420/1‑148 (MQ=255)
                            ttACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCca                                                                                                                       >  2:452611/1‑148 (MQ=255)
                                  gTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATcc                                                                                                                 <  1:1288096/148‑1 (MQ=255)
                                  gTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATcc                                                                                                                 <  2:3220865/148‑1 (MQ=255)
                                  gTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATcc                                                                                                                 >  2:715132/1‑148 (MQ=255)
                                    ataGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAg                                                                                                               >  1:407601/1‑148 (MQ=255)
                                      aGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCg                                                                                                             <  1:2569419/148‑1 (MQ=255)
                                          aCTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTaaa                                                                                                         >  2:1949971/1‑148 (MQ=255)
                                             ttCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATaa                                                                                                      <  1:178139/148‑1 (MQ=255)
                                                     tATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCAt                                                                                              <  2:3375221/148‑1 (MQ=255)
                                                     tATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCAt                                                                                              <  1:1866832/148‑1 (MQ=255)
                                                        cACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGc                                                                                           <  1:1020827/148‑1 (MQ=255)
                                                          cTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCtt                                                                                         >  1:3457615/1‑148 (MQ=255)
                                                            aaaGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTa                                                                                       <  2:1256186/148‑1 (MQ=255)
                                                            aaaGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTa                                                                                       >  1:1388113/1‑148 (MQ=255)
                                                                gCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAAACCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGc                                                                                   >  2:2933305/1‑148 (MQ=255)
                                                                 cTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGcc                                                                                  <  2:1833851/148‑1 (MQ=255)
                                                                 cTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGcc                                                                                  >  2:77960/1‑148 (MQ=255)
                                                                   cGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTg                                                                                >  1:1880306/1‑148 (MQ=255)
                                                                      tAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAgg                                                                             <  1:3368104/148‑1 (MQ=255)
                                                                      tAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAgg                                                                             <  2:3299299/148‑1 (MQ=255)
                                                                       aaaaaTCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGt                                                                            >  2:3166769/1‑148 (MQ=255)
                                                                          aaTCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTAt                                                                         >  1:3363356/1‑148 (MQ=255)
                                                                          aaTCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTAt                                                                         >  1:1450182/1‑148 (MQ=255)
                                                                              aTGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGtt                                                                     <  2:1040172/148‑1 (MQ=255)
                                                                              aTGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGtt                                                                     <  2:1545643/148‑1 (MQ=255)
                                                                                    aaTAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCa                                                               >  2:3546079/1‑148 (MQ=255)
                                                                                    aaTAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCa                                                               >  2:837312/1‑148 (MQ=255)
                                                                                       aGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGa                                                            >  1:1018828/1‑148 (MQ=255)
                                                                                        gATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGTCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGaa                                                           <  2:2237904/148‑1 (MQ=255)
                                                                                        gATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGaa                                                           >  2:328600/1‑148 (MQ=255)
                                                                                         aTGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGaaa                                                          >  1:28140/1‑148 (MQ=255)
                                                                                         aTGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGaaa                                                          <  2:892553/148‑1 (MQ=255)
                                                                                          tGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAAc                                                         >  1:1444684/1‑148 (MQ=255)
                                                                                          tGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAAc                                                         >  1:1445212/1‑148 (MQ=255)
                                                                                            aaaCAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCg                                                       >  2:2056139/1‑148 (MQ=255)
                                                                                             aaCAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGACGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGc                                                      >  1:3222866/1‑148 (MQ=255)
                                                                                              aCAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCt                                                     >  1:384585/1‑148 (MQ=255)
                                                                                                  tATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTtgtg                                                 >  2:1396281/1‑148 (MQ=255)
                                                                                                   aTGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTgtgt                                                <  1:2739562/148‑1 (MQ=255)
                                                                                                    tGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTc                                               <  1:325094/148‑1 (MQ=255)
                                                                                                       tCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCgg                                            <  1:893764/148‑1 (MQ=255)
                                                                                                       tCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCgg                                            <  2:2389734/148‑1 (MQ=255)
                                                                                                           ttAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGtttt                                        >  1:2371552/1‑148 (MQ=255)
                                                                                                           ttAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGtttt                                        >  1:3063778/1‑148 (MQ=255)
                                                                                                           ttAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGtttt                                        <  1:3166452/148‑1 (MQ=255)
                                                                                                             aGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGtttttt                                      <  1:950707/148‑1 (MQ=255)
                                                                                                              gCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACACTCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTa                                     <  1:2719994/148‑1 (MQ=255)
                                                                                                               cAATCCCGGGACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTAc                                    <  2:2078130/148‑1 (MQ=255)
                                                                                                                 aTCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACgg                                  <  1:2742995/148‑1 (MQ=255)
                                                                                                                    ccGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTa                               <  2:3418100/148‑1 (MQ=255)
                                                                                                                     cGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTAt                              <  2:1342843/148‑1 (MQ=255)
                                                                                                                     cGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTAt                              <  1:1185629/148‑1 (MQ=255)
                                                                                                                        gACTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTAtctt                           >  2:3390212/1‑148 (MQ=255)
                                                                                                                         aCTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTAtcttc                          <  1:3527400/148‑1 (MQ=255)
                                                                                                                         aCTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTAtcttc                          <  2:2869309/148‑1 (MQ=255)
                                                                                                                          cTTCTACCTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTAtcttct                         <  2:448365/148‑1 (MQ=255)
                                                                                                                            tCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTg                       >  1:2515340/1‑148 (MQ=255)
                                                                                                                               aCCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTGGTGTCCGGTTTTTTACGGCTATCTTCTTGaaa                    >  2:1502385/1‑148 (MQ=255)
                                                                                                                                ccTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGaaaa                   <  2:1646022/148‑1 (MQ=255)
                                                                                                                                 cTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAg                  <  1:250818/148‑1 (MQ=255)
                                                                                                                                 cTCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAg                  <  1:3564261/148‑1 (MQ=255)
                                                                                                                                  tCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGa                 <  1:458834/148‑1 (MQ=255)
                                                                                                                                  tCAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGa                 <  2:16102/148‑1 (MQ=255)
                                                                                                                                   cAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGAt                >  2:37184/1‑148 (MQ=255)
                                                                                                                                   cAACGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGAt                >  1:2260856/1‑148 (MQ=255)
                                                                                                                                    aaCGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGAtt               <  1:196196/148‑1 (MQ=255)
                                                                                                                                      cGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATtat             >  1:2412658/1‑148 (MQ=255)
                                                                                                                                      cGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATtat             >  1:1585421/1‑148 (MQ=255)
                                                                                                                                      cGGCGTGAAGATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATtat             >  1:1754046/1‑148 (MQ=255)
                                                                                                                                        gCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATaa           >  2:1418170/1‑148 (MQ=255)
                                                                                                                                            gAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGAt       <  2:3068026/148‑1 (MQ=255)
                                                                                                                                               gATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATGCg    <  2:2525980/148‑1 (MQ=255)
                                                                                                                                               gATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGATATCTTCTTGAAAAGATTATAAAGATGCg    >  2:2372310/1‑148 (MQ=255)
                                                                                                                                               gATACCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATGCg    >  2:673067/1‑148 (MQ=255)
                                                                                                                                                 ttCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATGCGtt  <  1:445448/148‑1 (MQ=255)
                                                                                                                                                 ttCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATGCGtt  >  1:113155/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TGCACTGGGCCTTTCATTAATCGTTAAGTTACCCGTATAGTCACTTTCCCATTTATCACTAAAGGCTCGGTAAAAATCATGGCCAATAGATGAAACAGTATGGTCAGTTAGCAATCCCGGGACTTCTACCTCAACGGCGTGAAGATTCCCTGCGGCGAACAGAAATTCTGCCGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATGCGTT  >  CP009273/1097797‑1098089

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: