Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 786,506 G→A 5.5% G292S (GGC→AGC)  aroG → 3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, phenylalanine repressible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913786,5060GA 5.5% 377.5 / 7.8 108G292S (GGC→AGC) aroG3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, phenylalanine repressible
Reads supporting (aligned to +/- strand):  ref base G (53/49);  new base A (3/3);  total (56/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 109 total aligned reads displayed.

CGCGAAGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACT  >  NC_000913/786358‑786654
                                                                                                                                                    |                                                                                                                                                    
cgcgAAGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTgg                                                                                                                                                     <  2:898552/150‑1 (MQ=255)
   gcaGCACGTTGCTGACGTGAAAGCAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCCGATGGATGTGTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGa                                                                                                                                                  <  1:600359/148‑1 (MQ=255)
       cACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAg                                                                                                                                              <  1:1070481/150‑1 (MQ=255)
          gTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGcc                                                                                                                                           >  1:1601567/1‑150 (MQ=255)
          gTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGcc                                                                                                                                           >  2:657285/1‑150 (MQ=255)
              cTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTAGCGAAAAGGCCatta                                                                                                                                       >  2:137328/1‑150 (MQ=255)
              cTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTAGCGAAAAGGCCatta                                                                                                                                       >  2:137315/1‑150 (MQ=255)
              cTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTAGCGAAAAGGCCatta                                                                                                                                       >  2:985372/1‑150 (MQ=255)
                gAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCattatt                                                                                                                                     <  1:1268983/150‑1 (MQ=255)
                 aaGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATTGATGTTTGTGCTGACGTTTGCCAGCCGATTGCCGGTGGCGAAAAGGCCATTATTg                                                                                                                                    <  1:1316419/150‑1 (MQ=255)
                 aaGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTg                                                                                                                                    <  2:55770/150‑1 (MQ=255)
                 aaGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTg                                                                                                                                    <  2:55768/150‑1 (MQ=255)
                 aaGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTg                                                                                                                                    <  1:258521/150‑1 (MQ=255)
                  aGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTgg                                                                                                                                   <  2:91828/150‑1 (MQ=255)
                   gTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGc                                                                                                                                  >  2:1172206/1‑150 (MQ=255)
                   gTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGc                                                                                                                                  >  2:1172209/1‑150 (MQ=255)
                   gTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGc                                                                                                                                  >  2:1577088/1‑150 (MQ=255)
                      aaaGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTg                                                                                                                               >  2:1354114/1‑150 (MQ=255)
                        agaagGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGAt                                                                                                                             <  1:103406/150‑1 (MQ=255)
                         gaagGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACTATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGAtg                                                                                                                            <  1:118651/150‑1 (MQ=255)
                           agGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGAtggt                                                                                                                          <  1:671910/150‑1 (MQ=255)
                             ggCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGtggtggtgg                                                                                                                        <  1:172970/150‑1 (MQ=255)
                             ggCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGtggtggtgg                                                                                                                        <  1:173052/150‑1 (MQ=255)
                             ggCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGAtggtgg                                                                                                                        >  2:1695513/1‑150 (MQ=255)
                             ggCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGAtggtgg                                                                                                                        >  2:1695521/1‑150 (MQ=255)
                             ggCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGAtggtgg                                                                                                                        >  2:570534/1‑150 (MQ=255)
                             ggCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGAtggtgg                                                                                                                        <  2:862051/150‑1 (MQ=255)
                                tGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCACAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGaaa                                                                                                                     <  2:147798/150‑1 (MQ=255)
                                tGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCTGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGaaa                                                                                                                     <  1:1035286/150‑1 (MQ=255)
                                tGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGaaa                                                                                                                     <  2:147800/150‑1 (MQ=255)
                                tGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGaaa                                                                                                                     <  2:939159/150‑1 (MQ=255)
                                tGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGaaa                                                                                                                     <  2:939162/150‑1 (MQ=255)
                                  aaCAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGc                                                                                                                   >  2:981683/1‑150 (MQ=255)
                                  aaCAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGc                                                                                                                   >  2:981679/1‑150 (MQ=255)
                                  aaCAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAAGTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGc                                                                                                                   >  2:981651/1‑150 (MQ=255)
                                  aaCAAAGCAGGCCTGCCAGCACACGAGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGGTTGTGCTGACGTTTGCCAGCAGATAGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGc                                                                                                                   >  2:981682/1‑150 (MQ=255)
                                    cAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCa                                                                                                                 >  1:1048022/1‑150 (MQ=255)
                                          aGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCtggt                                                                                                           >  1:434473/1‑150 (MQ=255)
                                            gCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCtggtgg                                                                                                         >  1:165030/1‑150 (MQ=255)
                                            gCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCtggtgg                                                                                                         >  2:1019976/1‑150 (MQ=255)
                                             ccTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGa                                                                                                        >  1:228588/1‑150 (MQ=255)
                                              cTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGaa                                                                                                       <  1:305860/150‑1 (MQ=255)
                                               tGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAg                                                                                                      >  2:903058/1‑150 (MQ=255)
                                               tGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAg                                                                                                      >  2:314414/1‑150 (MQ=255)
                                               tGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAg                                                                                                      <  1:1670259/150‑1 (MQ=255)
                                                   aGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCaa                                                                                                  >  2:68272/1‑150 (MQ=255)
                                                            gatgatCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCt                                                                                         <  1:29919/150‑1 (MQ=255)
                                                               gatCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCga                                                                                      <  1:1432819/150‑1 (MQ=255)
                                                                atCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCgag                                                                                     >  2:949400/1‑150 (MQ=255)
                                                                atCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCgag                                                                                     >  2:741542/1‑150 (MQ=255)
                                                                atCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCgag                                                                                     >  2:1397434/1‑150 (MQ=255)
                                                                atCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCgag                                                                                     >  2:949397/1‑150 (MQ=255)
                                                                   gATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGc                                                                                  >  2:974078/1‑150 (MQ=255)
                                                                         aGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAg                                                                            >  1:1460243/1‑150 (MQ=255)
                                                                         aGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAg                                                                            >  2:824211/1‑150 (MQ=255)
                                                                         aGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAg                                                                            >  1:748949/1‑150 (MQ=255)
                                                                              tGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCt                                                                       >  1:1600551/1‑150 (MQ=255)
                                                                              tGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCt                                                                       >  1:749090/1‑150 (MQ=255)
                                                                               gCTAACTCGTCCAAACAATTCAAACAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAACGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTg                                                                      <  2:1562861/150‑1 (MQ=255)
                                                                               gCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTg                                                                      <  2:1562832/150‑1 (MQ=255)
                                                                                cTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTgg                                                                     <  1:386618/150‑1 (MQ=255)
                                                                                cTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTgg                                                                     <  1:263727/150‑1 (MQ=255)
                                                                                 tAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGc                                                                    <  1:245542/150‑1 (MQ=255)
                                                                                    cTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTa                                                                 >  1:356904/1‑150 (MQ=255)
                                                                                     tCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTAc                                                                >  1:1220140/1‑150 (MQ=255)
                                                                                     tCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTAc                                                                >  2:23243/1‑150 (MQ=255)
                                                                                       gTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACgg                                                              <  2:1048022/150‑1 (MQ=255)
                                                                                        tCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGt                                                             >  1:705315/1‑150 (MQ=255)
                                                                                           aaaCAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAag                                                          >  1:1699408/1‑150 (MQ=255)
                                                                                           aaaCAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAag                                                          >  2:683552/1‑150 (MQ=255)
                                                                                             aCAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAagag                                                        >  2:676350/1‑150 (MQ=255)
                                                                                               aaTTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCa                                                      >  2:38568/1‑150 (MQ=255)
                                                                                               aaTTCAAAAAGCAGATGGAGGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTCAGAGCa                                                      >  2:1020527/1‑150 (MQ=255)
                                                                                                    aaaaaGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCAcc                                                 >  2:1065578/1‑150 (MQ=255)
                                                                                                     aaaaGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTAGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCg                                                <  1:137328/150‑1 (MQ=255)
                                                                                                     aaaaGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTAGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCg                                                <  1:137315/150‑1 (MQ=255)
                                                                                                     aaaaGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTAGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCg                                                <  1:985372/150‑1 (MQ=255)
                                                                                                         gCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGc                                            >  1:13467/1‑150 (MQ=255)
                                                                                                           aGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCt                                          <  1:1515306/150‑1 (MQ=255)
                                                                                                             atggatgTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGc                                        >  2:184667/1‑150 (MQ=255)
                                                                                                                gatgTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGGGCATCACCGATGCCTGCATc                                     >  2:845480/1‑150 (MQ=255)
                                                                                                                gatgTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATc                                     >  2:157184/1‑150 (MQ=255)
                                                                                                                    tttGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCt                                 >  2:684749/1‑150 (MQ=255)
                                                                                                                    tttGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCt                                 >  1:1534058/1‑150 (MQ=255)
                                                                                                                    tttGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCt                                 <  2:1143891/150‑1 (MQ=255)
                                                                                                                     ttGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTg                                <  1:560887/150‑1 (MQ=255)
                                                                                                                      tgtgCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCAGGCATCggtgat                               <  1:683552/150‑5 (MQ=255)
                                                                                                                         gCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGaa                            >  1:632631/1‑150 (MQ=255)
                                                                                                                         gCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGaa                            >  1:283440/1‑150 (MQ=255)
                                                                                                                         gCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGaa                            >  2:1138018/1‑150 (MQ=255)
                                                                                                                           tGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGa                          <  2:1601567/150‑1 (MQ=255)
                                                                                                                                          gATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGtt           <  1:1563548/150‑1 (MQ=255)
                                                                                                                                           aTTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTa          <  1:1569022/150‑1 (MQ=255)
                                                                                                                                             tGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACg        <  1:1397434/150‑1 (MQ=255)
                                                                                                                                              gCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGt       <  2:1709802/150‑1 (MQ=255)
                                                                                                                                               ccGGTGGAGAAAAGGCCATTATTGTCGTGATGGTGGATAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAATATACCGATGCTCTGTTACGTc      <  2:1716087/150‑1 (MQ=255)
                                                                                                                                                 ggTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCaa    <  1:974078/150‑1 (MQ=255)
                                                                                                                                                   tGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACt  <  1:890646/150‑1 (MQ=255)
                                                                                                                                                   tGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACt  <  1:890651/150‑1 (MQ=255)
                                                                                                                                                   tGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACt  <  1:924554/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGCGAAGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACT  >  NC_000913/786358‑786654

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: