Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,096,555 G→T 10.7% 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,5550GT10.7% 1191.6 / 66.7 329S53S (TCG→TCTfieFferrous iron and zinc transporter
Reads supporting (aligned to +/- strand):  ref base G (141/151);  new base T (14/23);  total (155/174)
Fisher's exact test for biased strand distribution p-value = 2.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 343 total aligned reads displayed.

TCTTATGGACGGCTGGTCAGTCGGGCGGCGATTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTATTCCTGTTTTTGACGGGTATTC  >  CP009273/4096406‑4096703
                                                                                                                                                     |                                                                                                                                                    
tcTTATGGACGGCTGGTCAGTCGGGCGGCGATTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGcgtcg                                                                                                                                                      >  1:1578371/1‑150 (MQ=255)
 cTTATGGACGGCTGGTCAGTCGGGCGGCGATTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGcgtcgt                                                                                                                                                     >  2:4598871/1‑150 (MQ=255)
  ttATGGACGGCTGGTCAGTCGGGCGGCGATTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCttt                                                                                                                                                    <  1:709055/150‑1 (MQ=255)
   tATGGACGGCTGGTCAGTCGGGCGGCGATTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTg                                                                                                                                                   <  1:4373389/150‑1 (MQ=255)
   tATGGACGGCTGGTCAGTCGGGCGGCGATTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTg                                                                                                                                                   <  1:4979465/150‑1 (MQ=255)
    aTGGACGGCTGGTCAGTCGGGCGGCGATTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGa                                                                                                                                                  <  2:1318466/150‑1 (MQ=255)
         cGGCTGGTCAGTCGGGCGGCGATTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAAt                                                                                                                                             <  2:1825289/150‑1 (MQ=255)
          ggCTGGTCAGTCGGGCGGCGATTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAAtt                                                                                                                                            >  1:2096507/1‑150 (MQ=255)
            cTGGTCTGTCGGTCGGCGATTGCTGCGACGGAGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTTAGTATTCTCGCCGCGCTGGTGGATTCGCTGTTGGATATCGGCGCGTCGTTGACGAATTTa                                                                                                                                          <  2:1700654/150‑1 (MQ=255)
            cTGGTCAGTCGGGCGGCGATTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTa                                                                                                                                          <  1:2820690/150‑1 (MQ=255)
             tGGTCAGTCGGGCGGCGATTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTAt                                                                                                                                         >  2:43170/1‑150 (MQ=255)
                 cAGTCGGGCGGCGATTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATtggt                                                                                                                                     >  2:4498425/1‑150 (MQ=255)
                    tCGGGCGGCGATTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTTGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTggtggt                                                                                                                                  <  1:3636332/150‑1 (MQ=255)
                     cGGGCGGCGATTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGgtggtg                                                                                                                                 <  1:3230114/150‑1 (MQ=255)
                       ggcggcGATTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCg                                                                                                                               >  2:825236/1‑150 (MQ=255)
                       ggcggcGATTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCg                                                                                                                               <  1:2224531/150‑1 (MQ=255)
                             gATTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTc                                                                                                                         >  2:928580/1‑150 (MQ=255)
                                tgctgcGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCTTTGACGAATTTATTGGTGGTGCGATATTCCCt                                                                                                                      >  1:4241114/1‑150 (MQ=255)
                                tgctgcGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCt                                                                                                                      >  1:2010329/1‑150 (MQ=255)
                                tgctgcGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCt                                                                                                                      <  1:1684112/150‑1 (MQ=255)
                                    gcGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCaa                                                                                                                  <  1:3118041/150‑1 (MQ=255)
                                       aCGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCt                                                                                                               <  1:4512005/150‑1 (MQ=255)
                                         ggCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGc                                                                                                             >  1:5014325/1‑150 (MQ=255)
                                           cGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCGGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCACCCTGCcg                                                                                                           >  2:1135462/1‑150 (MQ=255)
                                            gATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCcga                                                                                                          >  2:2819750/1‑150 (MQ=255)
                                            gATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCcga                                                                                                          >  2:4357114/1‑150 (MQ=255)
                                            gATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCcga                                                                                                          >  1:4869333/1‑150 (MQ=255)
                                            gATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCcga                                                                                                          <  1:3922896/150‑1 (MQ=255)
                                             aTGGCTTCGCTTCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCcgac                                                                                                         <  2:2316479/150‑1 (MQ=255)
                                              tGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCcgacg                                                                                                        >  1:196688/1‑150 (MQ=255)
                                               ggCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCcgacga                                                                                                       >  1:1087644/1‑150 (MQ=255)
                                                gCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGAt                                                                                                      >  1:238449/1‑150 (MQ=255)
                                                 cTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATa                                                                                                     >  1:4459542/1‑150 (MQ=255)
                                                  ttCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATaa                                                                                                    >  1:4059674/1‑150 (MQ=255)
                                                     gctgctATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCa                                                                                                 >  2:1001970/1‑150 (MQ=255)
                                                      ctgctATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCAc                                                                                                >  2:2459812/1‑150 (MQ=255)
                                                      ctgctATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGGCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCAc                                                                                                <  1:25635/150‑1 (MQ=255)
                                                       tgctATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACt                                                                                               >  2:406930/1‑150 (MQ=255)
                                                         ctATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCTTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCg                                                                                             <  1:1564060/150‑1 (MQ=255)
                                                            ttGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGttt                                                                                          <  2:58778/150‑1 (MQ=255)
                                                              gCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTgg                                                                                        >  2:2536698/1‑150 (MQ=255)
                                                               cTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGt                                                                                       >  1:1363890/1‑150 (MQ=255)
                                                                tGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTc                                                                                      <  2:2945516/150‑1 (MQ=255)
                                                                 gATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCa                                                                                     <  2:3331892/150‑1 (MQ=255)
                                                                  aTTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCAc                                                                                    >  1:620061/1‑150 (MQ=255)
                                                                  aTTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCAc                                                                                    <  2:3125858/150‑1 (MQ=255)
                                                                    tAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACgg                                                                                  >  2:2152055/1‑150 (MQ=255)
                                                                      aaaTTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTa                                                                                <  1:3503050/150‑1 (MQ=255)
                                                                           tttGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCt                                                                           >  2:5016545/1‑150 (MQ=255)
                                                                             tGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGa                                                                         >  2:4111179/1‑150 (MQ=255)
                                                                             tGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGa                                                                         <  2:1396034/150‑1 (MQ=255)
                                                                             tGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGa                                                                         >  2:4706959/1‑150 (MQ=255)
                                                                             tGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGa                                                                         <  2:3958932/150‑1 (MQ=255)
                                                                             tGCATGGTGGTATACCGGGTCGGTGAGTAGTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCTTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGa                                                                         <  2:2006266/150‑1 (MQ=255)
                                                                                aTGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTc                                                                      <  2:2787428/150‑1 (MQ=255)
                                                                                  ggtggtATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTccc                                                                    <  2:4983103/150‑1 (MQ=255)
                                                                                  ggtggtATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTccc                                                                    >  1:2337200/1‑150 (MQ=255)
                                                                                       tataCCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTcgcg                                                               <  2:4898680/150‑1 (MQ=255)
                                                                                            cGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCTGTGACGAATTTCTTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTCCAGCTGAGTCCCTCGCGGCGCt                                                          >  2:641561/1‑150 (MQ=255)
                                                                                               gTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGgc                                                       <  1:377853/150‑1 (MQ=255)
                                                                                                tCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCTTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGgcg                                                      <  1:1450763/150‑1 (MQ=255)
                                                                                                tCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGgcg                                                      <  1:527408/150‑1 (MQ=255)
                                                                                                  ggTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCa                                                    >  1:3518425/1‑150 (MQ=255)
                                                                                                     gAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAg                                                 <  2:3641750/150‑1 (MQ=255)
                                                                                                     gAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAg                                                 >  2:2744907/1‑150 (MQ=255)
                                                                                                     gAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAg                                                 <  2:1926263/150‑1 (MQ=255)
                                                                                                        tATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTAt                                              >  1:2796982/1‑150 (MQ=255)
                                                                                                         aTTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGCCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATg                                             <  1:535578/150‑1 (MQ=255)
                                                                                                          ttCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGt                                            <  2:1546701/150‑1 (MQ=255)
                                                                                                          ttCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGt                                            <  2:3010066/150‑1 (MQ=255)
                                                                                                           tctcGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGtt                                           >  1:219319/1‑150 (MQ=255)
                                                                                                            ctcGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGACGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGttt                                          <  2:3345493/150‑1 (MQ=255)
                                                                                                             tcGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCTCGTCTTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTa                                         <  1:3336759/150‑1 (MQ=255)
                                                                                                             tcGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCTTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTa                                         <  2:3768038/150‑1 (MQ=255)
                                                                                                               gccgcGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTAtc                                       >  2:4607920/1‑150 (MQ=255)
                                                                                                               gccgcGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTAtc                                       >  1:2515513/1‑150 (MQ=255)
                                                                                                                    gcTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCgg                                  >  2:3725549/1‑150 (MQ=255)
                                                                                                                       ggtggATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTc                               <  1:1056005/150‑1 (MQ=255)
                                                                                                                       ggtggATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTc                               <  2:3522706/150‑1 (MQ=255)
                                                                                                                        gtggATTCGCTGGTGGATATCGGCGCGTCTTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCg                              <  2:3610904/150‑1 (MQ=255)
                                                                                                                          ggATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGc                            >  1:1753773/1‑150 (MQ=255)
                                                                                                                             ttCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACt                         <  1:3275309/150‑1 (MQ=255)
                                                                                                                              tCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTa                        <  1:3849824/150‑1 (MQ=255)
                                                                                                                              tCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTa                        >  1:4762494/1‑150 (MQ=255)
                                                                                                                                gCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTAtt                      <  1:4681432/150‑1 (MQ=255)
                                                                                                                                gCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTAtt                      <  1:2992777/150‑1 (MQ=255)
                                                                                                                                  tggtggATATCGGCGCGTCTTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTATTcc                    <  2:4042769/150‑1 (MQ=255)
                                                                                                                                  tggtggATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTATTcc                    >  2:4275314/1‑150 (MQ=255)
                                                                                                                                    gtggATATCGGCGCGTCTTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTATTCCTg                  <  1:1014802/150‑1 (MQ=255)
                                                                                                                                       gATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTATTCCTGttt               <  2:5123199/150‑1 (MQ=255)
                                                                                                                                        atatCGGCGCGTCTTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTATTCCTGtttt              >  1:1918214/1‑150 (MQ=255)
                                                                                                                                        atatCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTATTCCTGtttt              <  2:3429996/150‑1 (MQ=255)
                                                                                                                                        atatCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTATTCCTGtttt              <  1:1676812/150‑1 (MQ=255)
                                                                                                                                        atatCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTATTCCTGtttt              >  1:787755/1‑150 (MQ=255)
                                                                                                                                          atCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTATTCCTGTTTTTg            <  2:3582650/150‑1 (MQ=255)
                                                                                                                                             ggCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTATTCCTGTTTTTGACg         >  1:184653/1‑150 (MQ=255)
                                                                                                                                               cgcgTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTATTCCTGTTTTTGACggg       >  2:1869694/1‑150 (MQ=255)
                                                                                                                                                 cgtcgtTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTATTCCTGTTTTTGACGGGTa     >  1:91746/1‑150 (MQ=255)
                                                                                                                                                   tcgtTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTATTCCTGTTTTTGACGGGTAtt   >  2:4244417/1‑150 (MQ=255)
                                                                                                                                                    cgtTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTATTCCTGTTTTTGACGGGTATTc  >  2:1982947/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCTTATGGACGGCTGGTCAGTCGGGCGGCGATTGCTGCGACGGCGATGGCTTCGCTGCTATTGCTGATTAAAATTTTTGCATGGTGGTATACCGGGTCGGTGAGTATTCTCGCCGCGCTGGTGGATTCGCTGGTGGATATCGGCGCGTCGTTGACGAATTTATTGGTGGTGCGATATTCCCTGCAACCTGCCGACGATAATCACTCGTTTGGTCACGGTAAAGCTGAGTCCCTCGCGGCGCTGGCGCAAAGTATGTTTATCTCCGGTTCGGCACTATTCCTGTTTTTGACGGGTATTC  >  CP009273/4096406‑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: