Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 743,797 G→T 100% R214R (CGC→CGA)  ybgO ← putative fimbrial protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273743,7970GT97.3% 625.0 / 6.2 203R214R (CGC→CGA) ybgOputative fimbrial protein
Reads supporting (aligned to +/- strand):  ref base G (4/2);  new base T (106/91);  total (110/93)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 204 total aligned reads displayed.

GTGCATTGATCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACGCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTATTTTG  >  CP009273/743651‑743941
                                                                                                                                                  |                                                                                                                                                
gTGCATTGATCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCAcgc                                                                                                                                                 >  2:1592210/1‑148 (MQ=255)
    aTTGATCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAAt                                                                                                                                             <  1:3148496/148‑1 (MQ=255)
      tGATCGTTTTGGCTCTTGGTCGTTTTAATAATGAACCTTTTAGACATTTTGTGGCGCCGTCTATACAGCAAAATAAATTTAAAAAAATCAAAAACAAGTGTTTACGGGGAAATTGAAAAAATGACGCTGCAAACCAGCACTCGAATGt                                                                                                                                           >  1:1528350/1‑148 (MQ=255)
         tCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGtttt                                                                                                                                        <  2:3132921/148‑1 (MQ=255)
         tCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGtttt                                                                                                                                        <  2:3446518/148‑1 (MQ=255)
            ttttGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTACa                                                                                                                                     >  2:693549/1‑148 (MQ=255)
             tttGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAg                                                                                                                                    >  1:1466847/1‑148 (MQ=255)
                 gCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATACAGCAAATTAAATTTACCAAAATCAATAACCTGTGTTTACGGTGAAATTGAAAAATTTACGCTTCAATCCAGCACTCGAATGTTTTACAGACCa                                                                                                                                >  2:311695/1‑148 (MQ=255)
                  ctctTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAg                                                                                                                               <  2:137135/148‑1 (MQ=255)
                    ctTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTa                                                                                                                             <  2:702510/148‑1 (MQ=255)
                        gTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAatat                                                                                                                         >  2:476073/1‑148 (MQ=255)
                        gTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTAACAAAATCAATAACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATACAGCACTCGAATGTTTTCCAGACCAGTAatat                                                                                                                         >  1:1136375/1‑148 (MQ=255)
                          cGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGa                                                                                                                       >  1:1340613/1‑148 (MQ=255)
                             tttAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGatat                                                                                                                    >  2:2816481/1‑148 (MQ=255)
                              ttAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATAtt                                                                                                                   <  1:2562840/148‑1 (MQ=255)
                              ttAATACTGAACGTTTTCGACATTCTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTTTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCAAGCACTCGAATGTTTTCCAGACCAGTAATATGATAtt                                                                                                                   >  1:2789138/1‑148 (MQ=255)
                                aaTACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTc                                                                                                                 <  2:3540756/148‑1 (MQ=255)
                                  tACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAg                                                                                                               >  2:1927738/1‑148 (MQ=255)
                                   aCTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAgg                                                                                                              <  1:3078304/148‑1 (MQ=255)
                                    cTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGt                                                                                                             <  2:2197282/148‑1 (MQ=255)
                                     tGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGtt                                                                                                            >  1:461147/1‑148 (MQ=255)
                                     tGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGtt                                                                                                            >  1:3126097/1‑148 (MQ=255)
                                      gAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGttt                                                                                                           <  1:1947915/148‑1 (MQ=255)
                                       aaCGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGtttt                                                                                                          >  2:660104/1‑148 (MQ=255)
                                        aCGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTg                                                                                                         >  1:3188867/1‑148 (MQ=255)
                                        aCGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTg                                                                                                         <  1:742299/148‑1 (MQ=255)
                                         cGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTgg                                                                                                        <  1:794979/148‑1 (MQ=255)
                                          gTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTggg                                                                                                       <  1:1547427/148‑1 (MQ=255)
                                            tttCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCa                                                                                                     >  2:2951624/1‑148 (MQ=255)
                                            tttCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCa                                                                                                     <  1:292220/148‑1 (MQ=255)
                                                gACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAg                                                                                                 <  1:2080374/148‑1 (MQ=255)
                                                   aTTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGt                                                                                              <  2:1491485/148‑1 (MQ=255)
                                                   aTTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGt                                                                                              >  2:3384908/1‑148 (MQ=255)
                                                    ttGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTc                                                                                             <  2:1486613/148‑1 (MQ=255)
                                                     tgtgtgGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCt                                                                                            >  2:760717/1‑148 (MQ=255)
                                                        gtgGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCg                                                                                         <  1:1860543/148‑1 (MQ=255)
                                                        gtgGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCg                                                                                         <  1:56596/148‑1 (MQ=255)
                                                         tgGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGt                                                                                        <  1:1116699/148‑1 (MQ=255)
                                                         tgGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGt                                                                                        >  2:620914/1‑148 (MQ=255)
                                                         tgGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGt                                                                                        <  2:892337/148‑1 (MQ=255)
                                                           gCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTa                                                                                      <  2:2949526/148‑1 (MQ=255)
                                                                 tataTCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTa                                                                                <  1:449317/148‑1 (MQ=255)
                                                                   tatCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTAcc                                                                              >  2:1390500/1‑148 (MQ=255)
                                                                    atCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCg                                                                             >  2:1767944/1‑148 (MQ=255)
                                                                    atCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCg                                                                             >  2:2301079/1‑148 (MQ=255)
                                                                           aaattaaattTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGt                                                                      >  1:2831515/1‑148 (MQ=255)
                                                                            aattaaattTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTc                                                                     >  1:545876/1‑148 (MQ=255)
                                                                             attaaattTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCa                                                                    >  1:777439/1‑148 (MQ=255)
                                                                              ttaaattTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCaa                                                                   >  1:1542638/1‑148 (MQ=255)
                                                                               taaattTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCaaa                                                                  >  1:971867/1‑148 (MQ=255)
                                                                                  attTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTg                                                               >  2:1155644/1‑148 (MQ=255)
                                                                                  attTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACGCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTg                                                               >  1:1086218/1‑148 (MQ=255)
                                                                                    ttACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCa                                                             <  2:821147/148‑1 (MQ=255)
                                                                                       ccAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCacga                                                          <  2:2589788/148‑1 (MQ=255)
                                                                                        cAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCacgtc                                                         >  2:1794127/1‑146 (MQ=255)
                                                                                        cAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCacgac                                                         >  2:1227090/1‑148 (MQ=255)
                                                                                         aaaaTCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCacgacg                                                        >  2:1997843/1‑148 (MQ=255)
                                                                                          aaaTCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCAcgacga                                                       >  2:687167/1‑148 (MQ=255)
                                                                                           aatcaatcACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGat                                                      >  1:2009021/1‑148 (MQ=255)
                                                                                             tcaatcaCCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGatat                                                    <  1:67698/148‑1 (MQ=255)
                                                                                              caatcaCCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGAtata                                                   >  2:1089226/1‑148 (MQ=255)
                                                                                               aatcaCCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACGCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAg                                                  <  1:1592210/148‑1 (MQ=255)
                                                                                                 tcaCCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTc                                                >  1:1699019/1‑148 (MQ=255)
                                                                                                     cTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTa                                            >  2:933059/1‑148 (MQ=255)
                                                                                                     cTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTa                                            >  2:2477460/1‑148 (MQ=255)
                                                                                                     cTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTa                                            >  2:2239739/1‑148 (MQ=255)
                                                                                                      tgtgTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTaa                                           >  1:3238858/1‑148 (MQ=255)
                                                                                                        tgtTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGc                                         >  2:1725264/1‑148 (MQ=255)
                                                                                                         gtTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCg                                        >  1:1818210/1‑148 (MQ=255)
                                                                                                           ttCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATATTCGCTAAGCGTa                                      >  1:3286748/1‑148 (MQ=255)
                                                                                                           ttCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTa                                      >  1:3011332/1‑148 (MQ=255)
                                                                                                             ccGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACt                                    <  2:1017254/148‑1 (MQ=255)
                                                                                                                 tGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGa                                >  1:1289791/1‑148 (MQ=255)
                                                                                                                    aaTGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATaa                             <  1:2441968/148‑1 (MQ=255)
                                                                                                                     aTGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAAt                            <  2:2675272/148‑1 (MQ=255)
                                                                                                                     aTGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAAt                            <  2:2066553/148‑1 (MQ=255)
                                                                                                                     aTGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAAt                            >  2:1379493/1‑148 (MQ=255)
                                                                                                                      tGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATc                           <  1:1238042/148‑1 (MQ=255)
                                                                                                                      tGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATc                           <  1:236141/148‑1 (MQ=255)
                                                                                                                        gAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCAt                         <  2:3197362/148‑1 (MQ=255)
                                                                                                                         aaaaaTTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATg                        <  2:2627417/148‑1 (MQ=255)
                                                                                                                         aaaaaTTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATg                        <  2:2759858/148‑1 (MQ=255)
                                                                                                                              ttGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATgcggcg                   >  1:577744/1‑148 (MQ=255)
                                                                                                                                 aCGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCa                >  1:503422/1‑148 (MQ=255)
                                                                                                                                  cGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCAt               >  1:1685447/1‑148 (MQ=255)
                                                                                                                                  cGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCAt               >  1:879545/1‑148 (MQ=255)
                                                                                                                                  cGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCAt               >  1:2382177/1‑148 (MQ=255)
                                                                                                                                    cTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCAttt             <  1:2619170/148‑1 (MQ=255)
                                                                                                                                       cAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTtctc          <  1:2221011/148‑1 (MQ=255)
                                                                                                                                       cAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTtctc          >  2:1702805/1‑148 (MQ=255)
                                                                                                                                         aTCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGt        <  1:693549/148‑1 (MQ=255)
                                                                                                                                          tCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTa       >  2:38423/1‑148 (MQ=255)
                                                                                                                                           ccAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTAt      <  2:3281364/148‑1 (MQ=255)
                                                                                                                                           ccAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTAt      >  1:2726685/1‑148 (MQ=255)
                                                                                                                                           ccAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTAt      >  1:357500/1‑148 (MQ=255)
                                                                                                                                            cAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTAtt     <  1:2659357/148‑1 (MQ=255)
                                                                                                                                            cAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTAtt     <  2:1290696/148‑1 (MQ=255)
                                                                                                                                            cAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTAtt     <  2:2403012/148‑1 (MQ=255)
                                                                                                                                             aGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCTCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTAttt    <  2:1231932/148‑1 (MQ=255)
                                                                                                                                               cACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTATTTTg  >  1:876098/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GTGCATTGATCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACGCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTATTTTG  >  CP009273/743651‑743941

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: ‑