Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,335,383 T→A 14.7% N756Y (AAC→TAC)  yfaL ← adhesin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,335,3830TA14.7% 812.2 / 113.5 285N756Y (AAC→TAC) yfaLadhesin
Reads supporting (aligned to +/- strand):  ref base T (128/115);  new base A (21/21);  total (149/136)
Fisher's exact test for biased strand distribution p-value = 8.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 289 total aligned reads displayed.

TGCAGGCTAACCGCACCGGACGTATTGGTCAGATTGCCATTAACCGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTGACTGAATATCCTGATCCGCGCCAGTAACGAACGTTGCCCCATCACCA  >  CP009273/2335235‑2335531
                                                                                                                                                    |                                                                                                                                                    
tGCAGGCTAACCGCACCGGACGTATTGGTCAGATTGCCATTAACCGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGtt                                                                                                                                                      <  1:3421424/149‑1 (MQ=255)
   aGGCTAACCGCACCGGACGTATTGGTCAGATTGCCATTAACCGTCACGTTGGAAAGATCAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTaaa                                                                                                                                                   >  2:1591439/1‑149 (MQ=255)
     gCTAACCGCACCGGACGTATTGGTCAGATTGCCATTAACCGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATa                                                                                                                                                 <  2:2047155/149‑1 (MQ=255)
      cTAACCGCACCGGACGTATTGGTCAGATTGCCATTAACCGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAg                                                                                                                                                <  2:154010/149‑1 (MQ=255)
      cTAACCGCACCGGACGTATTGGTCAGATTGCCATTAACCGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAg                                                                                                                                                >  1:4174076/1‑149 (MQ=255)
        aaCCGCACCGGACGTATTGGTCAGATTGCCATTAACCGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTAAAATAGTg                                                                                                                                              <  1:2570652/149‑1 (MQ=255)
          ccGCACCGGACGTATTGGTCAGATTGCCATTAACCGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTAAAATAGTGAg                                                                                                                                            >  2:2226772/1‑149 (MQ=255)
              aCCGGACGTATTGGTCAGATTGCCATTAACCGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCAt                                                                                                                                        >  2:4449004/1‑149 (MQ=255)
                  gACGTATTGGTCAGATTGCCATTAACCGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTAAAATAGTGAGGCATTAAg                                                                                                                                    <  2:3987535/149‑1 (MQ=255)
                   aCGTATTGGTCAGATTGCCATTAACCGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAgg                                                                                                                                   <  1:1540341/149‑1 (MQ=255)
                   aCGTATTGGTCAGATTGCCATTAACCGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAgg                                                                                                                                   <  1:167052/149‑1 (MQ=255)
                    cGTATTGGTCAGATTGCCATTAACCGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAggg                                                                                                                                  <  1:4502659/149‑1 (MQ=255)
                     gTATTGGTCAGATTGCCATTAACCGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGc                                                                                                                                 >  2:3417388/1‑149 (MQ=255)
                     gTATTGGTCAGATTGCCATTAACCGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTAAAATAGTGAGGCATTAAGGGc                                                                                                                                 <  1:4155858/149‑1 (MQ=255)
                      tATTGGTCAGATTGCCATTAACCGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCg                                                                                                                                >  1:981456/1‑149 (MQ=255)
                      tATTGGTCAGATTGCCATTAACCGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTAAAATAGTGAGGCATTAAGGGCg                                                                                                                                >  2:1365395/1‑149 (MQ=255)
                       aTTGGTCAGATTGCCATTAACCGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGa                                                                                                                               <  2:2759932/149‑1 (MQ=255)
                       aTTGGTCAGATTGCCATTAACCGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTAAAATAGTGAGGCATTAAGGGCGa                                                                                                                               <  2:2080161/149‑1 (MQ=255)
                        ttGGTCAGATTGCCATTAACCGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGAc                                                                                                                              >  2:4753780/1‑149 (MQ=255)
                        ttGGTCAGATTGCCATTAACCGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTAAAATAGTGAGGCATTAAGGGCGAc                                                                                                                              <  2:1944544/149‑1 (MQ=255)
                               gATTGCCATTAACCGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTACATAGTGAGGCATTAAGGGCGACGGAAGTa                                                                                                                       >  2:2471366/1‑149 (MQ=255)
                                agttccATTAACCGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTAt                                                                                                                      <  2:3304293/145‑1 (MQ=255)
                                  tGCCATTAACCGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATcc                                                                                                                    <  2:2782996/149‑1 (MQ=255)
                                    ccATTAACCGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTg                                                                                                                  <  1:3121152/149‑1 (MQ=255)
                                    ccATTAACCGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTg                                                                                                                  >  1:1426870/1‑149 (MQ=255)
                                     cATTAACCGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGc                                                                                                                 >  1:1486085/1‑149 (MQ=255)
                                     cATTAACCGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGc                                                                                                                 <  1:3784533/149‑1 (MQ=255)
                                       ttAACCGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTAAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGccc                                                                                                               <  2:1810205/149‑1 (MQ=255)
                                       ttAACCGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTAAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGccc                                                                                                               <  2:1456405/149‑1 (MQ=255)
                                         aaCCGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCgg                                                                                                             >  1:3232376/1‑149 (MQ=255)
                                          aCCGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGt                                                                                                            >  1:3266347/1‑149 (MQ=255)
                                          aCCGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGt                                                                                                            <  2:4662951/149‑1 (MQ=255)
                                           ccGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTc                                                                                                           <  2:4041406/149‑1 (MQ=255)
                                              tCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAgg                                                                                                        >  2:637922/1‑149 (MQ=255)
                                                aCGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGgcg                                                                                                      >  1:4686818/1‑149 (MQ=255)
                                                  gTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCa                                                                                                    >  1:2889127/1‑149 (MQ=255)
                                                   ttGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCaa                                                                                                   >  1:1907092/1‑149 (MQ=255)
                                                   ttGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCaa                                                                                                   >  1:1006710/1‑149 (MQ=255)
                                                   ttGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCaa                                                                                                   >  1:3014076/1‑149 (MQ=255)
                                                      gAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAAc                                                                                                <  1:4354754/149‑1 (MQ=255)
                                                       aaaGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAAcc                                                                                               >  1:2423362/1‑149 (MQ=255)
                                                            tAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTAcc                                                                                          >  1:3745854/1‑149 (MQ=255)
                                                            tAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTAcc                                                                                          <  1:524036/149‑1 (MQ=255)
                                                            tAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTAcc                                                                                          >  2:1161740/1‑149 (MQ=255)
                                                             aaGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCa                                                                                         >  2:4746360/1‑149 (MQ=255)
                                                             aaGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCa                                                                                         <  1:1225266/149‑1 (MQ=255)
                                                               gTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATATTGAGGCATTAAGTGCGACGGAAGTATCCTGACCGGTCAGGCGCAAAACCGTACCATc                                                                                       >  1:1150425/1‑149 (MQ=255)
                                                               gTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCCTTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATc                                                                                       >  2:1304578/1‑149 (MQ=255)
                                                                   ggCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTg                                                                                   >  1:2835405/1‑149 (MQ=255)
                                                                       gTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGt                                                                               >  2:778708/1‑149 (MQ=255)
                                                                        tCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTc                                                                              >  1:2983885/1‑149 (MQ=255)
                                                                         cAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCg                                                                             >  2:3533762/1‑149 (MQ=255)
                                                                           gTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGAt                                                                           <  1:3710657/149‑1 (MQ=255)
                                                                            tttCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATa                                                                          >  2:4514958/1‑149 (MQ=255)
                                                                                 aGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGcc                                                                     <  1:734845/149‑1 (MQ=255)
                                                                                 aGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGcc                                                                     >  2:1910728/1‑149 (MQ=255)
                                                                                    ttGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCt                                                                  <  2:2009762/149‑1 (MQ=255)
                                                                                    ttGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCt                                                                  >  2:2419731/1‑149 (MQ=255)
                                                                                     tGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTg                                                                 <  1:1678087/149‑1 (MQ=255)
                                                                                       gTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGa                                                               >  2:1998755/1‑149 (MQ=255)
                                                                                             aGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTATAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGc                                                         <  2:3045131/149‑1 (MQ=255)
                                                                                                    cTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAAtt                                                  <  2:230248/149‑1 (MQ=255)
                                                                                                    cTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAAtt                                                  <  1:4246517/149‑1 (MQ=255)
                                                                                                    cTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAAtt                                                  <  1:2093583/149‑1 (MQ=255)
                                                                                                    cTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAAtt                                                  <  1:4075495/149‑1 (MQ=255)
                                                                                                    cTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAAtt                                                  <  1:3660067/149‑1 (MQ=255)
                                                                                                     tGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTg                                                 <  2:362514/149‑1 (MQ=255)
                                                                                                      gTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTGa                                                >  1:907863/1‑149 (MQ=255)
                                                                                                       tCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTGAc                                               <  2:2368936/149‑1 (MQ=255)
                                                                                                         aaCGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTGACTg                                             <  2:4247449/149‑1 (MQ=255)
                                                                                                            gTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTGACTGAat                                          >  2:392582/1‑149 (MQ=255)
                                                                                                             tCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTAAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTGACTGAata                                         <  1:2012092/149‑1 (MQ=255)
                                                                                                              cacaCCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTGACTGAatat                                        <  1:2672405/149‑1 (MQ=255)
                                                                                                              cacaCCATCGGTGGCATTCACCAGCGTCCCATCGCCGTAAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTGACTGAatat                                        >  1:4550472/1‑149 (MQ=255)
                                                                                                               acacCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTGACTGAATATc                                       <  1:2282214/149‑1 (MQ=255)
                                                                                                                cacCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTGACTGAATATcc                                      <  2:3383612/149‑1 (MQ=255)
                                                                                                                cacCATCGGTGGCATTCACCAGCGTCCCATCGCCGTAAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTGACTGAATATcc                                      <  2:2750685/149‑1 (MQ=255)
                                                                                                                cacCATCGGTGGCATTCACCAGCGTCCCATCGCCGTAAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTGACTGAATATcc                                      <  1:1190517/149‑1 (MQ=255)
                                                                                                                  ccATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTGACTGAATATCCTg                                    <  1:103288/149‑1 (MQ=255)
                                                                                                                   cATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTGACTGAATATCCTGa                                   <  1:3338090/149‑1 (MQ=255)
                                                                                                                    aTCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTGACTGAATATCCTGAt                                  >  2:953736/1‑149 (MQ=255)
                                                                                                                    aTCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTGACTGAATATCCTGAt                                  <  2:892503/149‑1 (MQ=255)
                                                                                                                     tCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTGACTGAATATCCTGATc                                 <  1:1967703/149‑1 (MQ=255)
                                                                                                                     tCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTGACTGAATATCCTGATc                                 >  1:2642974/1‑149 (MQ=255)
                                                                                                                           gCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTGACTGAATATCCTGATCCGCGcc                           <  1:4034632/149‑1 (MQ=255)
                                                                                                                            cATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTGACTGAATATCCTGATCCGCGCCa                          >  2:1302072/1‑149 (MQ=255)
                                                                                                                              ttCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAAGTGACTGAATATCCTGATCCGCGCCAGt                        >  2:3742078/1‑149 (MQ=255)
                                                                                                                               tCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTGACTGAATATCCTGATCCGCGCCAGTa                       >  1:2752899/1‑149 (MQ=255)
                                                                                                                                 aCCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGACATAGCATCAATTGACTGAATATCCTGATCCGCGCCAGTaac                     >  1:691856/1‑149 (MQ=255)
                                                                                                                                  ccAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTGACTGAATATCCTGATCCGCGCCAGTaacg                    <  1:800437/149‑1 (MQ=255)
                                                                                                                                     gCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTGACTGAATATCCTGATCCGCGCCAGTaacgaac                 >  2:2854565/1‑149 (MQ=255)
                                                                                                                                     gCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTGACTGAATATCCTGATCCGCGCCAGTaacgaac                 <  1:3863749/149‑1 (MQ=255)
                                                                                                                                       gTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTGACTGAATATCCTGATCCGCGCCAGTAACGAACGt               <  1:3211262/149‑1 (MQ=255)
                                                                                                                                       gTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGGACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTGACTGAATATCCTGATCCGCGCCAGTAACGAACGt               <  1:216779/149‑1 (MQ=255)
                                                                                                                                          ccATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTGACTGAATATCCTGATCCGCGCCAGTAACGAACGTTGc            >  1:1374941/1‑149 (MQ=255)
                                                                                                                                           cATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTGACTGAATATCCTGATCCGCGCCAGTAACGAACGTTGcc           >  2:4268678/1‑149 (MQ=255)
                                                                                                                                             tCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTGACTGAATATCCTGATCCGCGCCAGTAACGAACGTTGcccc         >  2:708920/1‑149 (MQ=255)
                                                                                                                                               gCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTGACTGAATATCCTGATCCGCGCCAGTAACGAACGTTGCCCCAt       >  1:4607625/1‑149 (MQ=255)
                                                                                                                                                 cGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTGACTGAATATCCTGATCCGCGCCAGTAACGAACGTTGCCCCATCa     <  2:2600622/149‑1 (MQ=255)
                                                                                                                                                   ttAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTGACTGAATATCCTGATCCGCGCCAGTAACGAACGTTGCCCCATCAcc   <  2:2200839/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGCAGGCTAACCGCACCGGACGTATTGGTCAGATTGCCATTAACCGTCACGTTGGAAAGATAAGTCAGGCTGTCAGTTTCAAGGTTGGTATTAAGCTCACCTGTCAACGTCACACCATCGGTGGCATTCACCAGCGTCCCATCGCCGTTAAATAGTGAGGCATTAAGGGCGACGGAAGTATCCTGCCCGGTCAGGCGCAAAACCGTACCATCGCTGATGTCGATAGTGCCGCTGGAAATAGCATCAATTGACTGAATATCCTGATCCGCGCCAGTAACGAACGTTGCCCCATCACCA  >  CP009273/2335235‑2335531

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: