Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,789,381 C→A 100% L15I (CTC→ATC)  waaC → ADP‑heptose:LPS heptosyl transferase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,789,3810CA95.0% 507.6 / 15.7 171L15I (CTC→ATC) waaCADP‑heptose:LPS heptosyl transferase I
Reads supporting (aligned to +/- strand):  ref base C (4/5);  new base A (93/69);  total (97/74)
Fisher's exact test for biased strand distribution p-value = 5.03e-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 172 total aligned reads displayed.

AGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGG  >  CP009273/3789236‑3789524
                                                                                                                                                 |                                                                                                                                               
aGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTtatc                                                                                                                                               >  1:2937934/1‑148 (MQ=255)
 ggTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTcc                                                                                                                                              >  1:1708574/1‑148 (MQ=255)
  gTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCa                                                                                                                                             >  2:3453964/1‑148 (MQ=255)
   tGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCAt                                                                                                                                            >  2:982056/1‑148 (MQ=255)
   tGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCAt                                                                                                                                            >  1:846366/1‑148 (MQ=255)
      cGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACg                                                                                                                                         <  2:1988327/148‑1 (MQ=255)
      cGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACg                                                                                                                                         <  2:191906/148‑1 (MQ=255)
         tGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTg                                                                                                                                      <  1:1085704/148‑1 (MQ=255)
         tGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTg                                                                                                                                      <  1:1085466/148‑1 (MQ=255)
         tGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTg                                                                                                                                      <  2:1828504/148‑1 (MQ=255)
             gAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCg                                                                                                                                  >  2:966607/1‑148 (MQ=255)
              aGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGc                                                                                                                                 >  2:3131031/1‑148 (MQ=255)
              aGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGc                                                                                                                                 >  2:2998747/1‑148 (MQ=255)
               gggTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGca                                                                                                                                >  1:467664/1‑148 (MQ=255)
                ggTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGcac                                                                                                                               >  2:1324016/1‑148 (MQ=255)
                ggTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGcac                                                                                                                               >  2:1323869/1‑148 (MQ=255)
                     tCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGcactcact                                                                                                                          >  1:1514132/1‑148 (MQ=255)
                       aCCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGa                                                                                                                        <  2:2102462/148‑1 (MQ=255)
                       aCCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGa                                                                                                                        <  2:1449842/148‑1 (MQ=255)
                          agagCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGc                                                                                                                     >  1:3239193/1‑148 (MQ=255)
                              cTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCcag                                                                                                                 <  2:1064741/148‑1 (MQ=255)
                                tAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCcagca                                                                                                               <  2:727192/148‑1 (MQ=255)
                                 aaTCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCcagcag                                                                                                              >  2:761032/1‑148 (MQ=255)
                                 aaTCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCcagcag                                                                                                              >  2:2523309/1‑148 (MQ=255)
                                 aaTCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCcagcag                                                                                                              >  2:761002/1‑148 (MQ=255)
                                  aTCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAgg                                                                                                             >  1:464597/1‑148 (MQ=255)
                                   tCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGc                                                                                                            >  1:247059/1‑148 (MQ=255)
                                            cTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCTATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAgg                                                                                                   >  1:2414179/1‑148 (MQ=255)
                                             tCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAggg                                                                                                  <  2:1167771/148‑1 (MQ=255)
                                              ccccAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGa                                                                                                 <  1:2741732/148‑1 (MQ=255)
                                              ccccAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGa                                                                                                 >  1:3097883/1‑148 (MQ=255)
                                                ccAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGAtt                                                                                               >  1:729339/1‑148 (MQ=255)
                                                  aGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTaa                                                                                             <  1:1321818/148‑1 (MQ=255)
                                                   gcgcGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAg                                                                                            >  2:2492124/1‑148 (MQ=255)
                                                   gcgcGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAg                                                                                            <  1:2199932/148‑1 (MQ=255)
                                                     gcgTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGtt                                                                                          <  1:3353473/148‑1 (MQ=255)
                                                        tACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGa                                                                                       <  1:587905/148‑1 (MQ=255)
                                                        tACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGa                                                                                       <  2:2689897/148‑1 (MQ=255)
                                                        tACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGa                                                                                       >  2:1749980/1‑148 (MQ=255)
                                                        tACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGa                                                                                       >  2:2059580/1‑148 (MQ=255)
                                                           tGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTg                                                                                    <  2:207666/148‑1 (MQ=255)
                                                            ggAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTgg                                                                                   <  1:2799762/148‑1 (MQ=255)
                                                              aagaaCTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTCAGTTTGACTGggt                                                                                 >  1:886305/1‑148 (MQ=255)
                                                               agaaCTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGggtg                                                                                <  2:358336/148‑1 (MQ=255)
                                                                 aaCTCAACGCGCTATTGTTACAAGAGGAAGCCTGACTGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGggtggt                                                                              <  1:1560030/148‑1 (MQ=255)
                                                                 aaCTCAACGCGCTATTGTTACAAGAGGAAGACTGACTGATGCGGGTTTTGATCGTTAAAACATCGTCGATGTGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCATCAGGCAATCCCAGGGATTAAGTTTGACTTggtggt                                                                              <  2:2229378/148‑1 (MQ=255)
                                                                  aCTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGgtggtg                                                                             <  2:3069037/148‑1 (MQ=255)
                                                                   cTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGtggtgg                                                                            >  1:1604294/1‑148 (MQ=255)
                                                                       aCGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGgaaga                                                                        <  1:2059580/148‑1 (MQ=255)
                                                                       aCGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGgaaga                                                                        >  2:1820512/1‑148 (MQ=255)
                                                                        cgcgCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGgaagaa                                                                       >  2:2631439/1‑148 (MQ=255)
                                                                            cTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGt                                                                   <  1:3418389/148‑1 (MQ=255)
                                                                            cTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGGTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGt                                                                   <  1:3419876/148‑1 (MQ=255)
                                                                             tATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGtt                                                                  >  2:2011251/1‑148 (MQ=255)
                                                                                   tACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGcaca                                                            <  1:2947164/148‑1 (MQ=255)
                                                                                     cAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGa                                                          <  1:3184403/148‑1 (MQ=255)
                                                                                     cAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGa                                                          <  1:3184507/148‑1 (MQ=255)
                                                                                       agagGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAtt                                                        >  2:3225625/1‑148 (MQ=255)
                                                                                        gagGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAttc                                                       >  1:3301354/1‑148 (MQ=255)
                                                                                         agGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAttcc                                                      <  1:1068849/148‑1 (MQ=255)
                                                                                         agGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAttcc                                                      <  2:500607/148‑1 (MQ=255)
                                                                                           gAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAttcctt                                                    <  1:576983/148‑1 (MQ=255)
                                                                                            aaGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAttccttc                                                   <  1:2822626/148‑1 (MQ=255)
                                                                                            aaGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAttccttc                                                   <  1:1047233/148‑1 (MQ=255)
                                                                                             aGCCTGACGGATGCGGGTTTTGATCGTTAATACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAttccttcc                                                  <  2:1128494/148‑1 (MQ=255)
                                                                                             aGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGAttccttcc                                                  <  2:3275262/148‑1 (MQ=255)
                                                                                              gCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATtccttcct                                                 >  1:1746168/1‑148 (MQ=255)
                                                                                               ccTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTg                                                >  2:1196407/1‑148 (MQ=255)
                                                                                                    cGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACg                                           <  1:2520414/148‑1 (MQ=255)
                                                                                                    cGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACg                                           <  2:3097461/148‑1 (MQ=255)
                                                                                                    cGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACg                                           <  2:3097469/148‑1 (MQ=255)
                                                                                                     ggATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGc                                          >  1:3002649/1‑148 (MQ=255)
                                                                                                      gATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCt                                         >  2:393929/1‑148 (MQ=255)
                                                                                                       aTGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTg                                        >  1:1880567/1‑148 (MQ=255)
                                                                                                       aTGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTg                                        >  2:133692/1‑148 (MQ=255)
                                                                                                            ggTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGtt                                   <  2:1193136/148‑1 (MQ=255)
                                                                                                             gTTTTGATCGTTAAAACATCGTCGATGGGCGATTTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGTGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTg                                  <  1:3040352/148‑1 (MQ=255)
                                                                                                               tttGATCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAg                                >  1:2941411/1‑148 (MQ=255)
                                                                                                                   aTCGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAg                            >  1:1277395/1‑148 (MQ=255)
                                                                                                                   aTCGTTAAAACATCGTCGATGGGCGATGTTATCCATAAGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAg                            >  2:2591629/1‑148 (MQ=255)
                                                                                                                     cGTTAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGtt                          <  2:2937934/148‑1 (MQ=255)
                                                                                                                        tAAAACATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTAtt                       >  1:994641/1‑148 (MQ=255)
                                                                                                                        tAAAACATCGTCGATGGGCGATGTTATCCATAAGTTGCCCGCACTCACTGATGCCCAGCAGGCAATACCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGACGTTGAGCGAGTTAtt                       >  1:1836277/1‑148 (MQ=255)
                                                                                                                          aaaCATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTcc                     <  2:371627/148‑1 (MQ=255)
                                                                                                                           aaCATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCt                    >  1:3052384/1‑148 (MQ=255)
                                                                                                                            aCATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCtg                   <  2:125617/148‑1 (MQ=255)
                                                                                                                             cATCGTCGATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCtgt                  >  2:1456174/1‑148 (MQ=255)
                                                                                                                               tcgtcgATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTgg                >  2:1115018/1‑148 (MQ=255)
                                                                                                                                  tcgATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCaa             >  1:787911/1‑148 (MQ=255)
                                                                                                                                   cgATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAAt            >  1:1498529/1‑148 (MQ=255)
                                                                                                                                   cgATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAAt            >  2:400833/1‑148 (MQ=255)
                                                                                                                                   cgATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAAt            >  2:400834/1‑148 (MQ=255)
                                                                                                                                   cgATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAAt            >  1:409849/1‑148 (MQ=255)
                                                                                                                                   cgATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAAt            >  1:1050993/1‑148 (MQ=255)
                                                                                                                                    gATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATa           <  2:3529612/148‑1 (MQ=255)
                                                                                                                                    gATGGGCGATGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATa           <  2:1604294/148‑1 (MQ=255)
                                                                                                                                        ggCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTc       >  1:2904239/1‑148 (MQ=255)
                                                                                                                                            aTGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTg   >  2:2812984/1‑148 (MQ=255)
                                                                                                                                             tGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTgg  >  2:1030630/1‑148 (MQ=255)
                                                                                                                                             tGTTATCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTgg  >  1:516061/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
AGGTGACGCTGCGGAGGGTTATCACCAGAGCTTAATCGACATTACTCCCCAGCGCGTACTGGAAGAACTCAACGCGCTATTGTTACAAGAGGAAGCCTGACGGATGCGGGTTTTGATCGTTAAAACATCGTCGATGGGCGATGTTCTCCATACGTTGCCCGCACTCACTGATGCCCAGCAGGCAATCCCAGGGATTAAGTTTGACTGGGTGGTGGAAGAAGGGTTCGCACAGATTCCTTCCTGGCACGCTGCCGTTGAGCGAGTTATTCCTGTGGCAATACGTCGCTGG  >  CP009273/3789236‑3789524

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: