Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,978,003 C→A 18.8% E455D (GAG→GAT aslA ← putative Ser‑type periplasmic non‑aryl sulfatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,978,0030CA18.8% 876.1 / 130.8 292E455D (GAG→GATaslAputative Ser‑type periplasmic non‑aryl sulfatase
Reads supporting (aligned to +/- strand):  ref base C (132/105);  new base A (32/22);  total (164/128)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 296 total aligned reads displayed.

AACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGAAGAGATCTGCCAG  >  CP009273/3977854‑3978152
                                                                                                                                                     |                                                                                                                                                     
aaCACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGc                                                                                                                                                       >  2:3489195/1‑150 (MQ=255)
  cacCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTc                                                                                                                                                     <  2:1249966/150‑1 (MQ=255)
   acCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCg                                                                                                                                                    <  2:3269771/150‑1 (MQ=255)
     cgacgaTCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGc                                                                                                                                                  <  1:3756414/150‑1 (MQ=255)
       acgaTCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCt                                                                                                                                                <  2:3587533/150‑1 (MQ=255)
            cccGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACgg                                                                                                                                           <  1:2995985/150‑1 (MQ=255)
                    ttGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTg                                                                                                                                   <  2:2999338/150‑1 (MQ=255)
                        aTTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCa                                                                                                                               >  2:2276359/1‑150 (MQ=255)
                             tgtgCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGt                                                                                                                          >  2:4405406/1‑150 (MQ=255)
                              gtgCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGtt                                                                                                                         <  1:1379580/150‑1 (MQ=255)
                               tgCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTc                                                                                                                        >  2:4343496/1‑150 (MQ=255)
                                gCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTcc                                                                                                                       <  2:2175169/150‑1 (MQ=255)
                                gCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTcc                                                                                                                       <  2:2367002/150‑1 (MQ=255)
                                  cGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCa                                                                                                                     >  1:2464193/1‑150 (MQ=255)
                                  cGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCa                                                                                                                     >  1:3970529/1‑150 (MQ=255)
                                    gTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAgg                                                                                                                   >  1:2861172/1‑150 (MQ=255)
                                    gTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAgg                                                                                                                   <  2:1298602/150‑1 (MQ=255)
                                    gTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAgg                                                                                                                   <  2:1692323/150‑1 (MQ=255)
                                    gTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAgg                                                                                                                   >  1:1860743/1‑150 (MQ=255)
                                     tGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGga                                                                                                                  >  1:3348327/1‑150 (MQ=255)
                                       aaTCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGgaag                                                                                                                >  1:4624291/1‑150 (MQ=255)
                                          ccACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGaagaag                                                                                                             >  1:1314858/1‑150 (MQ=255)
                                            aCCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGa                                                                                                           <  2:549121/150‑1 (MQ=255)
                                             cccTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGAt                                                                                                          <  2:506011/150‑1 (MQ=255)
                                             cccTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGAt                                                                                                          >  1:3678789/1‑150 (MQ=255)
                                              ccTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATg                                                                                                         <  1:2694137/150‑1 (MQ=255)
                                                tGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTc                                                                                                       <  2:1952040/150‑1 (MQ=255)
                                                 gATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCt                                                                                                      <  2:1163700/150‑1 (MQ=255)
                                                 gATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCt                                                                                                      <  1:2783676/150‑1 (MQ=255)
                                                 gATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCt                                                                                                      <  1:3154930/150‑1 (MQ=255)
                                                  atatCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTg                                                                                                     <  2:4264588/150‑1 (MQ=255)
                                                   tatCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTgg                                                                                                    <  1:1902015/150‑1 (MQ=255)
                                                    atCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGt                                                                                                   >  1:450052/1‑150 (MQ=255)
                                                    atCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGt                                                                                                   <  2:4232328/150‑1 (MQ=255)
                                                    atCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGt                                                                                                   <  2:2021223/150‑1 (MQ=255)
                                                    atCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGt                                                                                                   >  1:2384794/1‑150 (MQ=255)
                                                    atCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGt                                                                                                   >  2:1630364/1‑150 (MQ=255)
                                                    atCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGt                                                                                                   <  2:1330828/150‑1 (MQ=255)
                                                     tCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTc                                                                                                  <  1:148068/150‑1 (MQ=255)
                                                        gCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCcac                                                                                               <  1:3539736/150‑1 (MQ=255)
                                                          tctGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCacac                                                                                             >  2:3292858/1‑150 (MQ=255)
                                                            tGGGTATAAGCGTAAGGTTGCAGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCATACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGTATGTCTGGTCCACACCa                                                                                           >  2:1038910/1‑150 (MQ=255)
                                                             gggTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCAt                                                                                          >  2:675934/1‑150 (MQ=255)
                                                               gTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCg                                                                                        >  1:1737205/1‑150 (MQ=255)
                                                               gTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCg                                                                                        >  1:985787/1‑150 (MQ=255)
                                                               gTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCg                                                                                        <  1:2281047/150‑1 (MQ=255)
                                                                tataAGCGTAAGGTTGCTGATTCAGGACGTGATAATTGAACTCATCAATACGCACAGCAGCGTGTTTACCGTTGAGGAAGTATTGCTCGGCCTTACGGTTAGACTGTCCATTTTTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGa                                                                                       <  2:1356935/150‑1 (MQ=255)
                                                                tataAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGa                                                                                       <  2:222389/150‑1 (MQ=255)
                                                                tataAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGa                                                                                       >  1:3218345/1‑150 (MQ=255)
                                                                 ataAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGAt                                                                                      <  2:4512109/150‑1 (MQ=255)
                                                                   aaGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATaa                                                                                    >  2:4454378/1‑150 (MQ=255)
                                                                    aGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATaaa                                                                                   <  2:3344140/150‑1 (MQ=255)
                                                                    aGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATaaa                                                                                   >  1:4554682/1‑150 (MQ=255)
                                                                        tAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAggtg                                                                               >  2:460517/1‑150 (MQ=255)
                                                                         aaGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAggtgg                                                                              >  1:813895/1‑150 (MQ=255)
                                                                          aGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGTATGTCTGGCCCACACCATCGATAAAggtggt                                                                             >  2:3997753/1‑150 (MQ=255)
                                                                           ggTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCCTCGATAAAGGTGGtt                                                                            >  2:3810766/1‑150 (MQ=255)
                                                                             ttGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGtttt                                                                          >  2:111190/1‑150 (MQ=255)
                                                                             ttGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGtttt                                                                          >  2:3944684/1‑150 (MQ=255)
                                                                             ttGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTAACAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGtttt                                                                          >  1:2137467/1‑150 (MQ=255)
                                                                              tGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTc                                                                         <  2:3114508/150‑1 (MQ=255)
                                                                                cTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCgg                                                                       <  1:1506377/150‑1 (MQ=255)
                                                                                  gAATCACGACGTTATACATGAACTCAACCATACGCACAGCAGCGAGTTTACCGTTTAGTACTAACTGCTCGGCCATACGGTTAGACTGACCATTTTTTCCCAGGAAGAAGGATGTCTCGTCCACACAATCGATAAAGGTGATTTTCGACa                                                                     >  2:2793837/1‑150 (MQ=255)
                                                                                     tCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTa                                                                  >  1:4564906/1‑150 (MQ=255)
                                                                                      cAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTaa                                                                 <  1:1582070/150‑1 (MQ=255)
                                                                                        ggACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAAt                                                               >  2:1477120/1‑150 (MQ=255)
                                                                                           cGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCg                                                            >  1:2741701/1‑150 (MQ=255)
                                                                                             tGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGcc                                                          >  2:723530/1‑150 (MQ=255)
                                                                                             tGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGcc                                                          >  2:1504296/1‑150 (MQ=255)
                                                                                              gATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCa                                                         <  1:2626480/150‑1 (MQ=255)
                                                                                               aTACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCAc                                                        >  2:149227/1‑150 (MQ=255)
                                                                                                 aCTAGAACTCATCCATACGAACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTAGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACtt                                                      >  2:418781/1‑150 (MQ=255)
                                                                                                  cTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACttt                                                     <  2:795071/150‑1 (MQ=255)
                                                                                                  cTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACttt                                                     >  2:2817967/1‑150 (MQ=255)
                                                                                                     gAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGc                                                  >  1:446411/1‑150 (MQ=255)
                                                                                                      aaCTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCt                                                 >  2:2664628/1‑150 (MQ=255)
                                                                                                         tCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCa                                              <  1:4031155/150‑1 (MQ=255)
                                                                                                          cATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAg                                             <  1:1614189/150‑1 (MQ=255)
                                                                                                           aTCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAgg                                            <  1:723530/150‑1 (MQ=255)
                                                                                                              cATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATg                                         <  2:2421749/150‑1 (MQ=255)
                                                                                                               aTACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAATGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGc                                        >  1:2514748/1‑150 (MQ=255)
                                                                                                               aTACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGc                                        >  2:1632099/1‑150 (MQ=255)
                                                                                                               aTACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGc                                        >  1:3378446/1‑150 (MQ=255)
                                                                                                                    cacaGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGc                                   >  1:55610/1‑150 (MQ=255)
                                                                                                                        gcagcGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGa                               >  2:4371388/1‑150 (MQ=255)
                                                                                                                          agcGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATc                             >  2:1473422/1‑150 (MQ=255)
                                                                                                                          agcGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATc                             >  2:2557853/1‑150 (MQ=255)
                                                                                                                            cGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCa                           >  2:3552837/1‑150 (MQ=255)
                                                                                                                                tttACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAgcgc                       >  1:4083020/1‑150 (MQ=255)
                                                                                                                                  tACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCgg                     >  2:1833199/1‑150 (MQ=255)
                                                                                                                                           ggAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGAaga            <  1:4405783/150‑1 (MQ=255)
                                                                                                                                            gAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGAagag           >  2:4259519/1‑150 (MQ=255)
                                                                                                                                            gAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGAagag           >  1:4293868/1‑150 (MQ=255)
                                                                                                                                            gAAGTAGTGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGAagag           >  1:21009/1‑150 (MQ=255)
                                                                                                                                                tagtGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGAAGAGATCt       >  1:1038803/1‑150 (MQ=255)
                                                                                                                                                 agtGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGAAGAGATCTg      <  2:4449100/150‑1 (MQ=255)
                                                                                                                                                 agtGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGAAGAGATCTg      <  2:332504/150‑1 (MQ=255)
                                                                                                                                                 agtGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGAAGAGATCTg      <  1:1154880/150‑1 (MQ=255)
                                                                                                                                                  gtGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGAAGAGATCTGc     >  1:2936706/1‑150 (MQ=255)
                                                                                                                                                  gtGATCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGAAGAGATCTGc     >  1:3442606/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AACACCGACGATCCCGCCGTTTGCATTACTGTGCCGGTGAATCCACCCTGATATCCGCTCTGGGTATAAGCGTAAGGTTGCTGAATCAGGACGTGATACTTGAACTCATCCATACGCACAGCAGCGAGTTTACCGTTGAGGAAGTAGTGCTCGGCCTTACGGTTAGACTGACCATTTGTTCCCAGGAAGAAGGATGTCTGGTCCACACCATCGATAAAGGTGGTTTTCGGCACTAAATTCGCCACTTTCGCTCCAGGATGCCCTGCCAGATCCAGCGCGGTAGGGAAGAGATCTGCCAG  >  CP009273/3977854‑3978152

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: