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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273743,7970GT84.4% 396.5 / 82.7 185R214R (CGC→CGAybgOputative fimbrial protein
Reads supporting (aligned to +/- strand):  ref base G (9/20);  new base T (63/93);  total (72/113)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 185 total aligned reads displayed.

CATTGATCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACGCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTATTTTGA  >  CP009273/743654‑743942
                                                                                                                                               |                                                                                                                                                 
cATTGATCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGaa                                                                                                                                               <  2:2032094/148‑1 (MQ=255)
cATTGATCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGaa                                                                                                                                               <  2:1661484/148‑1 (MQ=255)
cATTGATCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGaa                                                                                                                                               <  2:839497/148‑1 (MQ=255)
 aTTGATCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAAt                                                                                                                                              <  1:731947/148‑1 (MQ=255)
   tGATCTTTTTGTCTCTTGTTAGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGt                                                                                                                                            <  2:3170899/148‑1 (MQ=255)
   tGATCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGt                                                                                                                                            <  2:2624190/148‑1 (MQ=255)
     aTCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGttt                                                                                                                                          >  2:19117/1‑148 (MQ=255)
     aTCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGttt                                                                                                                                          <  1:1422120/148‑1 (MQ=255)
     aTCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACGCGAATGttt                                                                                                                                          <  2:1326863/148‑1 (MQ=255)
      tCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGtttt                                                                                                                                         <  1:807851/148‑1 (MQ=255)
      tCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATAGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGtttt                                                                                                                                         >  2:460357/1‑148 (MQ=255)
       cGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTc                                                                                                                                        >  2:2955272/1‑148 (MQ=255)
          tttGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAg                                                                                                                                     <  2:3300175/148‑1 (MQ=255)
             ggCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGAcc                                                                                                                                  <  1:86966/148‑1 (MQ=255)
                    ggTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAata                                                                                                                           <  2:73561/148‑1 (MQ=255)
                      tCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATg                                                                                                                         <  1:668743/148‑1 (MQ=255)
                      tCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATg                                                                                                                         <  1:3055531/148‑1 (MQ=255)
                      tCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATg                                                                                                                         <  1:2329912/148‑1 (MQ=255)
                       cGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGa                                                                                                                        >  1:1200527/1‑148 (MQ=255)
                          tttAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACGCGAATGTTTTCCAGACCAGTAATATGatat                                                                                                                     >  2:180580/1‑148 (MQ=255)
                           ttAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATAtt                                                                                                                    <  1:1017131/148‑1 (MQ=255)
                           ttAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATAtt                                                                                                                    <  1:1016358/148‑1 (MQ=255)
                             aaTACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTc                                                                                                                  <  1:1866903/148‑1 (MQ=255)
                             aaTACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTc                                                                                                                  >  2:1770740/1‑148 (MQ=255)
                              aTACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACGCGAATGTTTTCCAGACCAGTAATATGATATTTCa                                                                                                                 <  2:2247508/148‑1 (MQ=255)
                                aCTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAgg                                                                                                               <  1:1405051/148‑1 (MQ=255)
                                aCTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAgg                                                                                                               >  1:3022681/1‑148 (MQ=255)
                                aCTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACGCGAATGTTTTCCAGACCAGTAATATGATATTTCAgg                                                                                                               <  1:2501586/148‑1 (MQ=255)
                                  tGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGtt                                                                                                             <  2:251567/148‑1 (MQ=255)
                                  tGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACGCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGtt                                                                                                             <  2:3203587/148‑1 (MQ=255)
                                     aCGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACGCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTg                                                                                                          <  1:3144540/148‑1 (MQ=255)
                                         tttCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCa                                                                                                      >  2:1143410/1‑148 (MQ=255)
                                           tCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGt                                                                                                    >  1:2951798/1‑148 (MQ=255)
                                            cGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTa                                                                                                   >  2:1276187/1‑148 (MQ=255)
                                              aCATTGTGTGGCGACGTATATCCAGCAACTTAAATTTACCACAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAgg                                                                                                 <  2:3307024/148‑1 (MQ=255)
                                               cATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAggg                                                                                                <  2:3259875/148‑1 (MQ=255)
                                                 ttGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTc                                                                                              <  1:1245219/148‑1 (MQ=255)
                                                 ttGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACGCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTc                                                                                              <  1:1188745/148‑1 (MQ=255)
                                                      tgGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGt                                                                                         <  2:2483526/148‑1 (MQ=255)
                                                        gCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTa                                                                                       <  2:413665/148‑1 (MQ=255)
                                                         cGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTaa                                                                                      <  1:173417/148‑1 (MQ=255)
                                                         cGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTaa                                                                                      >  2:831721/1‑148 (MQ=255)
                                                         cGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTaa                                                                                      >  2:789247/1‑148 (MQ=255)
                                                              tataTCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTa                                                                                 <  1:1654416/148‑1 (MQ=255)
                                                              tataTCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTa                                                                                 >  1:3276158/1‑148 (MQ=255)
                                                                 atCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCg                                                                              >  2:2115894/1‑148 (MQ=255)
                                                                 atCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACGCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCg                                                                              <  1:2147292/148‑1 (MQ=255)
                                                                  tCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCgg                                                                             >  1:1890290/1‑148 (MQ=255)
                                                                    cAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCa                                                                           <  2:2828550/148‑1 (MQ=255)
                                                                     aGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACGCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCAc                                                                          >  2:1204462/1‑148 (MQ=255)
                                                                       cAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACGCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCg                                                                        >  1:2558683/1‑148 (MQ=255)
                                                                        aaattaaattTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGt                                                                       <  2:2054591/148‑1 (MQ=255)
                                                                          attaaattTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCa                                                                     >  2:3186593/1‑148 (MQ=255)
                                                                            taaattTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCaaa                                                                   >  2:943468/1‑148 (MQ=255)
                                                                            taaattTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCaaa                                                                   >  1:740185/1‑148 (MQ=255)
                                                                            taaattTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACGCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCaaa                                                                   >  2:986447/1‑148 (MQ=255)
                                                                                tttACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCCCCGTCAAATTGc                                                               >  2:2128057/1‑148 (MQ=255)
                                                                                 ttACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCa                                                              >  2:3053448/1‑148 (MQ=255)
                                                                                  tACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCac                                                             <  1:3208535/148‑1 (MQ=255)
                                                                                      aaaaTCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCacgacg                                                         >  2:3206033/1‑148 (MQ=255)
                                                                                          tcaatcaCCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGatat                                                     >  1:320639/1‑148 (MQ=255)
                                                                                          tcaatcaCCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGatat                                                     >  1:1986898/1‑148 (MQ=255)
                                                                                           caatcaCCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGAtata                                                    >  1:2632530/1‑148 (MQ=255)
                                                                                            aatcaCCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAg                                                   <  2:294091/148‑1 (MQ=255)
                                                                                            aatcaCCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAg                                                   <  2:67747/148‑1 (MQ=255)
                                                                                            aatcaCCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACGCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAg                                                   <  1:1754991/148‑1 (MQ=255)
                                                                                              tcaCCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTc                                                 >  2:2087146/1‑148 (MQ=255)
                                                                                               caCCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCg                                                <  2:2737382/148‑1 (MQ=255)
                                                                                                aCCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGc                                               <  2:3004920/148‑1 (MQ=255)
                                                                                                 ccTGTGTTTCCGGTGAAATGGAAAAATTGAAGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCt                                              >  2:3029123/1‑148 (MQ=255)
                                                                                                  cTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTa                                             <  2:907287/148‑1 (MQ=255)
                                                                                                  cTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTa                                             <  2:3135910/148‑1 (MQ=255)
                                                                                                     tgtTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGc                                          >  2:1467183/1‑148 (MQ=255)
                                                                                                        ttCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTa                                       <  1:2197043/148‑1 (MQ=255)
                                                                                                        ttCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACGCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTa                                       >  2:174711/1‑148 (MQ=255)
                                                                                                         tCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACGCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTAc                                      >  1:796882/1‑148 (MQ=255)
                                                                                                          ccGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACt                                     >  2:2133139/1‑148 (MQ=255)
                                                                                                          ccGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACt                                     <  1:147061/148‑1 (MQ=255)
                                                                                                               gAAATGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGAt                                >  2:2298743/1‑148 (MQ=255)
                                                                                                                  aTGGAAAAATTGACGCTGCAATCCAGCACGCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAAt                             <  1:1372177/148‑1 (MQ=255)
                                                                                                                   tGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATc                            >  2:2715086/1‑148 (MQ=255)
                                                                                                                   tGGAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATc                            >  2:2475128/1‑148 (MQ=255)
                                                                                                                     gAAAAATTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCAt                          <  1:2519121/148‑1 (MQ=255)
                                                                                                                      aaaaaTTGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATg                         >  1:408686/1‑148 (MQ=255)
                                                                                                                            tGACGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGcggcgg                   >  1:2698206/1‑148 (MQ=255)
                                                                                                                              aCGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCa                 <  1:2938508/148‑1 (MQ=255)
                                                                                                                              aCGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCa                 <  1:2476338/148‑1 (MQ=255)
                                                                                                                              aCGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCa                 <  2:696808/148‑1 (MQ=255)
                                                                                                                               cGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCAt                <  1:3061580/148‑1 (MQ=255)
                                                                                                                               cGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCAt                <  1:185389/148‑1 (MQ=255)
                                                                                                                               cGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCAt                <  2:2459904/148‑1 (MQ=255)
                                                                                                                               cGCTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCAt                <  2:522958/148‑1 (MQ=255)
                                                                                                                                 cTGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCAttt              >  2:1644767/1‑148 (MQ=255)
                                                                                                                                 cTGAAATCTATCACTCTATTGATTTCAAGACCAGTAATATGATTTTTCATGTTTTGGGCAGTATGTTCTTCGTTAACGCTACCGGCACCGTCAATTTGCACTACGATATAGTCGCTAAGCTTACTTTGATTATCATGAGGCGGCAttt              <  1:971342/148‑1 (MQ=255)
                                                                                                                                  tGCAATCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCTGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTtc             <  1:612199/148‑1 (MQ=255)
                                                                                                                                   gCAATCCAGCACGCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTtct            >  1:1585057/1‑148 (MQ=255)
                                                                                                                                     aaTCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCg          <  2:214342/148‑1 (MQ=255)
                                                                                                                                      aTCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGt         >  2:2303955/1‑148 (MQ=255)
                                                                                                                                       tCCAGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTa        >  2:1521681/1‑148 (MQ=255)
                                                                                                                                          aGCACTCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTAttt     <  2:876961/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
CATTGATCGTTTTGGCTCTTGGTCGTTTTAATACTGAACGTTTTCGACATTGTGTGGCGACGTATATCCAGCAAATTAAATTTACCAAAATCAATCACCTGTGTTTCCGGTGAAATGGAAAAATTGACGCTGCAATCCAGCACGCGAATGTTTTCCAGACCAGTAATATGATATTTCAGGTTTTGGGCAGTAGGGTCTTCGTTAACGCTACCGGCACCGTCAAATTGCACGACGATATAGTCGCTAAGCGTACTTTGATAATCATGCGGCGGCATTTCTCGTATTTTGA  >  CP009273/743654‑743942

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: