Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,995,335 C→A 6.7% V544V (GTC→GTA xdhA → xanthine dehydrogenase, molybdenum binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,995,3350CA 6.7% 1009.1 / 30.2 269V544V (GTC→GTAxdhAxanthine dehydrogenase, molybdenum binding subunit
Reads supporting (aligned to +/- strand):  ref base C (123/128);  new base A (3/15);  total (126/143)
Fisher's exact test for biased strand distribution p-value = 1.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 277 total aligned reads displayed.

CGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGAACGCGGCGGGCAGCTCTCTGCTGAAAG  >  CP009273/2995188‑2995484
                                                                                                                                                   |                                                                                                                                                     
cGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCAt                                                                                                                                                     <  1:1604071/150‑1 (MQ=255)
cGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCAt                                                                                                                                                     <  1:1546608/150‑1 (MQ=255)
cGACGTTCGCGTTATTTCAACTCAAGAAACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCAt                                                                                                                                                     >  2:4403031/1‑150 (MQ=255)
     ttCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTAc                                                                                                                                                >  1:1320949/1‑150 (MQ=255)
       cgcgTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACAt                                                                                                                                              <  2:2887198/150‑1 (MQ=255)
          gTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACAtcag                                                                                                                                           >  2:1571935/1‑150 (MQ=255)
           ttATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACAtcagt                                                                                                                                          <  1:1686685/150‑1 (MQ=255)
             aTTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACAtcagtca                                                                                                                                        >  2:2064730/1‑150 (MQ=255)
              tttCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGTGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACAtcagtcag                                                                                                                                       <  2:2020189/150‑1 (MQ=255)
               ttCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGc                                                                                                                                      >  2:2424685/1‑150 (MQ=255)
               ttCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGc                                                                                                                                      <  2:3229510/150‑1 (MQ=255)
               ttCAACTCAAGAGACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGc                                                                                                                                      <  1:3510154/150‑1 (MQ=255)
                 cAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGa                                                                                                                                    <  1:1439374/150‑1 (MQ=255)
                   aCTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATg                                                                                                                                  >  2:2790662/1‑150 (MQ=255)
                     tCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGaa                                                                                                                                <  1:356153/150‑1 (MQ=255)
                        agagACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCt                                                                                                                             <  1:3141411/150‑1 (MQ=255)
                        aGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCt                                                                                                                             <  1:860594/150‑1 (MQ=255)
                        aGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCt                                                                                                                             >  2:692673/1‑150 (MQ=255)
                        aGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCt                                                                                                                             >  1:2163788/1‑150 (MQ=255)
                        aGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGGTGCCGCGCCTGCGCTGCGGAGTGCGGCACTATTATTAAACGACAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCt                                                                                                                             <  1:4443053/150‑1 (MQ=255)
                          aTACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGa                                                                                                                           <  1:2341295/150‑1 (MQ=255)
                           tACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGAc                                                                                                                          >  1:1861491/1‑150 (MQ=255)
                            aCCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTAATGCTACATCAGTCAGCGATGAATCTGAcc                                                                                                                         <  2:1350917/150‑1 (MQ=255)
                              cGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCt                                                                                                                       >  1:1039917/1‑150 (MQ=255)
                               gACGTTACGCCGTTCGATCCCGGCGCTTTTGCCTCTCTCCAGAGCTATGTTGCCGCGACTGCTCTGCGCAGTGCGGCACTATTTTTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGTATCTGACCCTg                                                                                                                      <  2:35191/150‑1 (MQ=255)
                               gACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTg                                                                                                                      >  1:1743958/1‑150 (MQ=255)
                               gACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTg                                                                                                                      <  2:35002/150‑1 (MQ=255)
                                  gTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATa                                                                                                                   >  1:2244843/1‑150 (MQ=255)
                                   ttACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATaa                                                                                                                  >  1:1084792/1‑150 (MQ=255)
                                    tACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATaaa                                                                                                                 >  2:3848695/1‑150 (MQ=255)
                                     aCGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATaaaa                                                                                                                <  1:3423218/150‑1 (MQ=255)
                                       gCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAgg                                                                                                              >  2:4302239/1‑150 (MQ=255)
                                          gTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCa                                                                                                           >  1:4133293/1‑150 (MQ=255)
                                          gTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCa                                                                                                           >  1:2762158/1‑150 (MQ=255)
                                            tCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCata                                                                                                         <  1:692673/150‑1 (MQ=255)
                                             cGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCatat                                                                                                        >  2:3077329/1‑150 (MQ=255)
                                              gATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTAATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATc                                                                                                       >  1:3183341/1‑150 (MQ=255)
                                               aTCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTAATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCg                                                                                                      <  2:4018264/150‑1 (MQ=255)
                                                   cGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCt                                                                                                  <  1:10027/150‑1 (MQ=255)
                                                         aTTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGa                                                                                            >  2:1139421/1‑150 (MQ=255)
                                                         aTTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGa                                                                                            <  1:2644600/150‑1 (MQ=255)
                                                          tttGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTAATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGaa                                                                                           <  1:3809389/150‑1 (MQ=255)
                                                           ttGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCTCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTAATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAAc                                                                                          <  1:2910695/150‑1 (MQ=255)
                                                                       agagCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGc                                                                              >  2:3999388/1‑150 (MQ=255)
                                                                          gCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGt                                                                           <  1:1194861/150‑1 (MQ=255)
                                                                             aTGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTaa                                                                        >  2:1460944/1‑150 (MQ=255)
                                                                                 tGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTc                                                                    <  2:130313/150‑1 (MQ=255)
                                                                                  gCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCg                                                                   >  2:4196650/1‑150 (MQ=255)
                                                                                    cgcgCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGGGCTGGTTGAACGACCGGAAGAGCCGTTCATGTCGtt                                                                 >  2:4222978/1‑150 (MQ=255)
                                                                                    cgcgCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTAATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGtt                                                                 <  1:3450372/150‑1 (MQ=255)
                                                                                      cgcCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTaa                                                               >  2:1509974/1‑150 (MQ=255)
                                                                                       gcCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTaaa                                                              >  2:1578471/1‑150 (MQ=255)
                                                                                           gcgctgcgcAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTGAATGTCGTTAAAAGAt                                                          >  2:1736272/1‑150 (MQ=255)
                                                                                            cgctgcgcAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGAtt                                                         <  2:1971852/150‑1 (MQ=255)
                                                                                               tgcgcAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTgg                                                      <  1:1475515/150‑1 (MQ=255)
                                                                                               tgcgcAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTgg                                                      <  2:2495490/150‑1 (MQ=255)
                                                                                                    aGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATg                                                 <  2:716065/150‑1 (MQ=255)
                                                                                                     gTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATgg                                                <  1:3140816/150‑1 (MQ=255)
                                                                                                       gCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGAc                                              <  1:618518/150‑1 (MQ=255)
                                                                                                        cGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACg                                             <  2:1205892/150‑1 (MQ=255)
                                                                                                           cACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCTATGAATCTGACAATGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCtt                                          <  1:4439680/150‑1 (MQ=255)
                                                                                                           cACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCtt                                          <  1:2203363/150‑1 (MQ=255)
                                                                                                            aCTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCttt                                         <  1:1365595/150‑1 (MQ=255)
                                                                                                            aCTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCttt                                         <  1:1365146/150‑1 (MQ=255)
                                                                                                            aCTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCttt                                         >  2:2616154/1‑150 (MQ=255)
                                                                                                            aCTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCttt                                         >  1:1170980/1‑150 (MQ=255)
                                                                                                             cTATTATTACAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTc                                        <  2:943626/150‑1 (MQ=255)
                                                                                                             cTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTc                                        >  1:4315432/1‑150 (MQ=255)
                                                                                                             cTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTc                                        >  1:398755/1‑150 (MQ=255)
                                                                                                                ttattaAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTac                                     >  1:646933/1‑150 (MQ=255)
                                                                                                                 tattaAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTacc                                    <  2:1906998/150‑1 (MQ=255)
                                                                                                                  attaAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTacca                                   >  2:3826928/1‑150 (MQ=255)
                                                                                                                  attaAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTacca                                   >  2:1759310/1‑150 (MQ=255)
                                                                                                                  attaAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTacca                                   >  2:1759224/1‑150 (MQ=255)
                                                                                                                  attaAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTacca                                   >  2:866592/1‑150 (MQ=255)
                                                                                                                  attaAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGGTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTacca                                   >  2:3827020/1‑150 (MQ=255)
                                                                                                                   ttaAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTaccac                                  <  1:2369581/150‑1 (MQ=255)
                                                                                                                    taAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTaccacc                                 <  1:3716124/150‑1 (MQ=255)
                                                                                                                     aaaaGAGAAACTCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCAccc                                <  2:2968333/150‑1 (MQ=255)
                                                                                                                     aaaaGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCAccc                                >  1:3060307/1‑150 (MQ=255)
                                                                                                                     aaaaGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCAccc                                >  1:2505494/1‑150 (MQ=255)
                                                                                                                       aaGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTg                              <  1:449236/150‑1 (MQ=255)
                                                                                                                       aaGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTg                              <  1:1370239/150‑1 (MQ=255)
                                                                                                                         gagaAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGaa                            <  2:1787663/150‑1 (MQ=255)
                                                                                                                           gaAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGAAcg                          <  2:737203/150‑1 (MQ=255)
                                                                                                                            aaaaTCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCATGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGAAcgc                         >  2:3067091/1‑150 (MQ=255)
                                                                                                                             aaaTCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGAAcgcg                        <  2:2317606/150‑1 (MQ=255)
                                                                                                                             aaaTCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGAAcgcg                        >  1:1009405/1‑150 (MQ=255)
                                                                                                                             aaaTCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGAAcgcg                        >  1:1241202/1‑150 (MQ=255)
                                                                                                                             aaaTCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGAAcgcg                        >  1:3288001/1‑150 (MQ=255)
                                                                                                                                 catcGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGAACGcggcgg                    <  2:2164441/150‑1 (MQ=255)
                                                                                                                                    cGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGAACGCGGCGGGCa                 <  1:2082448/150‑1 (MQ=255)
                                                                                                                                           gccgcAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGAACGCGGCGGGCAGctctct          <  2:2645672/150‑1 (MQ=255)
                                                                                                                                           gccgcAGTAATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGAACGCGGCGGGCAGctctct          <  1:4161334/150‑1 (MQ=255)
                                                                                                                                            ccgcAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGAACGCGGCGGGCAGCTCTctg         >  1:4455842/1‑150 (MQ=255)
                                                                                                                                                aGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGAACGCGGCGGGCAGCTCTCTGCTGa     >  1:9567/1‑150 (MQ=255)
                                                                                                                                                 gTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGAACGCGGCGGGCAGCTCTCTGCTGaa    >  1:17592/1‑150 (MQ=255)
                                                                                                                                                  tCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGAACGCGGCGGGCAGCTCTCTGCTGaaa   <  1:3376099/150‑1 (MQ=255)
                                                                                                                                                   cATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGAACGCGGCGGGCAGCTCTCTGCTGAAAg  <  1:3342150/150‑1 (MQ=255)
                                                                                                                                                   cATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACAGGAAGAGCCGTTAATGTCGTTAAAAGATTGGGCGATGGACGCTTTCTACCACCCTGAACGCGGCGGGCAGCTCTCTGCTGAAAg  <  1:2928387/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTTGCCTCACGCCAGAGCTATGTTGCCGCGCCTGCGCTGCGCAGTGCGGCACTATTATTAAAAGAGAAAATCATCGCTCACGCCGCAGTCATGCTACATCAGTCAGCGATGAATCTGACCCTGATAAAAGGCCATATCGTGCTGGTTGAACGACCGGAAGAGCCGTTAATGTCGTTAAAAGATTTGGCGATGGACGCTTTCTACCACCCTGAACGCGGCGGGCAGCTCTCTGCTGAAAG  >  CP009273/2995188‑2995484

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: