Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,938,548 A→T T361S (ACG→TCG)  fucK → L‑fuculokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,938,5480AT98.1% 172.1 / NA 53T361S (ACG→TCG) fucKL‑fuculokinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (25/27);  minor base G (0/1);  total (25/28)
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.97e-01

AAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCACGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGATTAAAGCCAATATGCTTGATATTCCGGTAAAAGTTCTCGACGACGCCGAAACGACCGTCGCAGGAGCTGCGCTGTTCGGTTGGTATGGCGTAGGGG  >  NC_000913/2938304‑2938791
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
aaGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCgg                                                                                                                                                                                                                                                <  2:4111/250‑1 (MQ=255)
   tCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGCAGCTTAATACCTCGCGGGGGc                                                                                                                                                                                                                                            >  2:211111/1‑251 (MQ=255)
     aGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCAt                                                                                                                                                                                                                                          <  2:139277/251‑1 (MQ=255)
             ggTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATcg                                                                                                                                                                                                                                  >  1:216108/1‑251 (MQ=255)
                tCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGAGACGCTTAATACCTCGCGGGGGCATTTCTATcgcgc                                                                                                                                                                                                                               >  2:320525/1‑251 (MQ=255)
                        cGAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTgg                                                                                                                                                                                                                       >  2:15886/1‑251 (MQ=255)
                         gAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTgg                                                                                                                                                                                                                       >  2:202526/1‑250 (MQ=255)
                          aaCTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCCTTTCTATCGCGCGGCGCTGGa                                                                                                                                                                                                                      >  2:319546/1‑250 (MQ=255)
                                      aggcaggGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCGGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTgcg                                                                                                                                                                                                         <  1:481398/251‑1 (MQ=255)
                                                      cccAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATc                                                                                                                                                                                         >  1:556780/1‑251 (MQ=255)
                                                       ccAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGTGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCt                                                                                                                                                                                        <  2:65343/251‑1 (MQ=255)
                                                                  aTGGCTGGCATCCGGCGTGCTGGTATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACACTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAACAATGCAGTGTGATTTATTGTCGTGTCAGAAAGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGGGCGGCGCTGGAAGGGTTAACTGCGCAATTAAAGCGACATCTACAGATGCTgg                                                                                                                                                                             >  2:357364/1‑251 (MQ=255)
                                                                          cATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATc                                                                                                                                                                     <  2:38233/251‑1 (MQ=255)
                                                                              cGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGc                                                                                                                                                                 <  1:103718/251‑1 (MQ=255)
                                                                                 cGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACt                                                                                                                                                              >  1:281022/1‑251 (MQ=255)
                                                                                         aaTGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCCATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGcc                                                                                                                                                      <  2:424680/251‑1 (MQ=255)
                                                                                          aTGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCGCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCct                                                                                                                                                     <  2:228118/251‑1 (MQ=255)
                                                                                              gTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAAGCGGGCACTTTAAGGCCTCTGa                                                                                                                                                 <  1:584842/251‑1 (MQ=255)
                                                                                                 agaAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCTCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAatt                                                                                                                                              <  1:367480/251‑1 (MQ=255)
                                                                                                         gTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAGAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGtt                                                                                                                                        >  2:278672/1‑249 (MQ=255)
                                                                                                              ggACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTc                                                                                                                                    <  1:292283/248‑1 (MQ=255)
                                                                                                              ggACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGt                                                                                                                                 >  1:524006/1‑251 (MQ=255)
                                                                                                                  ggCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTgg                                                                                                                               >  1:582440/1‑249 (MQ=255)
                                                                                                                        aaCACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTc                                                                                                                       <  2:574863/251‑1 (MQ=255)
                                                                                                                         acacCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTAGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCg                                                                                                                      <  2:18771/251‑1 (MQ=255)
                                                                                                                                 gCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACAt                                                                                                              <  1:97034/251‑1 (MQ=255)
                                                                                                                                    aaTGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATtgt                                                                                                           >  1:10675/1‑251 (MQ=255)
                                                                                                                                     aTGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTCCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATtgtg                                                                                                          >  2:401008/1‑251 (MQ=255)
                                                                                                                                       gTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGa                                                                                                        <  2:115366/251‑1 (MQ=255)
                                                                                                                                          gattgaAGAAGCTCGTCTGATCGCGCCTGTCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTTTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTTTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATc                                                                                                     <  1:347446/251‑1 (MQ=255)
                                                                                                                                          gattgaAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAAt                                                                                                      <  2:251909/250‑1 (MQ=255)
                                                                                                                                              gaagaaGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGa                                                                                                  <  2:119830/250‑1 (MQ=255)
                                                                                                                                                agaagCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGAtt                                                                                                >  1:22247/1‑250 (MQ=255)
                                                                                                                                                       cGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGATTAAAGCCaa                                                                                        <  2:429091/251‑1 (MQ=255)
                                                                                                                                                         gcTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGATTAAAGCCAata                                                                                      <  2:74687/250‑1 (MQ=255)
                                                                                                                                                             aTCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCTGCGCTGGAAGGGGTAACTGCGCAATTACAGCGCAAGCTACAGATGCTGGAAAAAATCGGGCACGTTAAGGCCTCTGAATTATGGTTAGTCGGTGGAGGAAGTCTCAGCGCATTGTGGAATCAGATTAAAGGGAATATGCt                                                                                  <  2:288849/251‑1 (MQ=255)
                                                                                                                                                                   ccTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGATTAAAGCCAATATGCTTGatat                                                                            >  1:354771/1‑251 (MQ=255)
                                                                                                                                                                        cgcgGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGATTAAAGCCAATATGCTTGATATTCCgg                                                                       <  2:110002/251‑1 (MQ=255)
                                                                                                                                                                                   tAAAACTGCAGTGTGATTTATTGTCGTGTCAGACCGATGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGATTAAAGCCAATATGCTTGATATTCCGGTAAAAGTtctc                                                            >  2:292392/1‑251 (MQ=255)
                                                                                                                                                                                        aTGCAGTGTGATTTATTGTCGTGTAAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATAGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGATTAAAGCCAATATGCTTGATATTCCGGTAAAAGTTCTcgacg                                                        >  2:326012/1‑250 (MQ=255)
                                                                                                                                                                                                        tGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGCAGTCGCAACACATTGTGGAATCAGATTAAAGCCAATATGCTTGATATTCCGGTAACAGTTCTCGACGACTCCGAAACGACCGTc                                       >  2:357782/1‑251 (MQ=255)
                                                                                                                                                                                                         gTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGATTAAAGCCAATATGCTTGATATTCCGGTAAAAGTTCTCGACGACGCCGAAACGACCGTc                                       <  1:197479/250‑1 (MQ=255)
                                                                                                                                                                                                          tCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGATTAAAGCCAATATGCTTGATATTCCGGTAAAAGTTCTCGACGACGCCGAAACGACCGTCGc                                     >  1:103237/1‑251 (MQ=255)
                                                                                                                                                                                                            gtgtCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCCTTTCTATCGCGAGGCGCTGGCAGGGTTAACTGCGAAATTACAGCGCAATCTACAGATGATGGAAAAAGTCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCCGATTCAAGCCAATATGCTTGATATTCCGGTAAAAGTTCTCGACGACGCCGAAACGCCCGTCGc                                     >  1:371397/1‑249 (MQ=255)
                                                                                                                                                                                                                cAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGATTAAAGCCAATATGCTTGATATTCCGGTAAAAGTTCTCGACGACGCCGAAACGACCGTCGCAGGAg                                <  1:280107/250‑1 (MQ=255)
                                                                                                                                                                                                                 aGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGATTAAAGCCAATATGCTTGATATTCCGGTAAAAGTTCTCGACGACGCCGAAACGACCGTCGCAGGAGCt                              <  1:202526/251‑1 (MQ=255)
                                                                                                                                                                                                                  gAACGCTGGCTGGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGATTAAAGCCAATATGCTTGATATTCCGGTAAAAGTTCTCGACGACGCCGAAACGACCGTCGCAGGAGCTg                             >  1:85075/1‑251 (MQ=255)
                                                                                                                                                                                                                      gctggctgGCAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGATTAAAGCCAATATGCTTGATATTCCGGTAAAAGTTCTCGACGACGCCGAAACGACCGTCGCAGGAGCTGCGCt                         <  1:139442/251‑1 (MQ=255)
                                                                                                                                                                                                                        ttgctggaAAGGAGTGACGCTTAATACCTCGGGGGAGCCTTTCTATCGCGCGGCGCTGGGAGGGTTCACTGAGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTAGCAACACATTGTGGAATCAGATTAAAGCCAATATGCTTGATATTCCGGTAAAAGTTCTCGACGACGCCGAAACGACCGTCGCAGGCGCTGCGCTGt                       >  2:58273/3‑251 (MQ=255)
                                                                                                                                                                                                                          gctggcAAGGAGTGACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGACGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGATTAAAGCCAATATGCTTGATATTCCGGTAAAAGTTCTCGACGACGCCGAAACGACCGTCGCAGGAGCTGCGCTGTTc                     >  1:359237/1‑251 (MQ=255)
                                                                                                                                                                                                                                       gACGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGATTAAAGCCAATATGCTTGATATTCCGGTAAAAGTTCTCGACGACGCCGAAACGACCGTCGCAGGAGCTGCGCTGTTCGGTTGGTATGGc         <  2:305172/250‑1 (MQ=255)
                                                                                                                                                                                                                                         cGCTTAATACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGATTAAAGCCAATATGCTTGATATTCCGGTAAAAGTTCTCGACGACGCCGAAACGACCGTCGCAGGAGCTGCGCTGTTCGGTTGGTATGGCGTa      <  1:295978/251‑1 (MQ=255)
                                                                                                                                                                                                                                              aaTACCTCGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGGAGTCGCAACACATTGTGGAATCAGATTAAAGCCAATATGCTTGATATTCCGGTAAAAGTTCTCGACGACGCCGAAACGACCGTCGCAGGAGCTGCGCTGTTCGGTTGGTATGGCGTAgggg  >  1:41167/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
AAGTCAGTACGCCGGTTCCACCTGCGAACTGGATAGCCAGGCAGGGTTGTATAACCCAGGTATGCAATGGCTGGCATCCGGCGTGCTGGAATGGGTGAGAAAACTGTTCTGGACGGCTGAAACACCCTGGCAAATGTTGATTGAAGAAGCTCGTCTGATCGCGCCTGGCGCGGATGGCGTAAAAATGCAGTGTGATTTATTGTCGTGTCAGAACGCTGGCTGGCAAGGAGTGACGCTTAATACCACGCGGGGGCATTTCTATCGCGCGGCGCTGGAAGGGTTAACTGCGCAATTACAGCGCAATCTACAGATGCTGGAAAAAATCGGGCACTTTAAGGCCTCTGAATTATTGTTAGTCGGTGGAGGAAGTCGCAACACATTGTGGAATCAGATTAAAGCCAATATGCTTGATATTCCGGTAAAAGTTCTCGACGACGCCGAAACGACCGTCGCAGGAGCTGCGCTGTTCGGTTGGTATGGCGTAGGGG  >  NC_000913/2938304‑2938791

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