Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,009,705 G→A 100% intergenic (‑176/‑86) ysgA ← / → udp putative carboxymethylenebutenolidase/uridine phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,009,7050GA95.7% 858.0 / 20.3 256intergenic (‑176/‑86)ysgA/udpputative carboxymethylenebutenolidase/uridine phosphorylase
Reads supporting (aligned to +/- strand):  ref base G (3/8);  new base A (122/122);  total (125/131)
Fisher's exact test for biased strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.58e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 261 total aligned reads displayed.

GCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTGATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGTATATGTCCAAGTCTGATGTTTTTCATCTCGGCCTCACTAAAAACGATTTACAAGGGGCTACG  >  CP009273/4009559‑4009850
                                                                                                                                                  |                                                                                                                                                 
gcgcAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGtaatg                                                                                                                                                <  2:4001663/150‑1 (MQ=255)
gcgcAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGtaatg                                                                                                                                                <  2:2683089/150‑1 (MQ=255)
gcgcAAGGTTAACTATAGTCAGCATGCAACAACTCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTCTGCAACGCATTTGCGTCATGGtaatg                                                                                                                                                <  2:1533314/150‑1 (MQ=255)
 cgcAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGtaatga                                                                                                                                               <  1:1526927/150‑1 (MQ=255)
 cgcAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGtaatga                                                                                                                                               >  1:2460870/1‑150 (MQ=255)
  gcAATGTTAACTAGAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAg                                                                                                                                              <  2:2192212/150‑1 (MQ=255)
     aTGTGAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTAt                                                                                                                                           <  2:260351/150‑1 (MQ=255)
        ttAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCAc                                                                                                                                        <  1:4510025/150‑1 (MQ=255)
        ttAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCAc                                                                                                                                        >  2:1444600/1‑150 (MQ=255)
             tataGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGaaaa                                                                                                                                   <  1:822646/150‑1 (MQ=255)
              ataGTCAGCATGCAACAAATCACCTTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGaaaaa                                                                                                                                  <  2:107550/150‑1 (MQ=255)
                aGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAAt                                                                                                                                <  1:591211/150‑1 (MQ=255)
                   cAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGtt                                                                                                                             <  2:4131292/150‑1 (MQ=255)
                    aGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTa                                                                                                                            >  1:3460761/1‑150 (MQ=255)
                      cATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTaaa                                                                                                                          <  1:4157129/150‑1 (MQ=255)
                      cATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTaaa                                                                                                                          <  1:2246895/150‑1 (MQ=255)
                         gCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAAccc                                                                                                                       >  1:594276/1‑150 (MQ=255)
                           aacaaATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCtt                                                                                                                     >  1:1499806/1‑150 (MQ=255)
                            acaaATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTTTGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTc                                                                                                                    <  1:1141951/150‑1 (MQ=255)
                            acaaATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTc                                                                                                                    >  2:1168918/1‑150 (MQ=255)
                            acaaATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTc                                                                                                                    <  2:3262574/150‑1 (MQ=255)
                            acaaATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTc                                                                                                                    <  1:4211662/150‑1 (MQ=255)
                             caaATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCg                                                                                                                   <  1:74726/150‑1 (MQ=255)
                              aaaTCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCgg                                                                                                                  <  1:2326254/150‑1 (MQ=255)
                              aaaTCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCgg                                                                                                                  <  2:1375446/150‑1 (MQ=255)
                              aaaTCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCgg                                                                                                                  >  2:3466776/1‑150 (MQ=255)
                              aaaTCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCgg                                                                                                                  <  2:3868646/150‑1 (MQ=255)
                                 tCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTaa                                                                                                               <  1:2347037/150‑1 (MQ=255)
                                 tCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTaa                                                                                                               <  2:1921745/150‑1 (MQ=255)
                                  cacaTTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTaaa                                                                                                              >  2:910661/1‑150 (MQ=255)
                                    caTTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAgt                                                                                                            >  2:3685825/1‑150 (MQ=255)
                                     aTTGCCTGAATCGTCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAgtg                                                                                                           <  2:2795679/150‑1 (MQ=255)
                                       tGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTc                                                                                                         >  1:909137/1‑150 (MQ=255)
                                       tGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGAGTCATGGTACTGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTc                                                                                                         >  2:486450/1‑150 (MQ=255)
                                        gCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCt                                                                                                        <  1:2642648/150‑1 (MQ=255)
                                        gCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCt                                                                                                        >  2:3522915/1‑150 (MQ=255)
                                          cTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCttt                                                                                                      <  1:2415420/150‑1 (MQ=255)
                                             aaTCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTg                                                                                                   <  1:4407569/150‑1 (MQ=255)
                                             aaTCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTg                                                                                                   >  1:1210571/1‑150 (MQ=255)
                                                cGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGctt                                                                                                <  2:602721/150‑1 (MQ=255)
                                                 ggCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGcttc                                                                                               >  2:4480987/1‑150 (MQ=255)
                                                  gCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGcttct                                                                                              >  2:1177328/1‑150 (MQ=255)
                                                  gCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGcttct                                                                                              >  1:439032/1‑150 (MQ=255)
                                                     cATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTGATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTtcttct                                                                                           >  2:4196145/1‑150 (MQ=255)
                                                     cATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTtcttct                                                                                           >  1:4603771/1‑150 (MQ=255)
                                                     cATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTtcttct                                                                                           >  2:1243527/1‑150 (MQ=255)
                                                     cATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTtcttct                                                                                           <  2:566148/150‑1 (MQ=255)
                                                     cATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTtcttct                                                                                           <  1:1288871/150‑1 (MQ=255)
                                                      aTCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTg                                                                                          >  2:1952089/1‑150 (MQ=255)
                                                       tCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGa                                                                                         >  2:3391390/1‑150 (MQ=255)
                                                       tCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGa                                                                                         >  2:3716465/1‑150 (MQ=255)
                                                       tCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGa                                                                                         >  1:967090/1‑150 (MQ=255)
                                                        cTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGAc                                                                                        >  1:3596567/1‑150 (MQ=255)
                                                         ttttATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACt                                                                                       <  1:3654083/150‑1 (MQ=255)
                                                          tttATGCAGGCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTa                                                                                      <  1:3829988/150‑1 (MQ=255)
                                                           ttATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTaa                                                                                     >  1:1928612/1‑150 (MQ=255)
                                                           ttATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTaa                                                                                     <  2:3357732/150‑1 (MQ=255)
                                                               gCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCg                                                                                 <  2:4342842/150‑1 (MQ=255)
                                                                   tCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTc                                                                             <  2:2082452/150‑1 (MQ=255)
                                                                   tCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTc                                                                             <  2:4087922/150‑1 (MQ=255)
                                                                      tGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTcaca                                                                          >  1:1071917/1‑150 (MQ=255)
                                                                      tGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTcaca                                                                          >  2:1802764/1‑150 (MQ=255)
                                                                       gCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACag                                                                         <  1:4504648/150‑1 (MQ=255)
                                                                          gAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAgagg                                                                      <  1:767449/150‑1 (MQ=255)
                                                                              gAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGtt                                                                  <  1:2764698/150‑1 (MQ=255)
                                                                               aaGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTg                                                                 <  1:1789789/150‑1 (MQ=255)
                                                                                   gTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGtata                                                             >  1:629167/1‑150 (MQ=255)
                                                                                   gTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGtata                                                             >  2:4608519/1‑150 (MQ=255)
                                                                                        ttATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGTATATGTcc                                                        >  2:1961073/1‑150 (MQ=255)
                                                                                         tATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGTATATGTCCa                                                       <  1:394874/150‑1 (MQ=255)
                                                                                         tATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGTATATGTCCa                                                       >  1:2918487/1‑150 (MQ=255)
                                                                                          aTGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGTATATGTCCaa                                                      >  1:1668074/1‑150 (MQ=255)
                                                                                             tgATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGTATATGTCCAAGTc                                                   >  1:209023/1‑150 (MQ=255)
                                                                                               aTTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGTATATGTCCAAGTCTg                                                 <  2:3307212/150‑1 (MQ=255)
                                                                                               aTTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGTATATGTCCAAGTCTg                                                 >  2:1645113/1‑150 (MQ=255)
                                                                                                     aTCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGTATATGTCCAAGTCTGATGttt                                           >  1:3113929/1‑150 (MQ=255)
                                                                                                        aCTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGTATATGTCCAAGTCTGATGTTTTTc                                        <  2:3021318/150‑1 (MQ=255)
                                                                                                        aCTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGTATATGTCCAAGTCTGATGTTTTTc                                        <  2:706316/150‑1 (MQ=255)
                                                                                                                 gggTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGTATATGTCCAAGTCTGATGTTTTTCATCTCGGcc                               >  2:281760/1‑150 (MQ=255)
                                                                                                                  ggTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGTATATGTCCAAGTCTGATGTTTTTCATCTCGGCCt                              <  1:1471699/150‑1 (MQ=255)
                                                                                                                  ggTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGTATATGTCCAAGTCTGATGTTTTTCATCTCGGCCt                              <  1:1471304/150‑1 (MQ=255)
                                                                                                                  ggTAAATTTATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGTATATGTCCAAGTCTGATGTTTTTCATCTCGGCCt                              >  1:3217867/1‑150 (MQ=255)
                                                                                                                        tttATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGTATATGTCCAAGTCTGATGTTTTTCATCTCGGCCTCACTaa                        <  1:2016140/150‑1 (MQ=255)
                                                                                                                         ttATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGTATATGTCCAAGTCTGATGTTTTTCATCTCGGCCTCACTaaa                       >  1:1006488/1‑150 (MQ=255)
                                                                                                                         ttATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGTATATGTCCAAGTCTGATGTTTTTCATCTCGGCCTCACTaaa                       >  1:512698/1‑150 (MQ=255)
                                                                                                                         ttATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGTATATGTCCAAGTCTGATGTTTTTCATCTCGGCCTCACTaaa                       >  1:2392875/1‑150 (MQ=255)
                                                                                                                         ttATGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATACACATAGGAGTTGGATATGTCCAAGTCTGATGTTTTTCATCTCGGCCTCACTaaa                       >  1:3324093/1‑150 (MQ=255)
                                                                                                                            tGCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGTATATGTCCAAGTCTGATGTTTTTCATCTCGGCCTCACTAAAAAc                    >  2:500201/1‑150 (MQ=255)
                                                                                                                             gCAACGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGTATATGTCCAAGTCTGATGTTTTTCATCTCGGCCTCACTAAAAACg                   >  2:3438110/1‑150 (MQ=255)
                                                                                                                               aaCGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGTATATGTCCAAGTCTGATGTTTTTCATCTCGGCCTCACTAAAAACGAt                 >  1:47454/1‑150 (MQ=255)
                                                                                                                               aaCGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGTATATGTCCAAGTCTGATGTTTTTCATCTCGGCCTCACTAAAAACGAt                 <  1:3710004/150‑1 (MQ=255)
                                                                                                                                 cGCATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGTATATGTCCAAGTCTGATGTTTTTCATCTCGGCCTCACTAAAAACGAttt               >  1:3034307/1‑150 (MQ=255)
                                                                                                                                   cATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGTATATGTCCAAGTCTGATGTTTTTCATCTCGGCCTCACTAAAAACGATTTAc             <  2:4252439/150‑1 (MQ=255)
                                                                                                                                   cATTTGCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGTATATGTCCAAGTCTGATGTTTTTCATCTCGGCCTCACTAAAAACGATTTAc             <  2:2292562/150‑1 (MQ=255)
                                                                                                                                        gCGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGTATATGTCCAAGTCTGATGTTTTTCATCTCGGCCTCACTAAAAACGATTTACAAggg        <  1:2582025/150‑1 (MQ=255)
                                                                                                                                         cGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGTATATGTCCAAGTCTGATGTTTTTCATCTCGGCCTCACTAAAAACGATTTACAAgggg       <  2:4410048/150‑1 (MQ=255)
                                                                                                                                         cGTCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGTATATGTCCAAGTCTGATGTTTTTCATCTCGGCCTCACTAAAAACGATTTACAAgggg       >  1:428426/1‑150 (MQ=255)
                                                                                                                                           tCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGTATATGTCCAAGTCTGATGTTTTTCATCTCGGCCTCACTAAAAACGATTTACAAGGGGCt     >  1:3489752/1‑150 (MQ=255)
                                                                                                                                           tCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGTATATGTCCAAGTCTGATGTTTTTCATCTCGGCCTCACTAAAAACGATTTACAAGGGGCt     <  2:516236/150‑1 (MQ=255)
                                                                                                                                           tCATGGTAATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGTATATGTCCAAGTCTGATGTTTTTCATCTCGGCCTCACTAAAAACGATTTACAAGGGGCt     <  2:1750788/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTCATCTTTTATGCAGTCCTGCAGAATGAAGGGTGATTTATGTGATTTGCATCACTTTTGGTGGGTAAATTTATGCAACGCATTTGCGTCATGGTGATGAGTATCACGAAAAAATGTTAAACCCTTCGGTAAAGTGTCTTTTTGCTTCTTCTGACTAAACCGATTCACAGAGGAGTTGTATATGTCCAAGTCTGATGTTTTTCATCTCGGCCTCACTAAAAACGATTTACAAGGGGCTACG  >  CP009273/4009559‑4009850

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: