Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,732,141 G→C 6.9% intergenic (+18/‑60) xylF → / → xylG xylose ABC transporter periplasmic binding protein/xylose ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,732,1410GC 6.9% 308.8 / 9.3 96intergenic (+18/‑60)xylF/xylGxylose ABC transporter periplasmic binding protein/xylose ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base G (40/48);  new base C (4/3);  total (44/52)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

AGATACCACACTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGTTGAATGCTGG  >  NC_000913/3731997‑3732289
                                                                                                                                                |                                                                                                                                                    
aGATACCACACTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCg                                                                                                                                                 >  1:1224277/1‑150 (MQ=255)
aGATACCACACTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCg                                                                                                                                                 <  1:1310050/150‑1 (MQ=255)
 gATACCACACTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCgg                                                                                                                                                >  1:385308/1‑150 (MQ=255)
   tACCACACTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCgggg                                                                                                                                              <  1:946750/150‑1 (MQ=255)
   tACCACACTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCgggg                                                                                                                                              >  1:1527618/1‑150 (MQ=255)
     ccACACTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGc                                                                                                                                            >  2:1107331/1‑150 (MQ=255)
      cacaCTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCTAGCTGTAAGCGTTACGCTCCAGTGCTGTGCgggggtg                                                                                                                                           <  2:390963/150‑3 (MQ=255)
      cacaCTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCCGCTCCGCGCTGGGGCg                                                                                                                                           <  1:127408/150‑1 (MQ=255)
      cacaCTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCCGCTCCGCGCTGGGGCg                                                                                                                                           <  1:306549/150‑1 (MQ=255)
         acTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCCGCTCCGCGCtgggg                                                                                                                                             <  1:788627/145‑1 (MQ=255)
         acTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCCGCTCCGCGCtgggg                                                                                                                                             >  2:788627/1‑145 (MQ=255)
          cTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGCCGGATTCCACAAAGGGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGAt                                                                                                                                       <  1:1028483/150‑1 (MQ=255)
             aataatGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTtc                                                                                                                                    <  2:175992/150‑1 (MQ=255)
               taatGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTtctc                                                                                                                                  <  1:353700/150‑1 (MQ=255)
                   ggCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCa                                                                                                                              >  1:745806/1‑150 (MQ=255)
                   ggCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCa                                                                                                                              >  2:1176664/1‑150 (MQ=255)
                    gCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGa                                                                                                                                        >  1:923268/1‑139 (MQ=255)
                    gCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGa                                                                                                                                        <  2:923268/139‑1 (MQ=255)
                     ccTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATg                                                                                                                            >  2:1238732/1‑150 (MQ=255)
                     ccTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATg                                                                                                                            >  1:223392/1‑150 (MQ=255)
                     ccTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATg                                                                                                                            <  2:1269502/150‑1 (MQ=255)
                        gAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCcg                                                                                                                         >  2:427845/1‑150 (MQ=255)
                           aGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGt                                                                                                                      <  1:1477569/150‑1 (MQ=255)
                           aGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGt                                                                                                                      <  2:1665904/150‑1 (MQ=255)
                               gTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGtgaat                                                                                                                  >  1:897222/1‑150 (MQ=255)
                                tCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGtgaatg                                                                                                                 <  1:17576/150‑1 (MQ=255)
                                  cccTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGtgaatgaa                                                                                                               <  2:145966/150‑1 (MQ=255)
                                   ccTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCTGCGTgaatgaat                                                                                                              <  1:563229/150‑1 (MQ=255)
                                      cccGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTgg                                                                                                           <  2:745806/150‑1 (MQ=255)
                                       ccGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGc                                                                                                          <  1:1109071/150‑1 (MQ=255)
                                       ccGCCTCCTGACACAGATCGATGTTAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGAGCGGAGAGGGGGCTTGATTTCTATAAATGCAGCGTGAATGAATTGGc                                                                                                          >  2:838932/1‑150 (MQ=255)
                                         gCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCtt                                                                                                        >  1:746168/1‑150 (MQ=255)
                                          cctcctGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTa                                                                                                       <  1:196282/150‑1 (MQ=255)
                                          cctcctGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTa                                                                                                       <  2:1280771/150‑1 (MQ=255)
                                            tcctGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAgg                                                                                                     >  1:792510/1‑150 (MQ=255)
                                             cctGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGt                                                                                                    <  1:47268/150‑1 (MQ=255)
                                               tGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGCTAATTAAAGACGGATTCCACAAAGACAGCGAGCTGAAAGCGTTACGCCCCCGCGCCGAGCGGCGGCGTGATTTCTCTCCATGCCGCGTGAATGAAATGGCTTAGGtcc                                                                                                  >  2:633038/1‑148 (MQ=255)
                                                       atcgatGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTa                                                                                          <  2:968995/150‑1 (MQ=255)
                                                         cgatGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATg                                                                                        <  1:211370/150‑1 (MQ=255)
                                                           atGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGcc                                                                                      >  2:694489/1‑150 (MQ=255)
                                                           atGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGcc                                                                                      <  2:1527618/150‑1 (MQ=255)
                                                           atGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGcc                                                                                      <  2:164823/150‑1 (MQ=255)
                                                            tgtgAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGcc                                                                                      <  1:243290/149‑1 (MQ=255)
                                                            tgtgAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGcc                                                                                      >  2:243290/1‑149 (MQ=255)
                                                            tgtgAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCt                                                                                     <  2:385308/150‑1 (MQ=255)
                                                              tgAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTa                                                                                   <  2:68488/150‑1 (MQ=255)
                                                              tgAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTa                                                                                   <  2:245348/150‑1 (MQ=255)
                                                                aaTAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTtctctc                                                                                                                                >  1:584938/1‑103 (MQ=255)
                                                                aaTAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTtctctc                                                                                                                                <  2:584938/103‑1 (MQ=255)
                                                                aaTAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATc                                                                                 >  1:976278/1‑150 (MQ=255)
                                                                 aTAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGCGCGGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGTCGCGGCAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCt                                                                                <  2:1045501/150‑1 (MQ=255)
                                                                   aaaaaCAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTAc                                                                              >  2:120015/1‑150 (MQ=255)
                                                                     aaaCAACATCAACGATACCGTAATTAAACACGGAATCCACAAAGAGACCGAGCGGTAAGCGTTACGCCCCATCGAGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTTAATGAATTGGCTTAGGGTGAGTCGTTATCCCTTTTcttctt                                                                            <  2:1625533/150‑1 (MQ=255)
                                                                        caacaTCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTa                                                                               >  2:658482/1‑144 (MQ=255)
                                                                        caacaTCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTa                                                                               <  1:658482/144‑1 (MQ=255)
                                                                          acaTCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAAt                                                                       <  2:640881/150‑1 (MQ=255)
                                                                             tCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATgaa                                                                    >  1:1002342/1‑150 (MQ=255)
                                                                              cAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATgaag                                                                   >  1:1284352/1‑150 (MQ=255)
                                                                              cAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATgaag                                                                   >  1:396374/1‑150 (MQ=255)
                                                                              cAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATgaac                                                                   >  1:1290342/1‑149 (MQ=255)
                                                                               aaaGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATgaaga                                                                  >  1:309608/1‑150 (MQ=255)
                                                                                      cGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTAc                                                           <  1:1249838/150‑1 (MQ=255)
                                                                                       ggTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTAcc                                                          >  2:1209456/1‑150 (MQ=255)
                                                                                               aaGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTATCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTAGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCtt                                                  >  1:1153110/1‑150 (MQ=255)
                                                                                                 gACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCg                                                <  2:897222/150‑1 (MQ=255)
                                                                                                 gACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCg                                                >  2:312389/1‑150 (MQ=255)
                                                                                                  aCGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATTCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCgg                                               <  1:694489/150‑1 (MQ=255)
                                                                                                  aCGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCgg                                               <  1:427845/150‑1 (MQ=255)
                                                                                                     gATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCTCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAg                                            <  1:236809/150‑1 (MQ=255)
                                                                                                        tCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAgtgt                                         <  2:396374/150‑1 (MQ=255)
                                                                                                         ccACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGtgtg                                        >  1:826390/1‑150 (MQ=255)
                                                                                                          cacaAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGa                                       >  1:617504/1‑150 (MQ=255)
                                                                                                             aaaGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAgg                                    >  2:1271764/1‑150 (MQ=255)
                                                                                                              aaGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGc                                   <  2:903218/150‑1 (MQ=255)
                                                                                                                 agagCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGAt                                >  2:61408/1‑150 (MQ=255)
                                                                                                                    gcgagcTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGa                             >  2:1529307/1‑150 (MQ=255)
                                                                                                                    gcgagcTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGa                             >  2:697830/1‑150 (MQ=255)
                                                                                                                     cgagcTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGAt                            >  1:104638/1‑150 (MQ=255)
                                                                                                                      gagcTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATa                           <  2:1239939/150‑1 (MQ=255)
                                                                                                                      gagcTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATa                           <  2:534309/150‑1 (MQ=255)
                                                                                                                       agcTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATaa                          <  2:398362/150‑1 (MQ=255)
                                                                                                                           gTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTc                      <  2:104638/150‑1 (MQ=255)
                                                                                                                           gTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTc                      >  2:413497/1‑150 (MQ=255)
                                                                                                                               gCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCt                  <  2:223392/150‑1 (MQ=255)
                                                                                                                               aaatcaCGCCCCCGCTCCGCGCTGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTATGCt                  >  1:432590/6‑150 (MQ=255)
                                                                                                                                 atcaCGCCCCCGCTCCGCGCTGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTg                >  2:816570/4‑150 (MQ=255)
                                                                                                                                   tACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCg              <  2:527439/150‑1 (MQ=255)
                                                                                                                                    aCGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCgg             <  1:1209456/150‑1 (MQ=255)
                                                                                                                                    aCGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCgg             <  1:848331/150‑1 (MQ=255)
                                                                                                                                     cGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGt            >  2:417811/1‑150 (MQ=255)
                                                                                                                                     cGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGt            >  1:561069/1‑150 (MQ=255)
                                                                                                                                     cGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGt            >  1:834622/1‑150 (MQ=255)
                                                                                                                                       ccccAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGTTg          <  1:697830/150‑1 (MQ=255)
                                                                                                                                          cAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGTTGAAt       <  2:746168/150‑1 (MQ=255)
                                                                                                                                           aGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGt            >  1:1286274/1‑144 (MQ=255)
                                                                                                                                           aGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGt            <  2:1286274/144‑1 (MQ=255)
                                                                                                                                             cgcgGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGTTGAATGCt    >  2:1285055/1‑150 (MQ=255)
                                                                                                                                              gcgGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAgtg                                          >  2:764380/1‑111 (MQ=255)
                                                                                                                                              gcgGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAgtg                                          <  1:764380/111‑1 (MQ=255)
                                                                                                                                               cgGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGTTGAATGCTgg  <  1:120015/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
AGATACCACACTGAATAATGGCCTGAAAGATGTCCCCTCCCGCCTCCTGACACCGATCGATGTGAATAAAAACAACATCAAAGATACGGTAATTAAAGACGGATTCCACAAAGAGAGCGAGCTGTAAGCGTTACGCCCCAGCGCGGAGCGGGGGCGTGATTTCTCTCCATGCCGCGTGAATGAATTGGCTTAGGTGGAGTCGTTATGCCTTATCTACTTGAAATGAAGAACATTACCAAAACCTTCGGCAGTGTGAAGGCGATTGATAACGTCTGCTTGCGGTTGAATGCTGG  >  NC_000913/3731997‑3732289

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: