Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,162,618 C→A 5.6% H212N (CAT→AAT)  murB → UDP‑N‑acetylenolpyruvoylglucosamine reductase, FAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,162,6180CA 5.6% 1244.8 / 26.4 320H212N (CAT→AAT) murBUDP‑N‑acetylenolpyruvoylglucosamine reductase, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base C (146/156);  new base A (7/11);  total (153/167)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 323 total aligned reads displayed.

TTTTAAACATGAATACCAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTAAAGCATTACTGTCACAATTTCCAACAGCACCAAATTACCCCCAGGCGGATGGTTCAGT  >  CP009273/4162470‑4162766
                                                                                                                                                    |                                                                                                                                                    
  ttAAACATGAATACCAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCat                                                                                                                                                    <  1:2470732/149‑1 (MQ=255)
  ttAAACATGAATACCAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCat                                                                                                                                                    <  2:4199000/149‑1 (MQ=255)
  ttAAACATGAATACCAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCat                                                                                                                                                    >  1:4726911/1‑149 (MQ=255)
  ttAAACATGAATACCAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCat                                                                                                                                                    >  1:1761814/1‑149 (MQ=255)
     aaCATGAATACCAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATg                                                                                                                                                 <  2:4443423/149‑1 (MQ=255)
      aCATGAATACCAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATgc                                                                                                                                                <  2:872275/149‑1 (MQ=255)
      aCATGAATACCAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATgc                                                                                                                                                >  1:706211/1‑149 (MQ=255)
         tGAATACCAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGca                                                                                                                                             >  1:1403431/1‑149 (MQ=255)
         tGAATACCAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGca                                                                                                                                             >  1:1531860/1‑149 (MQ=255)
           aaTACCAGGACCGCTTCGCTATTGTATCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGcacc                                                                                                                                           >  1:4037300/1‑149 (MQ=255)
           aaTACCAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGcacc                                                                                                                                           >  2:2350326/1‑149 (MQ=255)
            aTACCAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGcacca                                                                                                                                          >  2:3161131/1‑149 (MQ=255)
            aTACCAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGcacca                                                                                                                                          >  2:1567488/1‑149 (MQ=255)
              aCCAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCaccacc                                                                                                                                        >  2:388955/1‑149 (MQ=255)
              aCCAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCaccacc                                                                                                                                        <  1:2858021/149‑1 (MQ=255)
               ccAGGACCGCTTCGCTATTGTATCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGAGCAccacca                                                                                                                                       >  1:3509007/1‑149 (MQ=255)
               ccAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCAccacca                                                                                                                                       >  1:1252454/1‑149 (MQ=255)
                cAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCaa                                                                                                                                      <  1:3432526/149‑1 (MQ=255)
                 aGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCaaa                                                                                                                                     >  2:899335/1‑149 (MQ=255)
                 aGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCaaa                                                                                                                                     >  2:778444/1‑149 (MQ=255)
                  ggACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAAc                                                                                                                                    >  1:764950/1‑149 (MQ=255)
                  ggACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAAc                                                                                                                                    >  1:4212473/1‑149 (MQ=255)
                     ccGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAAcccc                                                                                                                                 <  1:1682466/149‑1 (MQ=255)
                       gCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCt                                                                                                                               <  1:3270384/149‑1 (MQ=255)
                          tCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGAt                                                                                                                            <  1:4292728/149‑1 (MQ=255)
                          tCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGAt                                                                                                                            >  1:3706053/1‑149 (MQ=255)
                          tCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGAt                                                                                                                            <  1:3772985/149‑1 (MQ=255)
                          tCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGAt                                                                                                                            <  1:556876/149‑1 (MQ=255)
                           cGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATc                                                                                                                           <  2:3712389/149‑1 (MQ=255)
                           cGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATc                                                                                                                           >  1:882118/1‑149 (MQ=255)
                             cTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCa                                                                                                                         >  1:3124773/1‑149 (MQ=255)
                             cTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCa                                                                                                                         >  2:186259/1‑149 (MQ=255)
                                 tGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACGAAACTCCCTGATCCAAAAg                                                                                                                     <  1:2687499/149‑1 (MQ=255)
                                  gTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGt                                                                                                                    <  1:2207183/149‑1 (MQ=255)
                                  gTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTAATATGCGCTCCACCAAACTCCCTGATCCAAAAGt                                                                                                                    <  1:2352905/149‑1 (MQ=255)
                                    aGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGa                                                                                                                  >  1:1797293/1‑149 (MQ=255)
                                     gCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGaa                                                                                                                 <  1:3498556/149‑1 (MQ=255)
                                      ccGTAGGTCTGCGTCTGACAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAAt                                                                                                                <  1:3984685/149‑1 (MQ=255)
                                       cGTAGGTCTGCGTCTGCAAAAAGAGTGTCAACCTGTACTAACGTATGGTGAATTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTAATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATg                                                                                                               <  2:1098088/149‑1 (MQ=255)
                                           ggTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCaa                                                                                                           >  2:2179157/1‑149 (MQ=255)
                                            gtctgcgtctgcCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAAt                                                                                                          >  1:2448368/1‑149 (MQ=255)
                                             tctgcgtctgcCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATg                                                                                                         >  2:4640072/1‑149 (MQ=255)
                                               tgcgtctgcCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGcc                                                                                                       <  2:1165630/149‑1 (MQ=255)
                                               tgcgtctgcCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGcc                                                                                                       >  1:2026313/1‑149 (MQ=255)
                                                gcgtctgcCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCg                                                                                                      >  2:4157593/1‑149 (MQ=255)
                                                       ccAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGttt                                                                                               >  2:4586578/1‑149 (MQ=255)
                                                        cAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTttc                                                                                              <  2:764950/149‑1 (MQ=255)
                                                          aaaGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTttctt                                                                                            <  2:2758968/149‑1 (MQ=255)
                                                          aaaGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTttctt                                                                                            >  1:4703653/1‑149 (MQ=255)
                                                            agagTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCa                                                                                          >  1:4345818/1‑149 (MQ=255)
                                                            agagTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCa                                                                                          >  1:3032441/1‑149 (MQ=255)
                                                             gagTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCaa                                                                                         <  2:160233/149‑1 (MQ=255)
                                                              agTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTAATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCaaa                                                                                        <  1:4497823/149‑1 (MQ=255)
                                                                       cTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCtgtt                                                                               <  1:935580/149‑1 (MQ=255)
                                                                          tACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTa                                                                            <  1:3816634/149‑1 (MQ=255)
                                                                            cTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATc                                                                          <  2:3510084/149‑1 (MQ=255)
                                                                             tAACGTATGGTGACTTAACTCGTCTTGATCCTACAACAGTAACGCCTCAGCATGTATTTAATGCGGTGTGTCATATGCGCACCACCAATCTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCt                                                                         <  2:3616939/149‑1 (MQ=255)
                                                                               aCGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGc                                                                       <  2:4545740/149‑1 (MQ=255)
                                                                               aCGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGc                                                                       >  1:1467223/1‑149 (MQ=255)
                                                                                cGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGcc                                                                      <  2:3278334/149‑1 (MQ=255)
                                                                                    tGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGaaa                                                                  <  2:1769129/149‑1 (MQ=255)
                                                                                     ggTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAAc                                                                 <  2:4563326/149‑1 (MQ=255)
                                                                                      gTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACg                                                                >  2:3779378/1‑149 (MQ=255)
                                                                                      gTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACg                                                                >  2:2645923/1‑149 (MQ=255)
                                                                                        gACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGc                                                              <  1:2963154/149‑1 (MQ=255)
                                                                                         aCTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCt                                                             <  1:2930624/149‑1 (MQ=255)
                                                                                         aCTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTAATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCt                                                             >  1:4421144/1‑149 (MQ=255)
                                                                                          cTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTa                                                            >  1:1040129/1‑149 (MQ=255)
                                                                                           ttAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTaa                                                           >  2:2882285/1‑149 (MQ=255)
                                                                                            tAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTaaa                                                          <  2:231395/149‑1 (MQ=255)
                                                                                            tAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTaaa                                                          >  1:3533780/1‑149 (MQ=255)
                                                                                               cTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTAATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTAAAGCa                                                       >  2:806305/1‑149 (MQ=255)
                                                                                                  gTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTAAAGCATTa                                                    >  2:3094207/1‑149 (MQ=255)
                                                                                                    cTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTAAAGCATTACt                                                  <  2:4002983/149‑1 (MQ=255)
                                                                                                       gATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTAAAGCATTACTGTc                                               >  2:129839/1‑149 (MQ=255)
                                                                                                         tCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTAAAGCATTACTGTcac                                             <  2:4330320/149‑1 (MQ=255)
                                                                                                           cTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTAAAGCATTACTGTCACaa                                           >  2:2000811/1‑149 (MQ=255)
                                                                                                            tACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTAAAGCATTACTGTCACAAt                                          <  2:3911697/149‑1 (MQ=255)
                                                                                                              caacaGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTAAAGCATTACTGTCACAAttt                                        >  2:3047100/1‑149 (MQ=255)
                                                                                                                acaGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTAAATCATTACTGTCACAATTTcc                                      >  1:423998/1‑149 (MQ=255)
                                                                                                                   gTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTAAAGCATTACTGTCACAATTTCCAAc                                   >  2:2108649/1‑149 (MQ=255)
                                                                                                                    tAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTAAAGCATTACTGTCACAATTTCCAACa                                  >  2:2476201/1‑149 (MQ=255)
                                                                                                                      aCGCCACAGCACGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTAAAGCATTACTGTCACAATTTCCAACAGc                                >  1:4617873/1‑149 (MQ=255)
                                                                                                                      aCGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACTGCTAAAGCATTACTGTCACAATTTCAAACAga                                >  1:3700834/1‑148 (MQ=255)
                                                                                                                      aCGCCACAGCAAGTATTTAATGCGGTGTGTAATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTAAAGCATTACTGTCACAATTTCCAACAGc                                >  1:4145127/1‑149 (MQ=255)
                                                                                                                         ccACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTAAAGCATTACTGTCACAATTTCCAACAGCAcc                             <  2:1056351/149‑1 (MQ=255)
                                                                                                                         ccACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTAAAGCATTACTGTCACAATTTCCAACAGCAcc                             <  2:378780/149‑1 (MQ=255)
                                                                                                                          cacaGCAAGTATTTAATGCGGTGTGTAATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTAAAGCATTACTGTCACAATTTCCAACAGCACCa                            >  1:1373736/1‑149 (MQ=255)
                                                                                                                           acaGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTAAAGCATTACTGTCACAATTTCCAACAGCACCaa                           >  1:4376876/1‑149 (MQ=255)
                                                                                                                              gCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTAAAGCATTACTGTCACAATTTCCAACAGCACCAAAtt                        >  1:1556637/1‑149 (MQ=255)
                                                                                                                                  gTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTAAAGCATTACTGTCACAATTTCCAACAGCACCAAATTAccc                    >  1:690827/1‑149 (MQ=255)
                                                                                                                                   tATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTAAAGCATTACTGTCACAATTTCCAACAGCACCAAATTAcccc                   >  1:1862706/1‑149 (MQ=255)
                                                                                                                                    aTTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTAAAGCATTACTGTCACAATTTCCAACAGCACCAAATTAccccc                  >  1:1964288/1‑149 (MQ=255)
                                                                                                                                     tttAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTAAAGCATTACTGTCACAATTTCCAACAGCACCAAATTACCCCCa                 >  2:335606/1‑149 (MQ=255)
                                                                                                                                     tttAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTAAAGCATTACTGTCACAATTTCCAACAGCACCAAATTACCCCCa                 <  2:972387/149‑1 (MQ=255)
                                                                                                                                      ttAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTAAAGCATTACTGTCACAATTTCCAACAGCACCAAATTACCCCCAg                <  1:2783398/149‑1 (MQ=255)
                                                                                                                                              gtgtgtCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTAAAGCATTACTGTCACAATTTCCAACAGCACCAAATTACCCCCAGGCGGATgg        <  2:2961359/149‑1 (MQ=255)
                                                                                                                                              gtgtgtCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTAAAGCATTACTGTCACAATTTCCAACAGCACCAAATTACCCCCAGGCGGATgg        >  1:6513/1‑149 (MQ=255)
                                                                                                                                                 tgtCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTAAAGCATTACTGTCACAATTTCCAACAGCACCAAATTACCCCCAGGCGGATGGTTc     <  1:1567686/149‑1 (MQ=255)
                                                                                                                                                    cATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTAAAGCATTACTGTCACAATTTCCAACAGCACCAAATTACCCCCAGGCGGATGGTTCAGt  <  1:3584379/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTTTAAACATGAATACCAGGACCGCTTCGCTATTGTAGCCGTAGGTCTGCGTCTGCCAAAAGAGTGGCAACCTGTACTAACGTATGGTGACTTAACTCGTCTGGATCCTACAACAGTAACGCCACAGCAAGTATTTAATGCGGTGTGTCATATGCGCACCACCAAACTCCCTGATCCAAAAGTGAATGGCAATGCCGGTAGTTTCTTCAAAAACCCTGTTGTATCTGCCGAAACGGCTAAAGCATTACTGTCACAATTTCCAACAGCACCAAATTACCCCCAGGCGGATGGTTCAGT  >  CP009273/4162470‑4162766

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: