Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 166,662 C→A 100% S133* (TCG→TAG)  fhuC → iron‑hydroxamate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273166,6620CA99.2% 773.5 / ‑1.9 235S133* (TCG→TAG) fhuCiron‑hydroxamate transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (111/122);  total (112/123)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 239 total aligned reads displayed.

CGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGATATCGCC  >  CP009273/166517‑166807
                                                                                                                                                 |                                                                                                                                                 
 gCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTg                                                                                                                                               >  2:4033071/1‑149 (MQ=255)
    aaGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGtt                                                                                                                                            <  1:3753820/149‑1 (MQ=255)
      gTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTgg                                                                                                                                          <  2:1371955/149‑1 (MQ=255)
       tGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGc                                                                                                                                         <  1:1079488/149‑1 (MQ=255)
         gCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTCTCTAGCTGGTTGGCtt                                                                                                                                       >  2:3792849/1‑149 (MQ=255)
         gCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCtt                                                                                                                                       >  2:705707/1‑149 (MQ=255)
           ttATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGACGCTATCTAGCTGGTTGGCTTaa                                                                                                                                     >  2:1206953/1‑149 (MQ=255)
             aTTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTaaaa                                                                                                                                   <  1:1159549/149‑1 (MQ=255)
                 gCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGc                                                                                                                               <  1:2449191/149‑1 (MQ=255)
                  ccGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCt                                                                                                                              <  1:2754243/149‑1 (MQ=255)
                  ccGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCt                                                                                                                              <  2:1611995/149‑1 (MQ=255)
                          gcTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATc                                                                                                                      >  1:1739996/1‑149 (MQ=255)
                          gcTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATc                                                                                                                      <  1:1729718/149‑1 (MQ=255)
                          gcTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATc                                                                                                                      >  1:3292409/1‑149 (MQ=255)
                          gcTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATc                                                                                                                      >  1:1532746/1‑149 (MQ=255)
                            ttCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCACTGGCGCATCgg                                                                                                                    >  1:3427617/1‑149 (MQ=255)
                              cctccGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCt                                                                                                                  >  1:1848908/1‑149 (MQ=255)
                              cctccGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCt                                                                                                                  >  1:3113828/1‑149 (MQ=255)
                               ctccGGCAGAAGGGATGACCGTGCGTGAGCTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTg                                                                                                                 >  1:1111169/1‑149 (MQ=255)
                               ctccGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTg                                                                                                                 >  1:188168/1‑149 (MQ=255)
                               ctccGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTg                                                                                                                 >  1:2294118/1‑149 (MQ=255)
                                tccGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTgg                                                                                                                <  1:2050707/149‑1 (MQ=255)
                                tccGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCGGGTTGGCTTAAAACCGCTGGCGCATCGGCTgg                                                                                                                >  2:2175442/1‑149 (MQ=255)
                                 ccGGCAGTAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCGGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGt                                                                                                               <  2:440022/149‑1 (MQ=255)
                                 ccGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGt                                                                                                               >  2:909462/1‑149 (MQ=255)
                                 ccGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGt                                                                                                               >  1:1103140/1‑149 (MQ=255)
                                    gCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGa                                                                                                            >  1:2776990/1‑149 (MQ=255)
                                     cAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGAt                                                                                                           <  1:4042888/149‑1 (MQ=255)
                                      aGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATa                                                                                                          >  1:3127215/1‑149 (MQ=255)
                                         aGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGtc                                                                                                       <  2:3125770/149‑1 (MQ=255)
                                          gggATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGtct                                                                                                      >  2:3051342/1‑149 (MQ=255)
                                             aTGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGtctctc                                                                                                   <  1:1178472/149‑1 (MQ=255)
                                              tGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTctctct                                                                                                  <  1:2941609/149‑1 (MQ=255)
                                               gACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTg                                                                                                 <  1:1334463/149‑1 (MQ=255)
                                               gACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTg                                                                                                 >  2:3076299/1‑149 (MQ=255)
                                               gACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTg                                                                                                 <  1:2046140/149‑1 (MQ=255)
                                                aCCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTgg                                                                                                >  2:287337/1‑149 (MQ=255)
                                                aCCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTgg                                                                                                <  2:2843940/149‑1 (MQ=255)
                                                        tgAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAAc                                                                                        >  1:562779/1‑149 (MQ=255)
                                                        tgAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAAc                                                                                        >  1:2218235/1‑149 (MQ=255)
                                                           aCTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTc                                                                                     <  1:891307/149‑1 (MQ=255)
                                                           aCTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTc                                                                                     <  2:2235734/149‑1 (MQ=255)
                                                            cTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCa                                                                                    <  2:783577/149‑1 (MQ=255)
                                                            cTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCa                                                                                    <  2:112020/149‑1 (MQ=255)
                                                               gtggCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCg                                                                                 <  2:2213498/149‑1 (MQ=255)
                                                               gtggCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCg                                                                                 <  2:4085297/149‑1 (MQ=255)
                                                                tggCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGACAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCgg                                                                                <  2:554871/149‑1 (MQ=255)
                                                                        ggTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGAt                                                                        >  1:501434/1‑149 (MQ=255)
                                                                         gTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATc                                                                       >  1:3075718/1‑149 (MQ=255)
                                                                            gTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGcc                                                                    <  2:3127215/149‑1 (MQ=255)
                                                                            gTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGcc                                                                    >  1:3840387/1‑149 (MQ=255)
                                                                            gTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGcc                                                                    <  1:825086/149‑1 (MQ=255)
                                                                             ttACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCa                                                                   <  2:2850105/149‑1 (MQ=255)
                                                                             ttACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCa                                                                   >  1:2685506/1‑149 (MQ=255)
                                                                             ttACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCa                                                                   >  2:3396449/1‑149 (MQ=255)
                                                                              tACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCAt                                                                  <  1:3939191/149‑1 (MQ=255)
                                                                              tACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCAt                                                                  <  1:995976/149‑1 (MQ=255)
                                                                                 ccGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCt                                                               <  2:1836441/149‑1 (MQ=255)
                                                                                   gTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCtgg                                                             >  2:3658440/1‑149 (MQ=255)
                                                                                   gTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCtgg                                                             >  1:1027771/1‑149 (MQ=255)
                                                                                     ggCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCtggtg                                                           <  1:3555152/149‑1 (MQ=255)
                                                                                       cATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGgc                                                         <  2:3498105/149‑1 (MQ=255)
                                                                                         tGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGgcgc                                                       >  1:2041701/1‑149 (MQ=255)
                                                                                           gcgcgcTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAg                                                     <  2:3113809/149‑1 (MQ=255)
                                                                                              cgcTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGAt                                                  <  2:518957/149‑1 (MQ=255)
                                                                                              cgcTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGAt                                                  <  2:2703960/149‑1 (MQ=255)
                                                                                                  ggggCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGcc                                              >  1:222704/1‑149 (MQ=255)
                                                                                                   gggCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCg                                             >  1:3320249/1‑149 (MQ=255)
                                                                                                    ggCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGt                                            <  1:2879797/149‑1 (MQ=255)
                                                                                                    ggCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGt                                            <  1:3800161/149‑1 (MQ=255)
                                                                                                     gcgcTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGtt                                           <  1:3945644/149‑1 (MQ=255)
                                                                                                       gcTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGt                                         <  1:2627068/149‑1 (MQ=255)
                                                                                                         tttGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCt                                       <  2:4045000/149‑1 (MQ=255)
                                                                                                            ggggCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGtt                                    <  2:1353814/149‑1 (MQ=255)
                                                                                                            ggggCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGtt                                    <  2:3402916/149‑1 (MQ=255)
                                                                                                               gcggcAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCt                                 >  2:1869778/1‑149 (MQ=255)
                                                                                                               gcggcAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCt                                 <  2:3288843/149‑1 (MQ=255)
                                                                                                                  gcAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTcga                              <  1:540758/149‑1 (MQ=255)
                                                                                                                  gcAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTcga                              <  2:1791413/149‑1 (MQ=255)
                                                                                                                  gcAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTcga                              <  2:3523497/149‑1 (MQ=255)
                                                                                                                   cAGATCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTcgac                             <  2:2636139/149‑1 (MQ=255)
                                                                                                                       tCGCGAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAAc                         <  1:3051342/149‑1 (MQ=255)
                                                                                                                          cgAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGa                      >  2:1660931/1‑149 (MQ=255)
                                                                                                                           gAAAAAGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGAc                     >  2:3042570/1‑149 (MQ=255)
                                                                                                                                aGTCGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCgg                >  1:687556/1‑149 (MQ=255)
                                                                                                                                   cGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGgcgc             <  2:3588440/149‑1 (MQ=255)
                                                                                                                                   cGAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGgcgc             <  2:2264633/149‑1 (MQ=255)
                                                                                                                                    gAGGAAGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCt            >  1:2419235/1‑149 (MQ=255)
                                                                                                                                        aaGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGat        <  1:1000327/149‑1 (MQ=255)
                                                                                                                                        aaGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGat        >  1:590091/1‑149 (MQ=255)
                                                                                                                                        aaGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGat        <  2:1014607/149‑1 (MQ=255)
                                                                                                                                         aGCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGata       >  1:751321/1‑149 (MQ=255)
                                                                                                                                          gCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGatat      <  1:3353470/149‑1 (MQ=255)
                                                                                                                                          gCTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGatat      >  2:2485718/1‑149 (MQ=255)
                                                                                                                                           cTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGATATc     <  1:528228/149‑1 (MQ=255)
                                                                                                                                           cTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGATATc     <  1:2594451/149‑1 (MQ=255)
                                                                                                                                           cTATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGATATc     <  1:2579579/149‑1 (MQ=255)
                                                                                                                                            tATCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGATATCg    >  1:3262409/1‑149 (MQ=255)
                                                                                                                                             aTCTAGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGATATCGc   >  1:421859/1‑149 (MQ=255)
                                                                                                                                              tctaGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGATATCGcc  >  2:1073517/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGACCTCGGCGCTGGATATCGCC  >  CP009273/166517‑166807

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: