Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,774,215 A→G 5.6% intergenic (+21/‑209) mtlA → / → mtlD mannitol‑specific PTS enzyme: IIA, IIB and IIC components/mannitol‑1‑phosphate dehydrogenase, NAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,774,2150AG 5.6% 4718.8 / 165.4 1414intergenic (+21/‑209)mtlA/mtlDmannitol‑specific PTS enzyme: IIA, IIB and IIC components/mannitol‑1‑phosphate dehydrogenase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base A (734/595);  new base G (52/27);  total (788/626)
Fisher's exact test for biased strand distribution p-value = 7.98e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.69e-02

Only 100 of 1449 total aligned reads displayed.

CGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTGCGCTTCGGTGAAGAAGAAGATGACATTGCCCGTCTGGTGATTGGTATTGCTGCCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTAAAAACATTGATGAAGGTTAATACTATGAAAGCATTACATTTTGGCGCAGGTAATATCGGTCGTGGCTTTATCGGTAAACTGCTGGCA  >  NC_000913/3773942‑3774486
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                               
             gCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTGCGCTTCGGCGAAGAAGAAGATGACATTGCCCGTCTGGTGATTGGTATTGCTGCCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGACGCTCTGGCACACACCACCAGCGGGGACGAAGTGCTGGAACTGCTGGCAGGCCGTAAGTCAACCAATCCCACCCTCCCCACATGGAGACgggggg                                                                                                                                                                                                                                                                    >  4:184864/1‑274 (MQ=255)
                                                             gAAGGCGTGCGCTTCGGTGAAGAAGAAGATGACATTGCCCGTCTGGTGATTGGTATTGCTGCCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGGCCGTCATCGAGCGTCTGGCACCCACCACCAGCGTGGACGAAGTGCTGGAACTGCTGGCCGGTCGTAAGTAATCCAAGCCCACCCCCCCAAAATGGGGAAGGTGGGGCTAATTGCCTGAGTCCCTCCGCTTTACAGGACTACAGGGTACATcac                                                                                                                                                                                                                    >  4:355373/1‑274 (MQ=255)
                                                                    ggagcTTCGGTGAAGAATAAGATGACCTGGCCCGTCTGGTGAGTGCTATTGCTGCCCTTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAAGGCACTGGATGATGAGTCCGGCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGTAACGGCTTGCACGTCGTAAGTAATCCAATCCCACCCTCTCCCCATGGAGAAGGt                                                                                                                                                                                                                                                                       <  1:19499/213‑1 (MQ=255)
                                                                                        gaTGACATTGCCAGTCTGGTGATTAGTATTGCTGCCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCTAGCCTCTGGCAACCACCACCAGCGTGGTGGAAGTGCACGAACTGCTGGCAGGTCGTAAGTACTCCAATCCCACCATCCCCACATGGAGAAGGTGGGGTTAATTGCCtca                                                                                                                                                                                                                                                       >  4:559080/1‑210 (MQ=255)
                                                                                                                              tAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACAt                                                                                                                                                                                                                                                                                <  6:523300/149‑1 (MQ=255)
                                                                                                                              tAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACAt                                                                                                                                                                                                                                                                                <  5:1300414/149‑1 (MQ=255)
                                                                                                                                acaacGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATgg                                                                                                                                                                                                                                                                              <  7:803922/149‑1 (MQ=255)
                                                                                                                                acaacGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCa                                                                                                                                                                                                   <  2:546764/224‑1 (MQ=255)
                                                                                                                                  aacGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAAACCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCtctc                                                                                                                                              <  2:332202/275‑1 (MQ=255)
                                                                                                                                  aacGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGACATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGgag                                                                                                                                                                                                                                                                            >  8:1217182/1‑149 (MQ=255)
                                                                                                                                     gAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAg                                                                                                                                                                                                                                                                         >  5:22082/1‑149 (MQ=255)
                                                                                                                                      aGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAgg                                                                                                                                                                                                                                                                        >  7:88290/1‑149 (MQ=255)
                                                                                                                                        tacaTTCAGTTTATCACCAGCCTGACCGATGCACTGGAGGATGAGTCCGCCATCGAGCGTCTGGCACACACCACCAGCGTGGATGACGTGCTGGAACTGCTGTCATGTCGTAAGTAATCCAATCCCACCCACTCCATGTGGAGAAGGtt                                                                                                                                                                                                                                                                      <  6:647819/148‑2 (MQ=255)
                                                                                                                                        cacaTTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCTTCTCCATGTGGAGAGGGTg                                                                                                                                                                                                                                                                      <  5:1718091/149‑1 (MQ=255)
                                                                                                                                        cacaTTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCTTCTCCATGTGGAGAGGGTg                                                                                                                                                                                                                                                                      <  6:1424840/149‑1 (MQ=255)
                                                                                                                                         acaTTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTgg                                                                                                                                                                                                                                                                     <  5:865948/149‑1 (MQ=255)
                                                                                                                                         acaTTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTgg                                                                                                                                                                                                                                                                     <  7:1395923/149‑1 (MQ=255)
                                                                                                                                          caTTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTggg                                                                                                                                                                                                                                                                    <  5:1034257/149‑1 (MQ=255)
                                                                                                                                           aTTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTgggg                                                                                                                                                                                                                                                                   <  6:711597/149‑1 (MQ=255)
                                                                                                                                            ttCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGt                                                                                                                                                                                                                                                                  <  7:172121/149‑1 (MQ=255)
                                                                                                                                                 gTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAAtt                                                                                                                                                                                                                                                             >  8:758449/1‑149 (MQ=255)
                                                                                                                                                 attATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGGACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAAtt                                                                                                                                                                                                                                                             <  7:309588/148‑1 (MQ=255)
                                                                                                                                                       aCCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGa                                                                                                                                                                                                                                                       >  6:253660/1‑149 (MQ=255)
                                                                                                                                                        ccAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGAt                                                                                                                                                                                                                                                      <  7:1533707/149‑1 (MQ=255)
                                                                                                                                                        ccAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGAt                                                                                                                                                                                                                                                      >  6:910491/1‑149 (MQ=255)
                                                                                                                                                         cAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATg                                                                                                                                                                                                                                                     <  6:189382/149‑1 (MQ=255)
                                                                                                                                                         cAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGa                                                                                                                                                                                                                            <  2:85621/174‑1 (MQ=255)
                                                                                                                                                                  cAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCCCCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGc                                                                                                                                                                                                                                             <  7:1459905/148‑1 (MQ=255)
                                                                                                                                                                  cAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCt                                                                                                                                                                                                                                            >  5:588766/1‑149 (MQ=255)
                                                                                                                                                                  cAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCt                                                                                                                                                                                                                                            >  7:955718/1‑149 (MQ=255)
                                                                                                                                                                     tGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTAt                                                                                                                                                                                                                                         <  7:382459/149‑1 (MQ=255)
                                                                                                                                                                          tGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCTCACCCTCTCCACATGGAGATGGTGGGGGTAATTGCCTGCTGCGCTACNCTTATCCAGCCTACCGGCTGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGTCCGGATAAGGCGCTTCCTCCTCTTCNTGCCCTGCCCTTCCCCCCCCGGCGGNGGTGTGGGTGGGGGAAAAGCTTccccccccccc                                                                                                  >  3:174641/1‑274 (MQ=255)
                                                                                                                                                                           ggATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAg                                                                                                                                                                                                                                                                         <  5:656338/111‑1 (MQ=255)
                                                                                                                                                                              tgatgaGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTAc                                                                                                                                                                                                                                <  8:1546770/149‑1 (MQ=255)
                                                                                                                                                                               gatgaGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACa                                                                                                                                                                                                                               >  8:293663/1‑149 (MQ=255)
                                                                                                                                                                                atgaGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAg                                                                                                                                                                                                                              <  7:724919/149‑1 (MQ=255)
                                                                                                                                                                                     tCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGc                                                                                                                                                                                                                         <  8:881279/149‑1 (MQ=255)
                                                                                                                                                                                      ccGTCATCGAGCGTCTGGCACACCCCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAATTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTTATACGCTACGCTTATAAGGCCTAAAGGATGCATAACAATTTGTTGAATTTGCAGGTTCTTGTAGACCGGAAAAGGCCCTTACGCCGAATCCGGCGCTTCCCCTCGCCCGCAGTAGCGGGTTTTGGTGAGGGAAAGACCTCACCACCGCCCCCGCGGGt                                                                                          >  4:49696/1‑275 (MQ=255)
                                                                                                                                                                                        gTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATc                                                                                                                                                                                                                      <  2:708714/149‑1 (MQ=255)
                                                                                                                                                                                         tCATCGAGCGTCTGGCACACACCACCCGCGTGGATGAAGTGCTGTAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTaa                                                                                        <  2:221668/274‑1 (MQ=255)
                                                                                                                                                                                         tCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATca                                                                                                                                                                                                                     <  5:738575/149‑1 (MQ=255)
                                                                                                                                                                                          cATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGt                                                                                                                                                                                                                                                                       <  2:280345/98‑1 (MQ=255)
                                                                                                                                                                                           aTCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATcaca                                                                                                                                                                                                                   >  7:1034467/1‑149 (MQ=255)
                                                                                                                                                                                             cGAGCGTTGGGCACACACCCCCAGCGTGGGTTAAGGGGTGGAACGGTTGGCAGGTCGTAAGTAATCCACTCCCACCCGCTCCACATGGGGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTAGCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCTCCCCAGCCCTCTCGGGTAAAAAc                                                                                    <  4:41791/274‑1 (MQ=255)
                                                                                                                                                                                                 cGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATttgt                                                                                                                                                                                                             <  6:1106028/149‑1 (MQ=255)
                                                                                                                                                                                                       gCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAAtt                                                                                                                                                                                                       >  6:1334545/1‑149 (MQ=255)
                                                                                                                                                                                                       gCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAAtt                                                                                                                                                                                                       >  5:930473/1‑149 (MQ=255)
                                                                                                                                                                                                         acacacCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTg                                                                                                                                                                                                     >  8:385185/1‑149 (MQ=255)
                                                                                                                                                                                                             accaccAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAgg                                                                                                                                                                                                                             <  4:298228/121‑1 (MQ=255)
                                                                                                                                                                                                             accaccAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGTGTGAGGGAAAAGCCTCACCCCACCCCTCTCGGGCAAAACCATTGATGAAGGTTAATATTa                                                                >  1:289011/1‑278 (MQ=255)
                                                                                                                                                                                                                 ccaGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCt                                                                                                                                                                                             >  8:720923/1‑149 (MQ=255)
                                                                                                                                                                                                                      gTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAg                                                                                                                                                                                        >  6:431220/1‑149 (MQ=255)
                                                                                                                                                                                                                         gATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATttgt                                                                                                                                                                                                             <  5:1394800/125‑1 (MQ=255)
                                                                                                                                                                                                                         gATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGcc                                                                                                                                                                                     >  8:1467266/1‑149 (MQ=255)
                                                                                                                                                                                                                         gATGAAGTGCTGGAACTGCCGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAAt                                                                                                                                                                                                                 >  8:1148440/1‑121 (MQ=255)
                                                                                                                                                                                                                          aTGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCg                                                                                                                                                                                    >  5:1455510/1‑149 (MQ=255)
                                                                                                                                                                                                                           tGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGaa                                                                                                                   <  2:103940/213‑1 (MQ=255)
                                                                                                                                                                                                                            gAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCgg                                                                                                                                                                                   >  5:502287/1‑148 (MQ=255)
                                                                                                                                                                                                                             aaGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTATAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGAt                                                                                                                                                                                 <  8:1376285/149‑1 (MQ=255)
                                                                                                                                                                                                                             aaGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGAt                                                                                                                                                                                 <  7:1331556/149‑1 (MQ=255)
                                                                                                                                                                                                                               gTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGACCGGATaa                                                                                                                                                                               <  5:42806/149‑1 (MQ=255)
                                                                                                                                                                                                                                tGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAg                                                                                                                                                                              <  7:1595986/149‑1 (MQ=255)
                                                                                                                                                                                                                                tGCTGGAACTGCTGGCAGGTCGTAACTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAg                                                                                                                                                                              <  7:1102531/149‑1 (MQ=255)
                                                                                                                                                                                                                                 gCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCAt                                                                                                                                                                                                                       <  7:428713/107‑1 (MQ=255)
                                                                                                                                                                                                                                    ggAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGgcgc                                                                                                                                                                          <  8:1281508/149‑1 (MQ=255)
                                                                                                                                                                                                                                      aaCTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATttgtt                                                                                                                                                                                                            >  7:1253798/1‑113 (MQ=255)
                                                                                                                                                                                                                                         tgctgGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTAcg                                                                                                                                                                     >  6:1539550/1‑149 (MQ=255)
                                                                                                                                                                                                                                          gctgGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACg                                                                                                                                                                                                 >  3:61417/1‑120 (MQ=255)
                                                                                                                                                                                                                                          gctgGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTAcgc                                                                                                                                                                    <  5:333362/149‑1 (MQ=255)
                                                                                                                                                                                                                                           ctgGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAg                                                                                                                                                                              <  7:1175434/138‑1 (MQ=255)
                                                                                                                                                                                                                                           ctgGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTAcgcc                                                                                                                                                                   <  6:977364/149‑1 (MQ=255)
                                                                                                                                                                                                                                              gCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCa                                                                                                                                                                >  5:792551/1‑149 (MQ=255)
                                                                                                                                                                                                                                               cAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATaa                                                                                                                                                                               <  6:1203840/133‑1 (MQ=255)
                                                                                                                                                                                                                                               cAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCAt                                                                                                                                                               <  7:225879/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     gtaagtaaTCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGgcg                                                                                                                                                         >  8:442835/1‑149 (MQ=255)
                                                                                                                                                                                                                                                        aggcaattaacCCCACCTTCTCCATGTGGAGAGGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTg                                                                                                                                                      >  5:16245‑M2/12‑149 (MQ=255)
                                                                                                                                                                                                                                                             tccaatccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCt                                                                                                                                                                                                                                                         <  5:996889/45‑1 (MQ=255)
                                                                                                                                                                                                                                                             tccaatccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGt                                                                                                                                                                                          >  8:194510/1‑108 (MQ=255)
                                                                                                                                                                                                                                                             tccaatccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATcc                                                                                                                                                             >  2:488398/1‑137 (MQ=255)
                                                                                                                                                                                                                                                             tccaatccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCct                                                                                                                                                 >  7:623027/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               caatccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCa                                                                                                                                                                                                                                       >  8:1070658/1‑61 (MQ=255)
                                                                                                                                                                                                                                                                 atccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTctcc                                                                                                                                             >  6:595829/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  tccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGACCCTctccg                                                                                                                                            >  7:538133/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                   cccACCTTCTCCATGTGGAGAGGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATc                                                                                                                                                              <  6:1634618/130‑1 (MQ=25)
                                                                                                                                                                                                                                                                    ccACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCa                                                                                                                                          >  8:530850/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                       cccTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTg                                                                                                                                                                                                     >  4:770532/1‑87 (MQ=255)
                                                                                                                                                                                                                                                                        cttctcCATGTGGAGAGGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAAt                                                                                                                                                                                                                 <  6:806128/72‑1 (MQ=21)
                                                                                                                                                                                                                                                                          tctcCACATGGAGAAGGTGTGGTTAATTGCCTGATGCGCTACGCTTATCAGGACTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGA                                                                                                                                    >  8:1430666‑M1/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                          tctcCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTAAAAACATTGATGAAGGTTAATACTATGAAAGCATTACATTTTGGCGCAGGTAATATCGGTCGTGGCTTTATCGGTAAACTGCTGGCa  <  1:106586/279‑1 (MQ=255)
                                                                                                                                                                                                                                                                          tctcCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGA                                                                                                                                    >  5:1527480‑M1/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                          tctcCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGA                                                                                                                                    <  6:337159‑M1/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                          tctcCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGA                                                                                                                                    <  8:625983‑M1/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                          tctcCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCCGCGCTGCCCCTCTCCTCACGGAGA                                                                                                                                    >  6:1255911‑M1/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                            tcCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGTCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCTCAGCCCGCCCGGGTAAAAACATTGATGAAGGTTAATACTAGGAAAGCATTACATTTTGGCGCAGGTAATATCGGTCGTGGCTTTATCGGTAAActgctg     >  2:409659/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                            tcCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGTCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGACAAGCCTCACCCCAGCCCTCTCGGGTAAAAACATTGATGCAGGTTAATACTATGAACTCATTACATTTTGGTGCAGGTATTACCGCGCGTGGCTTTACCGGTAAActtctg     >  4:657234/1‑274 (MQ=39)
                                                                                                                                                                                                                                                                            tcCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTAAAAACATTGATGAAGGTTAATACTATGAAAGCATTACATTTTGGCGCAGGTAATATCGGTCGTGGCTTTATCGGTAAActgctg     >  2:616585/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                            tcCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGG                                                                                                                                  >  8:1624961‑M1/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                            gccACATGGAGAAGTTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTAcgcc                                                                                                                                                                   <  1:65467/115‑1 (MQ=255)
                                                                                                                                                                                                                                                                              cacaTGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGc                                                                                                                                >  5:1126846‑M1/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                               acaTGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAgg                                                                                                                                                                                       <  5:1204367/93‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                               
CGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTGCGCTTCGGTGAAGAAGAAGATGACATTGCCCGTCTGGTGATTGGTATTGCTGCCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTAAAAACATTGATGAAGGTTAATACTATGAAAGCATTACATTTTGGCGCAGGTAATATCGGTCGTGGCTTTATCGGTAAACTGCTGGCA  >  NC_000913/3773942‑3774486

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