Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,573,897 A→G 5.1% W137R (TGG→CGG)  eutT ← cobalamin adenosyltransferase involved in ethanolamine utilization

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,573,8970AG 5.1% 243.8 / 3.1 78W137R (TGG→CGG) eutTcobalamin adenosyltransferase involved in ethanolamine utilization
Reads supporting (aligned to +/- strand):  ref base A (41/33);  new base G (2/2);  total (43/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.80e-01

GTGGTCGAGATAGCGCAGCGGCTGATGAGAAAGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTG  >  NC_000913/2573754‑2574043
                                                                                                                                               |                                                                                                                                                  
gTGGTCGAGATAGCGCAGCGGCTGATGAGAAAGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGTGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGc                                                                                                                                                <  4:223901/148‑1 (MQ=255)
    tCGAGATAGCGCAGCGGCTGATGAGAAAGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCa                                                                                                                                           >  1:5585/1‑149 (MQ=255)
    tCGAGATAGCGCAGCGGCTGATGAGAAAGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCa                                                                                                                                           >  7:266023/1‑149 (MQ=255)
    tCGAGATAGCGCAGCGGCTGATGAGAAAGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCa                                                                                                                                           <  4:218500/149‑1 (MQ=255)
      gagaTAGCGCAGCGGCTGATGAGAAAGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCAc                                                                                                                                          <  7:106737/148‑1 (MQ=255)
            gcgcAGCGGCTGATGAGAAAGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTcc                                                                                                                                    <  2:293557/148‑1 (MQ=255)
              gcAGCGGCTGATGAGAAAGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGcc                                                                                                                                 <  2:11016/149‑1 (MQ=255)
               cAGCGGCTGATGAGAAAGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCa                                                                                                                                <  2:408986/149‑1 (MQ=255)
                  cGGCTGATGAGAAAGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGtt                                                                                                                             <  2:5585/149‑1 (MQ=255)
                  cGGCTGATGAGAAAGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGtt                                                                                                                             >  7:106041/1‑149 (MQ=255)
                      tgatgaGAAAGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCa                                                                                                                                <  3:77974/142‑1 (MQ=255)
                      tgatgaGAAAGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCa                                                                                                                                >  4:77974/1‑142 (MQ=255)
                       gatgaGAAAGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAAtt                                                                                                                        >  1:376408/1‑149 (MQ=255)
                       gatgaGAAAGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAAtt                                                                                                                        >  2:194766/1‑149 (MQ=255)
                           agaAAGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCa                                                                                                                    <  7:178275/149‑1 (MQ=255)
                             aaaGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGc                                                                                                                  >  2:339782/1‑149 (MQ=255)
                               aGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCca                                                                                                                <  2:376408/149‑1 (MQ=255)
                                      gCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCACCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGc                                                                                                          >  1:55723/1‑148 (MQ=255)
                                       cAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCGGCCGTGGCTggcgg                                                                                                                                          >  8:108874/1‑111 (MQ=255)
                                       cAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCCGCCGTGGCTggcgg                                                                                                                                          <  7:108874/115‑5 (MQ=255)
                                       cAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCCGCCGTGGCTggcgg                                                                                                                                          <  7:102991/115‑5 (MQ=255)
                                       cAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCGTGGCTggcgg                                                                                                                                          >  8:102991/1‑111 (MQ=255)
                                       cAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCa                                                                                                        <  4:86389/149‑1 (MQ=255)
                                       cAGAGCTTCGTCAGAAAGCCCGACACTCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCAGCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCa                                                                                                        <  2:342539/149‑1 (MQ=255)
                                         gATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAgc                                                                                                      >  1:379157/1‑149 (MQ=255)
                                                    gAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATg                                                                                                <  2:6807/144‑1 (MQ=255)
                                                    gAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATg                                                                                                >  1:6807/1‑144 (MQ=255)
                                                    gAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCt                                                                                           >  2:227203/1‑149 (MQ=255)
                                                     aaaGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTg                                                                                          <  3:288792/149‑1 (MQ=255)
                                                           cGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCgag                                                                                     >  7:130974/1‑148 (MQ=255)
                                                                 tCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACcgcg                                                                              >  5:70885/1‑149 (MQ=255)
                                                                  cGCCTGGCAGCCCCGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATAACCGACCGCCACGGCAGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCCGCGCAATGGTGCTGTCGAGAACcgcgc                                                                             >  1:30481/1‑149 (MQ=255)
                                                                  cGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACcgcgc                                                                             >  6:389740/1‑149 (MQ=255)
                                                                    ccTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGa                                                                           <  8:106041/149‑1 (MQ=255)
                                                                     cTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGa                                                                           >  4:355197/1‑148 (MQ=255)
                                                                       ggCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATTATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGGAATGGTGCTGGCGAGAACCGCGGGagaa                                                                        >  5:202142/1‑149 (MQ=255)
                                                                        gCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAAc                                                                       <  1:339782/149‑1 (MQ=255)
                                                                        gCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAAc                                                                       >  3:321857/1‑149 (MQ=255)
                                                                               aGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGAcg                                                                >  1:403026/1‑149 (MQ=255)
                                                                                  ggTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCgg                                                             <  5:346794/149‑1 (MQ=255)
                                                                                      cccccAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATc                                                          <  2:379157/148‑1 (MQ=255)
                                                                                        cccAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCt                                                       >  1:281670/1‑149 (MQ=255)
                                                                                         ccAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCtt                                                      >  2:179621/1‑149 (MQ=255)
                                                                                           aGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCtttt                                                    <  3:355197/149‑1 (MQ=255)
                                                                                                aTCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGa                                               <  2:403026/149‑1 (MQ=255)
                                                                                                aTCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGa                                               >  8:176248/1‑149 (MQ=255)
                                                                                                 tCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGAc                                              >  2:187237/1‑149 (MQ=255)
                                                                                                    gcgcgcATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCAt                                           >  1:48769/1‑149 (MQ=255)
                                                                                                      gcgcATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTAGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCAttt                                         >  8:266674/1‑149 (MQ=255)
                                                                                                        tcATAATGTTGCCAAGACGCGAACGGATATCAGCCAGCAACGGCTGCCACGGTTGCGCCAGATCAATTTGAAGCAACACTGCCAGCGCAAAGGTGCAGTCGAGAACCGAGCGAAAACCAAGACGCGGATCGCTTTTGGCTACCAtgttt                                       >  2:32760/2‑149 (MQ=255)
                                                                                                        gcATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCAttttt                                       >  7:137129/1‑149 (MQ=255)
                                                                                                         cATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACcgcgcg                                                                            <  4:178761/111‑1 (MQ=255)
                                                                                                         cATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACcgcgcg                                                                            >  3:178761/1‑111 (MQ=255)
                                                                                                          aTAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTc                                     <  8:130974/149‑1 (MQ=255)
                                                                                                          aTAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTc                                     >  6:306303/1‑149 (MQ=255)
                                                                                                           tAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGTCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGCCCATTTTTTcc                                    >  1:127029/1‑149 (MQ=255)
                                                                                                           tAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGTGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCAttt                                         >  4:8006/1‑144 (MQ=255)
                                                                                                           tAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGTGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCAttt                                         <  3:8006/144‑1 (MQ=255)
                                                                                                           tAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCttt                                                     >  1:81921/1‑132 (MQ=255)
                                                                                                           tAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCttt                                                     <  2:81921/132‑1 (MQ=255)
                                                                                                           tAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTcc                                    <  1:221950/149‑1 (MQ=255)
                                                                                                             aTGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTccgc                                  <  6:70885/149‑1 (MQ=255)
                                                                                                                  gCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGa                               >  4:244257/1‑147 (MQ=255)
                                                                                                                   ccAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACa                             <  3:206104/148‑1 (MQ=255)
                                                                                                                   ccAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAg                            >  7:102059/1‑149 (MQ=255)
                                                                                                                    cAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAgg                           <  4:80004/149‑1 (MQ=255)
                                                                                                                       gACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCTGTTGAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGCACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGAGCATTTTTTCCGCCGACAGGTGa                        <  8:139037/149‑1 (MQ=255)
                                                                                                                        aCGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAg                       <  7:266674/149‑1 (MQ=255)
                                                                                                                         cgcgAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGt                      <  8:266023/149‑1 (MQ=255)
                                                                                                                         cgcgAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGt                      >  5:250398/1‑149 (MQ=255)
                                                                                                                          gcgAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGt                      >  1:125117/1‑148 (MQ=255)
                                                                                                                            gAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAg                   >  5:5710/1‑149 (MQ=255)
                                                                                                                             aaCGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGc                  <  1:194766/149‑1 (MQ=255)
                                                                                                                                ggATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTa               >  6:32488/1‑149 (MQ=255)
                                                                                                                                    atCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCg           >  7:29934/1‑149 (MQ=255)
                                                                                                                                    atCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCg           >  8:160202/1‑149 (MQ=255)
                                                                                                                                     tCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCgg          <  2:281670/149‑1 (MQ=255)
                                                                                                                                       cGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCt        >  6:89444/1‑149 (MQ=255)
                                                                                                                                            gccaCGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACTCTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACATGTGAGTCAGCGTATCCGGCtttttt   <  5:389740/149‑1 (MQ=255)
                                                                                                                                            gccaCGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCtttttt   >  6:61045/1‑149 (MQ=255)
                                                                                                                                             ccaCGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTg  >  2:59029/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GTGGTCGAGATAGCGCAGCGGCTGATGAGAAAGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTG  >  NC_000913/2573754‑2574043

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