Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,096,555 G→T 6.1% S53S (TCG→TCT fieF → ferrous iron and zinc transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,096,5550GT 6.1% 1078.5 / 26.8 286S53S (TCG→TCTfieFferrous iron and zinc transporter
Reads supporting (aligned to +/- strand):  ref base G (129/139);  new base T (6/12);  total (135/151)
Fisher's exact test for biased strand distribution p-value = 3.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 296 total aligned reads displayed.

CTTATGGACGGCTGGTCAGTCGGGCGGCGATTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTATTCCTGTTTTTGACGGGTATTC  >  CP009273/4096407‑4096703
                                                                                                                                                    |                                                                                                                                                    
 ttATGGACGGCTGGTCAGTCGGGCGGCGATTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGtt                                                                                                                                                    <  2:3369785/150‑1 (MQ=255)
        cGGCTGGTCCGTCGGGAGGCGATTGCTGCGACGGCGATGGCATCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCTGCGAGTCGTTGACTAAt                                                                                                                                             >  2:3577788/1‑150 (MQ=255)
         ggCTGGTCAGTCGGGCGGCGATTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAAtt                                                                                                                                            <  2:2408518/150‑1 (MQ=255)
         ggCTGGTCAGTCGGGCGGCGATTGCTGCGACGGAGCTGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAAtt                                                                                                                                            >  2:936157/1‑150 (MQ=255)
          gCTGGTCAGTCGGGCGGCGATTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAAttt                                                                                                                                           >  1:3157953/1‑150 (MQ=255)
              gTCAGTCGGGCGGCGATTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATtg                                                                                                                                       <  2:726440/150‑1 (MQ=255)
               tCAGTCGGGCGGCGATTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATtgg                                                                                                                                      <  1:229172/150‑1 (MQ=255)
                 aGTCGGGCGGCGATTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCTTTGACGAATTTATtggtg                                                                                                                                    >  2:782506/1‑150 (MQ=255)
                 aGTCGGGCGGCGATTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCTTTGACGAATTTATtggtg                                                                                                                                    >  2:2894533/1‑150 (MQ=255)
                 aGTCGGGCGGCGATTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATtggtg                                                                                                                                    >  2:847811/1‑150 (MQ=255)
                     gggCGGCGATTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGc                                                                                                                                <  1:267395/150‑1 (MQ=255)
                          gcgATTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCTTTGACGAATTTATTGGTGGTGCGatat                                                                                                                           <  2:655319/150‑1 (MQ=255)
                           cgATTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATAtt                                                                                                                          <  1:1513434/150‑1 (MQ=255)
                            gATTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTc                                                                                                                         >  2:3363181/1‑150 (MQ=255)
                            gATTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTc                                                                                                                         >  1:3362937/1‑150 (MQ=255)
                             aTTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTcc                                                                                                                        <  1:1581490/150‑1 (MQ=255)
                               tgctgcGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATCTTCCCt                                                                                                                      >  2:4090144/1‑150 (MQ=255)
                               tgctgcGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCt                                                                                                                      >  2:2957824/1‑150 (MQ=255)
                               tgctgcGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCt                                                                                                                      <  1:2938541/150‑1 (MQ=255)
                                 ctgcGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGc                                                                                                                    <  1:512101/150‑1 (MQ=255)
                                   gcGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCaa                                                                                                                  <  2:1518671/150‑1 (MQ=255)
                                      aCGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCt                                                                                                               >  2:3603217/1‑150 (MQ=255)
                                       cGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTg                                                                                                              <  2:1817041/150‑1 (MQ=255)
                                        ggCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCTTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGc                                                                                                             >  1:1229760/1‑150 (MQ=255)
                                        ggCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGc                                                                                                             >  2:1185379/1‑150 (MQ=255)
                                            aTGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCcgac                                                                                                         <  1:2247916/150‑1 (MQ=255)
                                            aTGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCcgac                                                                                                         >  1:4400392/1‑150 (MQ=255)
                                              ggCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCcgacga                                                                                                       <  2:2107551/150‑1 (MQ=255)
                                               gCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGAt                                                                                                      <  2:594791/150‑1 (MQ=255)
                                               gCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGAt                                                                                                      >  1:2065227/1‑150 (MQ=255)
                                                cTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCTTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATa                                                                                                     <  1:2034663/150‑1 (MQ=255)
                                                  tCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAAt                                                                                                   <  2:3046326/150‑1 (MQ=255)
                                                    gctgctATTGCCGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCa                                                                                                 <  1:3241082/150‑1 (MQ=255)
                                                       gctATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTCTTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTc                                                                                              >  2:2576252/1‑150 (MQ=255)
                                                        ctATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCTGCGCGTCGTTGACGAATTTATTGGTGGTGCTATATTCCCTGCAACCTGCCGACGATAATCACTCg                                                                                             >  1:4197686/1‑150 (MQ=255)
                                                         tATTGCTGCTTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGTGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGt                                                                                            <  1:194414/150‑1 (MQ=255)
                                                         tATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGt                                                                                            >  1:3770222/1‑150 (MQ=255)
                                                          aTTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGtt                                                                                           >  1:1851999/1‑150 (MQ=255)
                                                           ttGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGttt                                                                                          <  1:578125/150‑1 (MQ=255)
                                                             gCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTgg                                                                                        >  2:1367410/1‑150 (MQ=255)
                                                              cTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGt                                                                                       >  1:328732/1‑150 (MQ=255)
                                                               tGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTc                                                                                      >  2:1236062/1‑150 (MQ=255)
                                                               tGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTc                                                                                      >  2:2323799/1‑150 (MQ=255)
                                                                gATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCa                                                                                     <  2:4464956/150‑1 (MQ=255)
                                                                 aTTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCAc                                                                                    >  2:476092/1‑150 (MQ=255)
                                                                 aTTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCAc                                                                                    >  1:1987165/1‑150 (MQ=255)
                                                                     aaaTTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTa                                                                                <  1:3546124/150‑1 (MQ=255)
                                                                         ttttGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGc                                                                            >  1:3027013/1‑150 (MQ=255)
                                                                         ttttGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGc                                                                            >  1:2617047/1‑150 (MQ=255)
                                                                          tttGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCt                                                                           <  2:3639836/150‑1 (MQ=255)
                                                                          tttGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATAGTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCt                                                                           <  2:365818/150‑1 (MQ=255)
                                                                            tGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGa                                                                         >  1:3975270/1‑150 (MQ=255)
                                                                            tGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGa                                                                         >  1:3974081/1‑150 (MQ=255)
                                                                            tGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGa                                                                         >  1:2209708/1‑150 (MQ=255)
                                                                             gCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAg                                                                        <  1:1347015/150‑1 (MQ=255)
                                                                              cATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGt                                                                       >  1:1570341/1‑150 (MQ=255)
                                                                               aTGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTc                                                                      >  1:2098175/1‑150 (MQ=255)
                                                                               aTGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTc                                                                      >  1:4123851/1‑150 (MQ=255)
                                                                               aTGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTc                                                                      >  1:2098301/1‑150 (MQ=255)
                                                                                tggtggTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTcc                                                                     <  1:684741/150‑1 (MQ=255)
                                                                                   tggtATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTc                                                                  <  1:852296/150‑1 (MQ=255)
                                                                                   tggtATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTc                                                                  >  2:1450628/1‑150 (MQ=255)
                                                                                   tggtATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCCGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCTTTTGGTCACGGTAAAGCTGAGTCCCTc                                                                  <  1:4439879/150‑1 (MQ=255)
                                                                                    ggtATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTcg                                                                 <  2:3138754/150‑1 (MQ=255)
                                                                                    ggtATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTcg                                                                 <  2:3138702/150‑1 (MQ=255)
                                                                                     gtATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTcgc                                                                <  2:1269342/150‑1 (MQ=255)
                                                                                     gtATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTcgc                                                                <  2:1269281/150‑1 (MQ=255)
                                                                                     gtATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTcgc                                                                <  2:3734935/150‑1 (MQ=255)
                                                                                        taCCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCTTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCgcggc                                                             <  2:4085146/150‑1 (MQ=255)
                                                                                            gggTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTg                                                         >  1:4495596/1‑150 (MQ=255)
                                                                                             ggTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTgg                                                        <  1:1412990/150‑1 (MQ=255)
                                                                                               tCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGgcg                                                      <  1:1731559/150‑1 (MQ=255)
                                                                                                   tGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCaaa                                                  <  1:3382045/150‑1 (MQ=255)
                                                                                                      gTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTa                                               <  1:847811/150‑1 (MQ=255)
                                                                                                       tATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTAt                                              <  2:2799876/150‑1 (MQ=255)
                                                                                                         ttCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGt                                            <  1:3825442/150‑1 (MQ=255)
                                                                                                           ctcGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGttt                                          <  1:3709018/150‑1 (MQ=255)
                                                                                                                 gcgcTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTcc                                    >  1:2831900/1‑150 (MQ=255)
                                                                                                                  cgcTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTTCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCg                                   >  2:2802189/1‑150 (MQ=255)
                                                                                                                   gcTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCgg                                  <  2:2555744/150‑1 (MQ=255)
                                                                                                                      ggtggATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTc                               <  1:4322244/150‑1 (MQ=255)
                                                                                                                       gtggATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCg                              >  2:2981032/1‑150 (MQ=255)
                                                                                                                        tggATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCgg                             >  1:1100410/1‑150 (MQ=255)
                                                                                                                         ggATTCGCTGGTGGATATCGGCGCGTCTTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGc                            <  2:224791/150‑1 (MQ=255)
                                                                                                                          gATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCa                           >  1:1051338/1‑150 (MQ=255)
                                                                                                                          gATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCAGCa                           <  2:4222756/150‑1 (MQ=255)
                                                                                                                               gCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTAtt                      <  2:1854626/150‑1 (MQ=255)
                                                                                                                                 tggtggATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTATTcc                    >  1:340927/1‑150 (MQ=255)
                                                                                                                                    tggATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTATTCCTGt                 <  1:3949422/150‑1 (MQ=255)
                                                                                                                                    tggATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTATTCCTGt                 <  1:932161/150‑1 (MQ=255)
                                                                                                                                      gATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTATTCCTGttt               >  1:635921/1‑150 (MQ=255)
                                                                                                                                       atatCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTATTCCTGtttt              <  1:1100184/150‑1 (MQ=255)
                                                                                                                                         atCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTATTCATGTTTTTg            >  2:3016179/1‑150 (MQ=255)
                                                                                                                                            ggCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTATTCCTGTTTTTGACg         >  1:580427/1‑150 (MQ=255)
                                                                                                                                            ggCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTATTCCTGTTTTTGACg         <  1:1236062/150‑1 (MQ=255)
                                                                                                                                            ggCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTAGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTATTCCTGTTTTTGACg         <  1:3459175/150‑1 (MQ=255)
                                                                                                                                              cgcgTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTATTCCTGTTTTTGACggg       <  1:3042440/150‑1 (MQ=255)
                                                                                                                                                  tcgtTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTATTCCTGTTTTTGACGGGTAtt   <  2:707499/150‑1 (MQ=255)
                                                                                                                                                   cgtTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTATTCCTGTTTTTGACGGGTATTc  <  2:3891355/150‑1 (MQ=255)
                                                                                                                                                   cgtTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTATTCCTGTTTTTGACGGGTATTc  >  2:1944554/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTTATGGACGGCTGGTCAGTCGGGCGGCGATTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTATTCCTGTTTTTGACGGGTATTC  >  CP009273/4096407‑4096703

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: