Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 784,482 G→T 12.6% I318I (ATC→ATA)  galK ← galactokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273784,4820GT12.6% 996.8 / 78.8 285I318I (ATC→ATA) galKgalactokinase
Reads supporting (aligned to +/- strand):  ref base G (130/119);  new base T (20/16);  total (150/135)
Fisher's exact test for biased strand distribution p-value = 7.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 292 total aligned reads displayed.

TCTTTAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATACGTTTCAGGTCGCCTTGCTCCAGCGCGCTGGCAGCTTCAACGGTGCGGGCGTTTTCAGTCAG  >  CP009273/784333‑784631
                                                                                                                                                     |                                                                                                                                                     
tctTTAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACg                                                                                                                                                       <  1:356480/150‑1 (MQ=255)
tctTTAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACg                                                                                                                                                       <  1:3484279/150‑1 (MQ=255)
 ctTTAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGa                                                                                                                                                      >  1:148073/1‑150 (MQ=255)
 ctTTAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGa                                                                                                                                                      <  1:969291/150‑1 (MQ=255)
  tttAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGAt                                                                                                                                                     <  1:902721/150‑1 (MQ=255)
   ttAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGAtt                                                                                                                                                    >  1:2806391/1‑150 (MQ=255)
    tAATACCTGTTTTTGCTTCATATTGTTCAGCGTCAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCAACTTTGTCGCCAATCACAGCTTTGACGAttt                                                                                                                                                   <  2:1516295/150‑1 (MQ=255)
    tAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGAttt                                                                                                                                                   >  2:3459455/1‑150 (MQ=255)
     aaTACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTc                                                                                                                                                  >  1:1987995/1‑150 (MQ=255)
       tACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTa                                                                                                                                                >  2:2772377/1‑150 (MQ=255)
       tACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTa                                                                                                                                                <  2:803416/150‑1 (MQ=255)
        aCCTGTTTTTTCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACTATTTCTAc                                                                                                                                               >  1:819691/1‑150 (MQ=255)
        aCCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTAc                                                                                                                                               <  2:245729/150‑1 (MQ=255)
           tGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCag                                                                                                                                            <  2:472652/150‑1 (MQ=255)
              ttttGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAgt                                                                                                                                         <  1:4119869/150‑1 (MQ=255)
               tttGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAgtg                                                                                                                                        >  1:3767411/1‑150 (MQ=255)
                ttGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACTATTTCTACCAGAgtgt                                                                                                                                       >  2:2925677/1‑150 (MQ=255)
                 tGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTc                                                                                                                                      >  1:724779/1‑150 (MQ=255)
                 tGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTc                                                                                                                                      >  1:604063/1‑150 (MQ=255)
                 tGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTc                                                                                                                                      >  2:2820029/1‑150 (MQ=255)
                   cTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCaa                                                                                                                                    >  2:3404023/1‑150 (MQ=255)
                   cTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCaa                                                                                                                                    <  1:407912/150‑1 (MQ=255)
                   cTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCaa                                                                                                                                    <  1:3301083/150‑1 (MQ=255)
                   cTTCATATTGTTCAGCGACAGCTTGCTGTAAGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCaa                                                                                                                                    >  1:4290954/1‑150 (MQ=255)
                    ttCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACTATTTCTACCAGAGTGTCAAt                                                                                                                                   >  2:4272262/1‑150 (MQ=255)
                    ttCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAAt                                                                                                                                   >  1:1072197/1‑150 (MQ=255)
                    ttCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAAt                                                                                                                                   >  1:1151347/1‑150 (MQ=255)
                     tCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAAtt                                                                                                                                  <  1:1049725/150‑1 (MQ=255)
                       atatTGTTCAGCGACAGCTTGTTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTg                                                                                                                                >  2:2575394/1‑150 (MQ=255)
                       atatTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACTATTTCTACCAGAGTGTCAATTTg                                                                                                                                >  2:2089346/1‑150 (MQ=255)
                       atatTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTg                                                                                                                                >  2:3017336/1‑150 (MQ=255)
                       atatTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTg                                                                                                                                >  2:2753242/1‑150 (MQ=255)
                         atTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCg                                                                                                                              >  1:2983604/1‑150 (MQ=255)
                          ttGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCgg                                                                                                                             <  1:2671568/150‑1 (MQ=255)
                              tCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACg                                                                                                                         <  2:1852078/150‑1 (MQ=255)
                                 gCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACTATTTCTACCAGAGTGTCAATTTGCGGCACGGTg                                                                                                                      <  1:2297242/150‑1 (MQ=255)
                                 gCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTg                                                                                                                      <  1:4220304/150‑1 (MQ=255)
                                  cGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGa                                                                                                                     <  1:472436/150‑1 (MQ=255)
                                   gACAGCTTGATGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGAt                                                                                                                    >  2:4349877/1‑150 (MQ=255)
                                     cAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGAttt                                                                                                                  >  1:2924663/1‑150 (MQ=255)
                                      aGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTc                                                                                                                 >  1:5355/1‑150 (MQ=255)
                                        cTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGCTTTCGa                                                                                                               >  2:3102685/1‑150 (MQ=255)
                                        cTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGa                                                                                                               >  2:2163224/1‑150 (MQ=255)
                                          tGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGaaa                                                                                                             >  1:2879716/1‑150 (MQ=255)
                                               tACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACTATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAatcatc                                                                                                        >  1:3992112/1‑150 (MQ=255)
                                               tACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAatcatc                                                                                                        >  2:3739737/1‑150 (MQ=255)
                                                 cGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATcgc                                                                                                      <  1:2458839/150‑1 (MQ=255)
                                                    caggcaCCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCa                                                                                                   >  2:2212197/1‑150 (MQ=255)
                                                      ggcaCCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATa                                                                                                 <  2:1739178/150‑1 (MQ=255)
                                                       gcaCCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACTATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATag                                                                                                >  2:1072825/1‑150 (MQ=255)
                                                       gcaCCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATag                                                                                                <  2:2781397/150‑1 (MQ=255)
                                                          ccAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACTATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAgg                                                                                             >  2:261826/1‑150 (MQ=255)
                                                           cAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGc                                                                                            >  2:18901/1‑150 (MQ=255)
                                                            aGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCa                                                                                           <  1:1242074/150‑1 (MQ=255)
                                                              ctctTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATg                                                                                         >  2:1287648/1‑150 (MQ=255)
                                                               tctTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACTATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATga                                                                                        <  1:3649478/150‑1 (MQ=255)
                                                                   ccGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACt                                                                                    >  2:3649922/1‑150 (MQ=255)
                                                                   ccGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACt                                                                                    >  1:3197591/1‑150 (MQ=255)
                                                                   ccGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACt                                                                                    >  1:4580477/1‑150 (MQ=255)
                                                                        aTCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACTATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGcc                                                                               <  1:3141260/150‑1 (MQ=255)
                                                                        aTCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACTATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGcc                                                                               <  1:4175962/150‑1 (MQ=255)
                                                                               cgacgaTACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACt                                                                        <  2:3153491/150‑1 (MQ=255)
                                                                                gacgaTACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTc                                                                       >  2:656442/1‑150 (MQ=255)
                                                                                gacgaTACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTc                                                                       <  1:4397004/150‑1 (MQ=255)
                                                                                gacgaTACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTc                                                                       <  1:3017336/150‑1 (MQ=255)
                                                                                   gaTACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGcc                                                                    >  2:2896457/1‑150 (MQ=255)
                                                                                   gaTACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGcc                                                                    >  1:949147/1‑150 (MQ=255)
                                                                                    aTACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGccc                                                                   >  2:1388869/1‑150 (MQ=255)
                                                                                     tACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCa                                                                  >  1:1644180/1‑150 (MQ=255)
                                                                                        aGCCGCCAAATCCGCCGCCGGTTATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATAc                                                               <  1:1559686/150‑1 (MQ=255)
                                                                                           cgccAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATACGtt                                                            >  1:3731095/1‑150 (MQ=255)
                                                                                            gccAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATACGttt                                                           <  2:1469922/150‑1 (MQ=255)
                                                                                              cAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATACGTTTCa                                                         >  2:55324/1‑150 (MQ=255)
                                                                                                  tCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATACGTATCAGGTc                                                     >  2:3733466/1‑150 (MQ=255)
                                                                                                     gccgccgGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATACGTTTCAGGTCGcc                                                  >  1:1867539/1‑150 (MQ=255)
                                                                                                        gccgGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATACGTTTCAGGTCGCCTTg                                               >  1:3692527/1‑150 (MQ=255)
                                                                                                             tCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATACGTTTCAGGTCGCCTTGCTCCa                                          >  2:4514144/1‑150 (MQ=255)
                                                                                                             tCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATACGTTTCAGGTCGCCTTGCTCCa                                          >  1:2140836/1‑150 (MQ=255)
                                                                                                              cATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATACGTTTCAGGTCGCCTTGCTCCAg                                         >  2:1985430/1‑150 (MQ=255)
                                                                                                               aTGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATACGTTTCAGGTCGCCTTGCTCCAgc                                        <  2:2316078/150‑1 (MQ=255)
                                                                                                                tGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATACGTTTCAGGTCGCCTTGCTCCAgcg                                       <  2:3182831/150‑1 (MQ=255)
                                                                                                                   gTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCTCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATACGTTTCAGGTCGCCTTGCTCCAgcgcgc                                    >  1:2537808/1‑150 (MQ=255)
                                                                                                                        ccACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATACGTTTCAGGTCGCCTTGCTCCAGCGCGCTGGCa                               <  2:2997881/150‑1 (MQ=255)
                                                                                                                        ccACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGGGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATACGTTTCAGGTCGCCTTGCTCCAGCGCGCTGGCa                               <  1:3856691/150‑1 (MQ=255)
                                                                                                                         cACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATACGTTTCAGGTCGCCTTGCTCCAGCGCGCTGGCAg                              <  1:1297532/150‑1 (MQ=255)
                                                                                                                         cACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATACGTTTCAGGTCGCCTTGCTCCAGCGCGCTGGCAg                              >  2:90736/1‑150 (MQ=255)
                                                                                                                            cTTTGTCGCCAATCACAGCTTTGACTATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATACGTTTCAGGTCGCCTTGCTCCAGCGCGCTGGCAGCtt                           <  2:1931875/150‑1 (MQ=255)
                                                                                                                               tGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATACGTTTCTGGTCGCCTTGCTCCAGCGCGCTGGCAGCTTCaa                        >  2:363143/1‑150 (MQ=255)
                                                                                                                                     cAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATACGTTTCAGGTCGCCTTGCTCCAGCGCGCTGGCAGCTTCAACGGTGc                  >  1:2469332/1‑150 (MQ=255)
                                                                                                                                       aTCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATACGTTTCAGGTCGCCTTGCTCCAGCGCGCTGGCAGCTTCAACGGTGCgg                >  1:3377166/1‑150 (MQ=255)
                                                                                                                                       aTCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATACGTTTCAGGTCGCCTTGCTCCAGCGCGCGGGCAGCTTCAACGGTGCgg                <  2:4614341/150‑1 (MQ=255)
                                                                                                                                        tCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATACGTTTCAGGTCGCCTTGCTCCAGCGCGCTGGCAGCTTCAACGGTGCggg               >  1:3427536/1‑150 (MQ=255)
                                                                                                                                        tCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATACGTTTCAGGTCGCCTTGCTCCAGCGCGCTGGCAGCTTCAACGGTGCggg               <  2:2852015/150‑1 (MQ=255)
                                                                                                                                        tCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATACGTTTCAGGTCGCCTTGCTCCAGCGCGCTGGCAGCTTCAACGGTGCggg               >  1:4134034/1‑150 (MQ=255)
                                                                                                                                         cacaGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATACGTTTCAGGTCGCCTTGCTCCAGCGCGCTGGCAGCTTCAACGGTGCGGGc              >  1:1600217/1‑150 (MQ=255)
                                                                                                                                               tttGACTATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATACGTTTCAGGTCGCCTTGCTCCAGCGCGCTGGCAGCTTCAACGGTGCGGGCGTTTtc        >  2:2373784/1‑150 (MQ=255)
                                                                                                                                                  gACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATACGTTTCAGGTCGCCTTGCTCCAGCGCGCTGGCAGCTTCAACGGTGCGGGCGTTTtcagt     >  2:4154211/1‑150 (MQ=255)
                                                                                                                                                   aCGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATACGTTTCAGGTCGCCTTGCTCCAGCGCGCTGGCAGCTTCAACGGTGCGGGCGTTTtcagtc    <  2:4489454/150‑1 (MQ=255)
                                                                                                                                                    cGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATACGTTTCAGGTCGCCTTGCTCCAGCGCGCTGGCAGCTTCAACGGTGCGGGCGTTTtcagtca   >  1:2448457/1‑150 (MQ=255)
                                                                                                                                                     gATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATACGTTTCAGGTCGCCTTGCTCCAGCGCGCTGGCAGCTTCAACGGTGCGGGCGTTTtcagtcag  >  2:3534679/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCTTTAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATACGTTTCAGGTCGCCTTGCTCCAGCGCGCTGGCAGCTTCAACGGTGCGGGCGTTTTCAGTCAG  >  CP009273/784333‑784631

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: ‑