Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,145,072 C→T P360S (CCG→TCG)  pflD → putative glycine radical domain‑containing pyruvate formate‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,145,0720CT94.8% 179.4 / NA 59P360S (CCG→TCG) pflDputative glycine radical domain‑containing pyruvate formate‑lyase
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base T (21/34);  total (21/38)
Fisher's exact test for biased strand distribution p-value = 5.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTCCACCAGTAGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAACCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTG  >  NC_000913/4144925‑4145218
                                                                                                                                                    |                                                                                                                                                  
cgcTCCACCAGTAGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGa                                                                                                                                                 <  2:52627/151‑1 (MQ=255)
   tCCACCAGTAGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGAAcc                                                                                                                                              >  2:552784/1‑151 (MQ=255)
    ccaccaGTAGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGAAc                                                                                                                                               <  1:765323/149‑1 (MQ=255)
    ccaccaGTAGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGAAc                                                                                                                                               <  1:499018/149‑1 (MQ=255)
    ccaccaGTAGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGGGAACGTGCTTTCGTTCCTTTGCCTTTACGCCTATCAAAGCG‑TGCAATTACCGCAATCGaa                                                                                                                                                >  2:906545/1‑148 (MQ=255)
      accaGTAGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGTACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGAACCTCg                                                                                                                                           <  2:229884/151‑1 (MQ=255)
       ccaGTAGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGAACCTCgg                                                                                                                                          >  2:76521/1‑151 (MQ=255)
        caGTAGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGAACCTCGGc                                                                                                                                         >  1:1004612/1‑151 (MQ=255)
         aGTAGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGAACCTCGGCg                                                                                                                                        <  1:879572/151‑1 (MQ=255)
          gTAGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGCCTTAGCGCGGTTAACGTGCTTTCGGTCCTTTGCCTTGACGCCTCTCAAAGCG‑TGCAATTACCGCAATCCAACCTCTGCg                                                                                                                                        >  2:770001/1‑150 (MQ=255)
          gTAGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGa                                                                                                                                                 <  2:17618/141‑1 (MQ=255)
          gTAGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGa                                                                                                                                                 >  1:17618/1‑141 (MQ=255)
              cgcgcgTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGACCCTCGGCGTgcgc                                                                                                                                   <  2:866509/151‑1 (MQ=255)
                      aTTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGc                                                                                                                           <  2:903724/151‑1 (MQ=255)
                          cGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAAAGCCTTg                                                                                                                       >  2:195504/1‑151 (MQ=255)
                            ccGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGAt                                                                                                                     <  2:843179/151‑1 (MQ=255)
                                  tCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTGCCTTGATGCCTATCAAAGCG‑TGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGacac                                                                                                               <  2:27749/150‑1 (MQ=255)
                                   cccGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACg                                                                                                              <  2:719562/151‑1 (MQ=255)
                                     cGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGcc                                                                                                            <  1:845906/151‑1 (MQ=255)
                                     cGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGcc                                                                                                            <  2:245519/151‑1 (MQ=255)
                                     cGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGcc                                                                                                            >  2:281768/1‑151 (MQ=255)
                                           gCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCt                                                                                                      <  2:781894/151‑1 (MQ=255)
                                             tataCCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCtg                                                                                                     >  2:865257/1‑150 (MQ=255)
                                                  cGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGaaa                                                                                               <  1:895408/151‑1 (MQ=255)
                                                        gCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAAccgcc                                                                                         <  1:132370/151‑1 (MQ=255)
                                                           cggcggGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGa                                                                                       <  1:743983/150‑1 (MQ=255)
                                                            ggcggGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGaa                                                                                      >  2:877581/1‑150 (MQ=255)
                                                                  ttAACCGAGAACGGCCGTAGCGTGGTGAACGTGCTTTCGTTCCTTGGCTTTGACGCCTACCAAGGCG‑TGCAATTCCCGCATTCAAACCTCGGCGTGCGCAATATCGCTTGGCACGACACGCCGTGCCTGATGAAAACCGCCGAAACCAtt                                                                                <  2:608095/150‑1 (MQ=255)
                                                                   tAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTc                                                                               >  2:443453/1‑150 (MQ=255)
                                                                    accgcagaACGGGCGTCGCGGGGGGAACGTGTTTTCGTCCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTCCCTCAACCGACCCGGGCCGTGCGCTCTAACGCCTTGATCGACACGCCGTTCTTGATGAACACCGCCGAACCCAtt                                                                                <  2:39457/143‑1 (MQ=255)
                                                                      cagagaACGGCCGTAGCGCGGTGAACGTGTTTTCGTTCCTTGGCCTTGCCGCCTACCAAGGCT‑TGAAATCACCGCATTCGAACCTCGGCGTGCGCACTAACGCTTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCtt                                                                           <  2:830937/149‑1 (MQ=255)
                                                                              ggACGTAGCGCGGTGAACGTGCTTTCGTTGCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATGCCCGCAATCGAACCTCGGCGTGCGCACTAACGCTTTGAGCGACACGCCGTTCCTGATAAAAACCGCCGAAACCATTCGCTTCGGTACCg                                                                   <  2:23722/151‑1 (MQ=255)
                                                                              ggACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTCCCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTAcc                                                                    <  1:1060539/150‑1 (MQ=255)
                                                                              ggACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTCCCGCAATCGAACCGCGGCGTGCGCACTACCGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTAcc                                                                    <  1:562255/150‑1 (MQ=255)
                                                                              ggACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCg                                                                   >  1:420176/1‑151 (MQ=255)
                                                                                  gTAGCGCGGTTAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCgg                                                                  <  1:164321/148‑1 (MQ=255)
                                                                                     gcgcGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTcc                                                            >  2:199280/1‑151 (MQ=255)
                                                                                      cgcgTTGAACTTGCTTTCGTTCCTTTCCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGACCCTCGGCGTGCGCCCTAACCCTTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTcc                                                            <  1:991652/150‑1 (MQ=255)
                                                                                        cgGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTCCCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCg                                                           <  1:723764/149‑1 (MQ=255)
                                                                                         ggTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGACCCTCGGCGTGCGCACTAACGCCTTGATCGACCCGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCa                                                         <  2:818044/150‑1 (MQ=255)
                                                                                            gAACTTGCTTTCGTTCCTTTGCTTTGACGCCTATC‑AAGCGTTGAAATTGCCGCATGCGAACCTCGCCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAATCCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAAt                                                      <  1:910053/150‑1 (MQ=255)
                                                                                             aaCGTGCTTGCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGAACCTGGGCGTGCGCCCTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATGCGCTTCGGTACCGGTATTCCGCAAAt                                                      <  2:561932/149‑1 (MQ=255)
                                                                                                 tGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTa                                                 >  1:8589/1‑150 (MQ=255)
                                                                                                    tttCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGAAATTACCGCAATCGACCCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTa                                                 <  2:1053892/147‑1 (MQ=255)
                                                                                                        gTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGa                                          >  1:52683/1‑150 (MQ=255)
                                                                                                         ttCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGa                                          >  2:697160/1‑149 (MQ=255)
                                                                                                         ttCCTTGGCTTTGACGCCTATCAAAGCG‑GGCAATTCCCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGa                                          <  1:868690/149‑1 (MQ=255)
                                                                                                          tCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAgt                                       >  1:628212/1‑151 (MQ=255)
                                                                                                           ccTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAgtg                                      <  1:562186/151‑1 (MQ=255)
                                                                                                              ttGCCTTGACGCCTATCAAGGCG‑TGCAATTACCGCAATCGACCCCGGGCGTGCGCACTAACTCTTGGATCGCCACGCCGTTCTTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAgtgg                                     <  1:836003/149‑1 (MQ=255)
                                                                                                               tGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAgtggtgg                                  >  1:290323/1‑151 (MQ=255)
                                                                                                                        tccTATAAAAGCG‑TGCAATTACCGAAACCGAACCTCGGCGTGCGCTCTAACGCTTTGACCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGc                           <  2:116753/148‑1 (MQ=255)
                                                                                                                         ccTATCAAAGCG‑TGCAATTCCCGCAATCGACCCCCGGCGTGCGCACTACCGCTTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGt                         <  2:470084/150‑1 (MQ=255)
                                                                                                                          cTATCAAAGCG‑TGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTc                       >  1:477872/1‑151 (MQ=255)
                                                                                                                                    g‑TGCAATTCCCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGt              <  1:717415/150‑1 (MQ=255)
                                                                                                                                         aaTTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGt         >  1:937028/1‑151 (MQ=255)
                                                                                                                                            tcccGCAATCGACCCTCGGCGTGCGCCCTAACGCCTGGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGttt       <  1:462328/148‑1 (MQ=255)
                                                                                                                                              ccGCAACCGACCCTCGGCGTGCGCCCTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGt         <  2:995382/146‑1 (MQ=255)
                                                                                                                                                gCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTg  >  2:230236/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CGCTCCACCAGTAGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCG‑TGCAATTACCGCAACCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGTTTCGCTG  >  NC_000913/4144925‑4145218

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑