Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,850,839 G→T 30.0% P309Q (CCG→CAG)  yidJ ← sulfatase/phosphatase superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,850,8390GT30.0% 610.6 / 275.7 322P309Q (CCG→CAG) yidJsulfatase/phosphatase superfamily protein
Reads supporting (aligned to +/- strand):  ref base G (109/116);  new base T (53/44);  total (162/160)
Fisher's exact test for biased strand distribution p-value = 3.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 327 total aligned reads displayed.

GTGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTGGTCA  >  CP009273/3850690‑3850988
                                                                                                                                                     |                                                                                                                                                     
gTGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCg                                                                                                                                                       >  2:4366731/1‑150 (MQ=255)
gTGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCg                                                                                                                                                       >  1:4904501/1‑150 (MQ=255)
gTGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCg                                                                                                                                                       >  2:2294066/1‑150 (MQ=255)
  ggCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGa                                                                                                                                                     >  1:1382721/1‑150 (MQ=255)
   gCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGaa                                                                                                                                                    <  1:3354578/150‑1 (MQ=255)
     tctTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATg                                                                                                                                                  <  2:1348797/150‑1 (MQ=255)
     tctTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATg                                                                                                                                                  <  1:4042818/150‑1 (MQ=255)
       tttCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCg                                                                                                                                                >  2:2766773/1‑150 (MQ=255)
       tttCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCg                                                                                                                                                <  1:267577/150‑1 (MQ=255)
           aCGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTg                                                                                                                                            >  1:2122489/1‑150 (MQ=255)
            cGGCAAGGATATTTTCCCCCGGCATAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGa                                                                                                                                           <  2:4666278/150‑1 (MQ=255)
            cGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGa                                                                                                                                           <  2:4127909/150‑1 (MQ=255)
             ggCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGAt                                                                                                                                          <  2:564328/150‑1 (MQ=255)
             ggCAAGGATATGTTCCCCCGGCAGAATCTCTGGATTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGAt                                                                                                                                          <  2:3281596/150‑1 (MQ=255)
              gCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATg                                                                                                                                         >  1:752575/1‑150 (MQ=255)
              gCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTGGCGGCGAACGGATGATCAGCGGAATGCGGGTGATg                                                                                                                                         >  1:1767916/1‑150 (MQ=255)
               cAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATgt                                                                                                                                        <  2:4976728/150‑1 (MQ=255)
                 aGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATgtcg                                                                                                                                      >  1:861888/1‑150 (MQ=255)
                   gATATTTTCCCCCGGCAGAATGTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATgtcgtc                                                                                                                                    >  2:4475827/1‑150 (MQ=255)
                   gATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATgtcgtc                                                                                                                                    >  2:63321/1‑150 (MQ=255)
                    atatTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCa                                                                                                                                   >  2:319168/1‑150 (MQ=255)
                     tatTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTACATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCAt                                                                                                                                  >  1:4767904/1‑150 (MQ=255)
                     tatTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCAt                                                                                                                                  <  2:2988348/150‑1 (MQ=255)
                       ttttCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATAc                                                                                                                                >  2:2610890/1‑150 (MQ=255)
                           cccccGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATcg                                                                                                                            <  2:1053438/150‑1 (MQ=255)
                              ccGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATcgccg                                                                                                                         <  1:3124824/150‑1 (MQ=255)
                               cGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATcgccgc                                                                                                                        >  1:1441753/1‑150 (MQ=255)
                               cGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATcgccgc                                                                                                                        >  2:3799473/1‑150 (MQ=255)
                                  cAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGcccc                                                                                                                     >  1:815282/1‑150 (MQ=255)
                                   aGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCt                                                                                                                    <  1:3864441/150‑1 (MQ=255)
                                    gAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCtt                                                                                                                   <  2:1832977/150‑1 (MQ=255)
                                     aaTCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCttt                                                                                                                  >  2:1107044/1‑150 (MQ=255)
                                     aaTCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCttt                                                                                                                  >  2:3724145/1‑150 (MQ=255)
                                     aaTCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCttt                                                                                                                  <  1:4963913/150‑1 (MQ=255)
                                     aaTCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCttt                                                                                                                  >  2:1428908/1‑150 (MQ=255)
                                      aTCTCTGGTTTTTCTATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTa                                                                                                                 <  1:5161786/150‑1 (MQ=255)
                                      aTCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTa                                                                                                                 >  1:2242218/1‑150 (MQ=255)
                                          ctGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGa                                                                                                             >  2:2202172/1‑150 (MQ=255)
                                          ctGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGa                                                                                                             <  1:2947510/150‑1 (MQ=255)
                                          cnGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGa                                                                                                             >  1:3314369/1‑150 (MQ=255)
                                           tGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGAt                                                                                                            >  1:4877173/1‑150 (MQ=255)
                                           tGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGAt                                                                                                            <  1:4292434/150‑1 (MQ=255)
                                            ggTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATc                                                                                                           <  2:2037076/150‑1 (MQ=255)
                                             gTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGACCGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCa                                                                                                          >  1:538540/1‑150 (MQ=255)
                                               ttttCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGc                                                                                                        >  1:4687187/1‑150 (MQ=255)
                                                   cAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTAt                                                                                                    <  2:4433035/150‑1 (MQ=255)
                                                   cAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTAt                                                                                                    <  1:58671/150‑1 (MQ=255)
                                                     atatCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTAtgt                                                                                                  >  1:2307061/1‑150 (MQ=255)
                                                     atatCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTAtgt                                                                                                  >  1:4471263/1‑150 (MQ=255)
                                                      tatCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTAtgtg                                                                                                 <  1:1214912/150‑1 (MQ=255)
                                                         cTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTgcgc                                                                                              <  2:1668883/150‑1 (MQ=255)
                                                              aGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCcatc                                                                                         >  1:1992280/1‑150 (MQ=255)
                                                                cgcCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCcatcat                                                                                       <  1:3620615/150‑1 (MQ=255)
                                                                  ccATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCAttt                                                                                     <  1:3813844/150‑1 (MQ=255)
                                                                    atcatTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCg                                                                                   >  2:268412/1‑150 (MQ=255)
                                                                     tcatTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGc                                                                                  >  1:4352586/1‑150 (MQ=255)
                                                                       atTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCg                                                                                >  2:4444020/1‑150 (MQ=255)
                                                                        ttGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGt                                                                               <  1:1654530/150‑1 (MQ=255)
                                                                        ttGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGAGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTACGCCGt                                                                               >  2:5232098/1‑150 (MQ=255)
                                                                         tGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTg                                                                              <  1:1096265/150‑1 (MQ=255)
                                                                           tCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGAt                                                                            <  2:1763040/150‑1 (MQ=255)
                                                                              gCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCgg                                                                         <  2:3781756/150‑1 (MQ=255)
                                                                                aGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAg                                                                       <  1:299256/150‑1 (MQ=255)
                                                                                   aaaTCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAggaa                                                                    >  1:4587236/1‑148 (MQ=255)
                                                                                   aaaTCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGta                                                                    >  2:4254986/1‑150 (MQ=255)
                                                                                    aaTCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGtat                                                                   <  2:1281096/150‑1 (MQ=255)
                                                                                      tCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTataa                                                                 >  1:3783204/1‑150 (MQ=255)
                                                                                         atatGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTataaata                                                              >  1:1172329/1‑150 (MQ=255)
                                                                                           atGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAAc                                                            >  1:3208704/1‑150 (MQ=255)
                                                                                           atGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAAc                                                            >  1:1525737/1‑150 (MQ=255)
                                                                                             gactgactGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAAccc                                                          >  2:3184608/1‑150 (MQ=255)
                                                                                               ctgactgGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCAc                                                        <  1:4976940/150‑1 (MQ=255)
                                                                                               ctgactgGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCAc                                                        >  1:3358999/1‑150 (MQ=255)
                                                                                                tgactgGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACg                                                       >  2:3424512/1‑150 (MQ=255)
                                                                                                tgactgGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCATTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACg                                                       <  2:364449/150‑1 (MQ=255)
                                                                                                   ctgGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTAt                                                    <  2:3898839/150‑1 (MQ=255)
                                                                                                    tgGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTAtt                                                   <  1:401282/150‑1 (MQ=255)
                                                                                                     ggCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTAttt                                                  <  1:3534325/150‑1 (MQ=255)
                                                                                                              aCCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGc                                         >  2:2396585/1‑150 (MQ=255)
                                                                                                                 tGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGctct                                      >  2:850076/1‑150 (MQ=255)
                                                                                                                 tGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTCTGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGctct                                      <  2:4922737/150‑1 (MQ=255)
                                                                                                                 tGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGctct                                      >  2:1588608/1‑150 (MQ=255)
                                                                                                                 tGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGctct                                      <  2:714336/150‑1 (MQ=255)
                                                                                                                  gTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTg                                     <  2:67939/150‑1 (MQ=255)
                                                                                                                      ccgcTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGt                                 >  1:612284/1‑150 (MQ=255)
                                                                                                                          accccTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAg                             <  2:1144892/149‑1 (MQ=255)
                                                                                                                                tGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATtg                       >  2:2965301/1‑150 (MQ=255)
                                                                                                                                 gcggcgAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATtga                      >  1:102657/1‑150 (MQ=255)
                                                                                                                                 gcggcgAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATtga                      <  2:2888787/150‑1 (MQ=255)
                                                                                                                                  cggcgAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATtgat                     <  1:996757/150‑1 (MQ=255)
                                                                                                                                     cgAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGAc                  <  2:3583988/150‑1 (MQ=255)
                                                                                                                                     cgAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGAc                  <  1:3457258/150‑1 (MQ=255)
                                                                                                                                       aaCGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGAccc                <  1:1857039/150‑1 (MQ=255)
                                                                                                                                          ggATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTc             >  1:4223628/1‑150 (MQ=255)
                                                                                                                                            atgatCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCg           >  1:2088371/1‑150 (MQ=255)
                                                                                                                                             tgatCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGa          >  2:5115999/1‑150 (MQ=255)
                                                                                                                                              gatCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGAt         >  2:4760904/1‑150 (MQ=255)
                                                                                                                                                 cAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTg      <  2:4594357/150‑1 (MQ=255)
                                                                                                                                                 cAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTg      <  2:1988936/150‑1 (MQ=255)
                                                                                                                                                  aGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTgg     >  2:1390435/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTGGTCA  >  CP009273/3850690‑3850988

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: