Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,574,070 C→T 5.0% C79Y (TGC→TAC)  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,574,0700CT 5.0% 258.0 / 2.2 79C79Y (TGC→TAC) eutTcobalamin adenosyltransferase involved in ethanolamine utilization
Reads supporting (aligned to +/- strand):  ref base C (37/38);  new base T (2/2);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.51e-01

ATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTC  >  NC_000913/2573923‑2574216
                                                                                                                                                   |                                                                                                                                                  
aTTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCgcagc                                                                                                                                                    <  2:54037/148‑1 (MQ=255)
 tttGCAGCCACACTGCCAGCGCAATGGGGCTGGCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCgcagca                                                                                                                                                   <  7:199400/148‑1 (MQ=255)
  ttGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGCCAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGcaccac                                                                                                                                                  >  5:7952/1‑148 (MQ=255)
  ttGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCAc                                                                                                                                                  >  3:12935/1‑148 (MQ=255)
         cacaCTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCg                                                                                                                                           >  1:342020/1‑148 (MQ=255)
         cacaCTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTCGCGACACAGTTCGCAGCACGCCTGCg                                                                                                                                           >  1:318013/1‑148 (MQ=255)
           cacTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGa                                                                                                                                         >  5:298011/1‑148 (MQ=255)
            acTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGAt                                                                                                                                        <  4:186057/148‑1 (MQ=255)
             cTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATg                                                                                                                                       <  6:298011/148‑1 (MQ=255)
                ccAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGtt                                                                                                                                     <  2:187934/147‑1 (MQ=255)
                 cAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCa                                                                                                                                   <  8:99030/148‑1 (MQ=255)
                 cAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGAGGATCGCTTTTGGCGAACATTTTTTCATCCGACAGGTGATTCAGCGTATCCGGCTTTATTGCCACCGGCTGGCGACACAGTTAGAAGCACGCCTGCGGATGTTc                                                                                                                                    >  6:177045/1‑147 (MQ=255)
                   gcgcAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATc                                                                                                                                 >  6:67871/1‑148 (MQ=255)
                    cgcAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCg                                                                                                                                >  6:129539/1‑148 (MQ=255)
                      cAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCt                                                                                                                              <  1:248498/148‑1 (MQ=255)
                      cAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCt                                                                                                                              <  6:121777/148‑1 (MQ=255)
                        aTGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTAc                                                                                                                            >  3:209785/1‑148 (MQ=255)
                         tGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACt                                                                                                                           >  4:345058/1‑148 (MQ=255)
                          ggTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTg                                                                                                                          >  1:328142/1‑148 (MQ=255)
                          ggTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTg                                                                                                                          >  4:329907/1‑148 (MQ=255)
                          ggTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTg                                                                                                                          >  3:127076/1‑148 (MQ=255)
                                 tCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGcc                                                                                                                   <  4:12935/148‑1 (MQ=255)
                                    agaACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCAt                                                                                                                >  3:49714/1‑148 (MQ=255)
                                       aCCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGTCCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTGGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGaa                                                                                                             <  5:177045/148‑1 (MQ=255)
                                               aaaaCCCAGCCGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTTCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGCTCAGCCCATGAACGggctgc                                                                                                     >  6:312310/1‑148 (MQ=255)
                                               aaaaCCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGa                                                                                                              <  5:176643/139‑1 (MQ=255)
                                               aaaaCCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGa                                                                                                              >  6:176643/1‑139 (MQ=255)
                                               aaaaCCCAGACGCGGATCGCTTTTGGCGAACATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGggctgc                                                                                                     >  5:379843/1‑148 (MQ=255)
                                                aaaCCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGggctgcg                                                                                                    >  4:262288/1‑148 (MQ=255)
                                                 aaCCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGggctgcgg                                                                                                   <  5:67871/148‑1 (MQ=255)
                                                 aaCCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGggctgcgg                                                                                                   <  4:127076/148‑1 (MQ=255)
                                                     cAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGggctgcggctgc                                                                                               >  2:283289/1‑148 (MQ=255)
                                                           cgGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGt                                                                                         >  6:268728/1‑148 (MQ=255)
                                                                   ttttGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGt                                                                                  >  4:265828/1‑147 (MQ=255)
                                                                    tttGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTc                                                                                 >  1:266252/1‑147 (MQ=255)
                                                                    tttGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCAGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTc                                                                                 >  8:339955/1‑147 (MQ=255)
                                                                          gACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCaacaaac                                                                          <  5:246252/148‑1 (MQ=255)
                                                                          gACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCaacaaac                                                                          <  1:300296/148‑1 (MQ=255)
                                                                          gACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCaacaaa                                                                           >  6:155389/1‑147 (MQ=255)
                                                                            ccATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTc                                                                                       >  7:13429/1‑133 (MQ=255)
                                                                            ccATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTc                                                                                       <  8:13429/133‑1 (MQ=255)
                                                                             cATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAgg                                                                       <  4:209785/148‑1 (MQ=255)
                                                                                 ttttCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATg                                                                                                                                       >  8:54308/1‑80 (MQ=255)
                                                                                 ttttCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATg                                                                                                                                       <  7:54308/80‑1 (MQ=255)
                                                                                     ccgccgACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCt                                                                <  1:283289/147‑1 (MQ=255)
                                                                                                tGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGTGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATa                                                     >  8:22378/1‑147 (MQ=255)
                                                                                                tGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATaa                                                    <  5:129539/148‑1 (MQ=255)
                                                                                                 gAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGCCACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTAAATaa                                                    >  6:237344/1‑147 (MQ=255)
                                                                                                 gAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATaa                                                    >  4:318777/1‑147 (MQ=255)
                                                                                                   gTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAACt                                                 <  3:262288/148‑1 (MQ=255)
                                                                                                    tCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAACtt                                                <  5:268728/148‑1 (MQ=255)
                                                                                                      aGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAACtag                                               <  3:345058/147‑3 (MQ=255)
                                                                                                      aGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTg                                               <  2:266252/147‑1 (MQ=255)
                                                                                                      aGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGGCAATAAACTTGa                                              <  8:41695/148‑1 (MQ=255)
                                                                                                          tATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCaa                                                       >  8:410043/1‑135 (MQ=255)
                                                                                                          tATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCaa                                                       <  7:410043/135‑1 (MQ=255)
                                                                                                           aTCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACt                                                                                                                           >  2:74371/1‑66 (MQ=255)
                                                                                                           aTCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACt                                                                                                                           <  1:74371/66‑1 (MQ=255)
                                                                                                           aTCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATgcg                                          >  2:81038/1‑147 (MQ=255)
                                                                                                           aTCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTAATCGTCAACAAACAGGCGACCCTGCTCGTCAATAAAATTGATgcg                                          >  6:233932/1‑147 (MQ=255)
                                                                                                            tCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAAc                                                  <  3:123493/138‑1 (MQ=255)
                                                                                                            tCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAAc                                                  >  4:123493/1‑138 (MQ=255)
                                                                                                             ccGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATgcg                                          >  7:55966/1‑145 (MQ=255)
                                                                                                              cGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGGCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAg                                       <  3:329907/147‑1 (MQ=255)
                                                                                                                tcTTTTTTTCCACCGGCTTGCGACCAAGTTTGGAGCACGCCTGCGGATGTTCATCGCTACTGGTAATCCCATTAACGTGCTGAGGCTGCTTCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGTTGCGCAGATg                                    <  6:379843/147‑1 (MQ=255)
                                                                                                                    ttttGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATggcg                                 <  8:55966/147‑1 (MQ=255)
                                                                                                                     tttGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATggcgg                                <  2:328142/147‑1 (MQ=255)
                                                                                                                       tGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCtt                             >  4:9232/1‑148 (MQ=255)
                                                                                                                         ccACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTc                           <  3:318777/148‑1 (MQ=255)
                                                                                                                           aCCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGACTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTc                           >  7:37118/1‑146 (MQ=255)
                                                                                                                                   gCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTc                  <  5:155389/147‑1 (MQ=255)
                                                                                                                                   gCGACACAGTTCGAAGCACGCCTGCGGATGTTCATCGCTACTGGGCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTc                  <  4:65814/147‑1 (MQ=255)
                                                                                                                                     gACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGGCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGTCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCg               <  5:233932/148‑1 (MQ=255)
                                                                                                                                        acaGTTCGCAGCACGCCTGCGTATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGc            <  6:19390/148‑1 (MQ=255)
                                                                                                                                         caGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCggg             >  7:172552/1‑146 (MQ=255)
                                                                                                                                         caGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCa           <  7:258406/148‑1 (MQ=255)
                                                                                                                                           gTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCaga         <  8:172552/148‑1 (MQ=255)
                                                                                                                                           gTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCaga         <  3:265828/148‑1 (MQ=255)
                                                                                                                                             tCGCAGTACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGc            >  7:55385/1‑143 (MQ=255)
                                                                                                                                             tCGCAGTACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGc            >  7:61905/1‑143 (MQ=255)
                                                                                                                                             tCGCAGTACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGc            <  8:55385/143‑1 (MQ=255)
                                                                                                                                             tCGCAGTACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGc            <  8:61905/143‑1 (MQ=255)
                                                                                                                                               gcaacaCGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAggg      <  6:7952/147‑1 (MQ=255)
                                                                                                                                                  gcaCGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTc  <  5:237344/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATCCGGCTTTTTTGCCACCGGCTGGCGACACAGTTCGCAGCACGCCTGCGGATGTTCATCGCTACTGGTCAGCCCATGAACGGGCTGCGGCTGCTGCTGTTCATCGTCAACAAACAGGCGGCCCTGCTCGTCAATAAACTTGATGCGCAGATGGCGGCTTTCCAGCAACTCCCGGGCAGAGGGCGTC  >  NC_000913/2573923‑2574216

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: