Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,834,943 C→T 8.1% S407F (TCC→TTC)  ynjE → molybdopterin synthase sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,834,9430CT 8.1% 616.3 / 31.3 198S407F (TCC→TTC) ynjEmolybdopterin synthase sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base C (99/83);  new base T (5/11);  total (104/94)
Fisher's exact test for biased strand distribution p-value = 1.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 204 total aligned reads displayed.

TGGCACCATGCGCAGCGCCGATGATATTACCGCTATGTGGAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAA  >  CP009273/1834794‑1835092
                                                                                                                                                     |                                                                                                                                                     
tGGCACCATGCGCAGCGCCGATGATATTACCGCTATGTGGAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTc                                                                                                                                                       >  1:3507595/1‑150 (MQ=255)
tGGCACCATGCGCAGCGCCGATGATATTACCGCTATGTGGAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTc                                                                                                                                                       >  1:336827/1‑150 (MQ=255)
    aCCATGCGCAGCGCCGATGATATTACCGCTATGTGGAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCgtg                                                                                                                                                   >  1:116906/1‑150 (MQ=255)
     ccATGCGCAGCGCCGATGATATTACCGCTATGTGGAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCgtgt                                                                                                                                                  >  1:3294792/1‑150 (MQ=255)
      cATGCGCAGCGCCGATGATATTACCGCTATGTGGAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTa                                                                                                                                                 >  1:2213611/1‑150 (MQ=255)
       aTGCGCAGCGCCGATGATATTACCGCTATGTGGAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTAt                                                                                                                                                >  2:2256750/1‑150 (MQ=255)
       aTGCGCAGCGCCGATGATATTACCGCTATGTGGAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTAt                                                                                                                                                <  1:4116280/150‑1 (MQ=255)
             aGCGCCGATGATATTACCGCTATGTGGAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGAcggc                                                                                                                                          >  1:949843/1‑150 (MQ=255)
             aGCGCCGATGATATTACCGCTATGTGGAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGAcggc                                                                                                                                          >  2:138749/1‑150 (MQ=255)
                 ccGATGATATTACCGCTATGTGGAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCt                                                                                                                                      >  2:3580615/1‑150 (MQ=255)
                  cGATGATATTACCGCTATGTGGAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTg                                                                                                                                     >  2:1597548/1‑150 (MQ=255)
                   gatgatATTACCGCTATGTGGAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTTCGTGTATGACGGCGGCTgg                                                                                                                                    <  1:1526546/150‑1 (MQ=255)
                   gatgatATTACCGCTATGTGGAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTTCGTGTATGACGGCGGCTgg                                                                                                                                    <  2:1143847/150‑1 (MQ=255)
                   gatgatATTACCGCTATGTGGAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTTCGTGTATGACGGCGGCTgg                                                                                                                                    <  1:3772112/150‑1 (MQ=255)
                   gatgatATTACCGCTATGTGGAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTgg                                                                                                                                    >  2:2935634/1‑150 (MQ=255)
                   gatgatATTACCGCTATGTGGAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTgg                                                                                                                                    >  2:2948960/1‑150 (MQ=255)
                    atgatATTACCGCTATGTGGAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGt                                                                                                                                   <  1:2683540/150‑1 (MQ=255)
                           tACCGCTATGTGGAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATg                                                                                                                            >  2:3447078/1‑150 (MQ=255)
                             ccGCTATGTGGAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGa                                                                                                                          >  2:2397848/1‑150 (MQ=255)
                             ccGCTATGTGGAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGa                                                                                                                          >  1:2087369/1‑150 (MQ=255)
                              cGCTATGTGGAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGag                                                                                                                         <  2:434198/150‑1 (MQ=255)
                              cGCTATGTGGAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGag                                                                                                                         <  2:3039905/150‑1 (MQ=255)
                               gCTATGTGGAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGagc                                                                                                                        >  1:841444/1‑150 (MQ=255)
                                cTATGTGGAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGagca                                                                                                                       >  1:4983779/1‑150 (MQ=255)
                                cTATGTGGAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGagca                                                                                                                       >  1:4983770/1‑150 (MQ=255)
                                 tATGTGGAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTCCGAATGGagcag                                                                                                                      >  2:500999/1‑150 (MQ=255)
                                 tATGTGGAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGagcag                                                                                                                      >  1:730771/1‑150 (MQ=255)
                                 tATGTGGAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGagcag                                                                                                                      <  1:4614827/150‑1 (MQ=255)
                                  aTGTGGAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGagcagc                                                                                                                     >  2:2955531/1‑150 (MQ=255)
                                   tgtgGAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCAGTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCg                                                                                                                    <  2:4422565/150‑1 (MQ=255)
                                    gtgGAAAGCGTGTAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGa                                                                                                                   >  1:1783525/1‑150 (MQ=255)
                                    gtgGAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGa                                                                                                                   >  2:1964699/1‑150 (MQ=255)
                                       gAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATcc                                                                                                                >  2:4687226/1‑150 (MQ=255)
                                        aaaGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCa                                                                                                               >  2:798544/1‑150 (MQ=255)
                                          aGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCaaa                                                                                                             >  1:1616102/1‑150 (MQ=255)
                                            cGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCaaaaa                                                                                                           >  1:3393406/1‑150 (MQ=255)
                                             gTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCaaaaaa                                                                                                          >  1:4375958/1‑150 (MQ=255)
                                               ggAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATc                                                                                                        >  2:1719643/1‑150 (MQ=255)
                                                  atatCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCgg                                                                                                     >  1:3034671/1‑150 (MQ=255)
                                                  atatCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCgg                                                                                                     <  1:2855430/150‑1 (MQ=255)
                                                      cAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGc                                                                                                 >  1:74615/1‑150 (MQ=255)
                                                       aaaCCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCa                                                                                                >  1:4318710/1‑150 (MQ=255)
                                                       aaaCCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCa                                                                                                <  2:4502597/150‑1 (MQ=255)
                                                        aaCCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCaa                                                                                               >  2:4972464/1‑150 (MQ=255)
                                                         aCCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAAc                                                                                              <  2:2435788/150‑1 (MQ=255)
                                                            agagCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCgg                                                                                           >  1:618154/1‑150 (MQ=255)
                                                            agagCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCgg                                                                                           <  1:1351397/150‑1 (MQ=255)
                                                             gagCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGt                                                                                          >  1:2822425/1‑150 (MQ=255)
                                                             gagCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGt                                                                                          <  1:1954837/150‑1 (MQ=255)
                                                               gcagcaAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGa                                                                                        >  2:5118003/1‑150 (MQ=255)
                                                               gcagaaAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGa                                                                                        <  1:1108518/150‑1 (MQ=255)
                                                                cagcaAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGaa                                                                                       >  1:2979109/1‑150 (MQ=255)
                                                                cagcaAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGaa                                                                                       <  1:2716873/150‑1 (MQ=255)
                                                                cagcaAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGaa                                                                                       <  1:4761808/150‑1 (MQ=255)
                                                                   caAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCACCCGGTGAAcgc                                                                                    >  2:4084707/1‑150 (MQ=255)
                                                                   caAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAAcgc                                                                                    <  1:1791920/150‑1 (MQ=255)
                                                                    aaGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAAcgcg                                                                                   <  1:1457130/150‑1 (MQ=255)
                                                                        ttCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGccc                                                                               >  2:3345561/1‑150 (MQ=255)
                                                                        ttCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGccc                                                                               <  2:3819688/150‑1 (MQ=255)
                                                                          cATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCgg                                                                             >  1:4428616/1‑150 (MQ=255)
                                                                           aTTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGa                                                                            <  1:5183103/150‑1 (MQ=255)
                                                                           aTTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGa                                                                            <  2:4375927/150‑1 (MQ=255)
                                                                           aTTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGa                                                                            <  1:3171282/150‑1 (MQ=255)
                                                                             tCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACa                                                                          >  1:651839/1‑150 (MQ=255)
                                                                             tCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCATTATGTACGCACGCGACATGGTTTGGAAGAATGTTTTCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACa                                                                          <  2:3942319/150‑1 (MQ=255)
                                                                              cTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAg                                                                         <  1:1737801/150‑1 (MQ=255)
                                                                                    cGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCaa                                                                   <  2:3911265/150‑1 (MQ=255)
                                                                                       cACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATa                                                                >  2:4878450/1‑150 (MQ=255)
                                                                                          cGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTTCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACa                                                             <  2:3383949/150‑1 (MQ=255)
                                                                                          cGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACa                                                             <  2:664316/150‑1 (MQ=255)
                                                                                           ggCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACAt                                                            <  1:3614058/150‑1 (MQ=255)
                                                                                           ggCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACAt                                                            <  1:3456542/150‑1 (MQ=255)
                                                                                           ggCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACAt                                                            >  1:1574377/1‑150 (MQ=255)
                                                                                               ggCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCa                                                        >  1:2916392/1‑150 (MQ=255)
                                                                                                 cgcgcgTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAAt                                                      >  1:2545169/1‑150 (MQ=255)
                                                                                                 cgcgcgTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAAt                                                      <  2:1922366/150‑1 (MQ=255)
                                                                                                    gcgTCCGAAACATTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCTATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTAc                                                   <  2:1086745/150‑1 (MQ=255)
                                                                                                     cgTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTTCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATCACATTGCAATTACt                                                  >  2:4196424/1‑150 (MQ=255)
                                                                                                       tCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGa                                                <  1:3271750/150‑1 (MQ=255)
                                                                                                        ccGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGAc                                               >  1:5121729/1‑150 (MQ=255)
                                                                                                         cGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACg                                              >  1:2556817/1‑150 (MQ=255)
                                                                                                           aaaCCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCt                                            >  2:576384/1‑150 (MQ=255)
                                                                                                           aaaCCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCt                                            >  1:3759871/1‑150 (MQ=255)
                                                                                                                tttATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGc                                       <  1:4496549/150‑1 (MQ=255)
                                                                                                                tttATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGc                                       >  2:2058443/1‑150 (MQ=255)
                                                                                                                    tGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACt                                   <  2:4446492/150‑1 (MQ=255)
                                                                                                                       acgcacgcGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCa                                >  1:2017773/1‑150 (MQ=255)
                                                                                                                       acgcacgcGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCa                                <  1:238020/150‑1 (MQ=255)
                                                                                                                        cgcacgcGCCATGTGTTGGAAGAATGTTTTCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGGAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAg                               <  2:1211529/150‑1 (MQ=255)
                                                                                                                               gcCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGa                        <  2:2625562/150‑1 (MQ=255)
                                                                                                                                ccATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGAt                       <  1:398550/150‑1 (MQ=255)
                                                                                                                                      gTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCAc                 >  1:4602990/1‑150 (MQ=255)
                                                                                                                                       ttGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACt                >  1:261164/1‑150 (MQ=255)
                                                                                                                                          gaagaaTGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCca             >  2:3475897/1‑150 (MQ=255)
                                                                                                                                           aagaaTGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCcac            <  1:1527061/150‑1 (MQ=255)
                                                                                                                                           aagaaTGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCcac            >  2:2000391/1‑150 (MQ=255)
                                                                                                                                            agaaTGTTTTCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCcaca           <  1:2115076/150‑1 (MQ=255)
                                                                                                                                                 gTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCcgc      <  2:4902629/150‑1 (MQ=255)
                                                                                                                                                  ttttCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCcgcg     <  2:4104515/150‑1 (MQ=255)
                                                                                                                                                    tCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTa   >  1:3673156/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGCACCATGCGCAGCGCCGATGATATTACCGCTATGTGGAAAGCGTGGAATATCAAACCAGAGCAGCAAGTTTCATTCTACTGCGGCACCGGCTGGCGCGCGTCCGAAACCTTTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAA  >  CP009273/1834794‑1835092

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: