Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,422,131 G→A 5.6% A142T (GCT→ACT)  yfcH → putative NAD‑dependent nucleotide‑sugar epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,422,1310GA 5.6% 228.0 / 2.0 71A142T (GCT→ACT) yfcHputative NAD‑dependent nucleotide‑sugar epimerase
Reads supporting (aligned to +/- strand):  ref base G (30/36);  new base A (2/2);  total (32/39)
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.66e-01

AAAAACTGGTCGATTTGATTAATGCCAGCGACACGCCACCGTCGGTACTCATTTCCGGCTCGGCAACGGGCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATT  >  NC_000913/2421985‑2422274
                                                                                                                                                  |                                                                                                                                               
aaaaaCTGGTCGATTTGATTAATGCCAGCGCCACGCCACCGTCGGTACTCATTTCCGGCTCGGCAACGGGCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGcgct                                                                                                                                               >  5:310080/1‑149 (MQ=255)
aaaaaCTGGTCGATTTGATTAATGCCAGCGACACGCCACCGTCGGTACTCATTTCCGGCTCGGCAACGGGCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGcgct                                                                                                                                               <  1:339782/149‑1 (MQ=255)
     cTGGTCGATTTGATTAATGCCAGCGACACGCCACCGTCGGTACTCATTTCCGGCTCGGCAACGGGCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTg                                                                                                                                          <  4:88177/149‑1 (MQ=255)
     cTGGGCGATTTGATTAATGCCAGCGACACGCCACCGTCGGTACTCATTTCCGGCTCGGCAACGGGCTATTATGGTGACTTAGGTGAAGTGGTGGGTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTg                                                                                                                                          <  7:346634/149‑1 (MQ=255)
         tCGATTTGATTAATGCCAGCGACACGCCACCGTCGGTACTCATTTCCGGCTCGGCAACGGGCTATTATGGTGACTTAGGTGAAGTTGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGgaa                                                                                                                                      <  7:240598/149‑1 (MQ=255)
           gATTTGATTAATGCCAGCGACACGCCACCGTCGGTACTCATTTCCGGCTCGGCAACGGGCTATTATGGTGTCTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGgaaga                                                                                                                                    <  4:9009/149‑1 (MQ=255)
            aTTTGATTAATGCCAGCGACACGCCACCGTCGGTACTCATTTCCGGCTCGGCAACGGGCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGgaaga                                                                                                                                    >  8:327873/1‑148 (MQ=255)
               tGATTAATGCCAGCGACACGCCACCGTCGGTACTCATTTCCGGCTCGGCAACGGGCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAAt                                                                                                                                 <  5:331120/148‑1 (MQ=255)
                 aTTAATGCCAGCGACACGCCACCGTCGGTACTCATTTCCGGCTCGGCAACGGGCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATtgc                                                                                                                              >  3:64265/1‑149 (MQ=255)
                         cAGCGACACGCCACCGTCGGTACTCATTTCCGGCTCGGCAACGGGCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCggg                                                                                                                      >  7:42437/1‑149 (MQ=255)
                          aGCGACACGCCACCGTCGGTACTCATTTCCGGCTCGGCAACGGGCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGgc                                                                                                                     <  3:329433/149‑1 (MQ=255)
                                acGCCACCGTCGGTACTCATTTCCGGCTCGGCAACGGGCTATTATGGTGACATAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTTCCTGCCGGGCGCAAAg                                                                                                               >  1:261922/1‑149 (MQ=255)
                                        gTCGGTACTCATTTCCGGCTCGGCAACGGGCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACaaaa                                                                                                       <  2:17282/149‑1 (MQ=255)
                                           ggTACTCATTTCCCGATCGGCAACGGGCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAAcgc                                                                                                    <  2:72579/149‑1 (MQ=255)
                                             tACTCATTTCCGGCTCGGCAACGGGCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAAGCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGa                                                                                                  <  1:33043/149‑1 (MQ=255)
                                             tACTCATTTCCGGCTCGGCAACGGGCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGa                                                                                                  <  7:340889/149‑1 (MQ=255)
                                             tACTCATTTCCGGCTCGGCAACGGGCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGa                                                                                                  <  7:21771/149‑1 (MQ=255)
                                             tACTCATTTCCGGCTCGGCAACGGGCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGa                                                                                                  <  7:16875/149‑1 (MQ=255)
                                             tACTCAATTCCGGCGCGGCAACCGGCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGa                                                                                                  <  8:320876/149‑1 (MQ=255)
                                               cTCATTTGGGGCTCGGCTTCGGGCTATTTTGGTGACTTAGGTGACGTGGTGGTTACCGAAGAGGAACCGGCGCATAACGATTTTAGCCATAATCTCTGCGCTGGCTGGGAAGTAATTGCCTGCGGGGGGCAAAGTGTCAAAACGCGAg                                                                                                 <  1:131752/148‑1 (MQ=255)
                                               cTCATTTCCGGCTCGGCAACGGGCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAg                                                                                                 <  8:129204/148‑1 (MQ=255)
                                                  aTTTCCGGCTCGGCAACGGGCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAgtgt                                                                                              <  7:327873/148‑1 (MQ=255)
                                                  aTTTCCGGCTCGGCAACGGGCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGtgtg                                                                                             >  7:307771/1‑149 (MQ=255)
                                                   tttCCGGCTCGGCAACGGGCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGc                                                                                            <  4:104907/149‑1 (MQ=255)
                                                       cGGCTCGGCAACGGGCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCctgc                                                                                        <  4:64265/149‑1 (MQ=255)
                                                          cgcGGCAACGTGCTATTATGTTGACTTAGGTTAAGTGGTCGTTACCGAAGAGGAACCGCAGCATAACGAATTTACCCATAAACTCTTCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGATTGTGCCtgctgc                                                                                     <  2:276399/147‑1 (MQ=255)
                                                           tCGGCAACGGGCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCg                                                                                    >  8:66611/1‑149 (MQ=255)
                                                              gCAACGGGCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTAc                                                                                 <  2:50541/149‑1 (MQ=255)
                                                                  cGGGCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGgt                                                                             >  8:132822/1‑149 (MQ=255)
                                                                   gggCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGgtg                                                                            >  7:318265/1‑149 (MQ=255)
                                                                   gggCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCAGCGTACCGgtg                                                                            >  8:226820/1‑149 (MQ=255)
                                                                     gCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGgtg                                                                            >  2:77708/1‑147 (MQ=255)
                                                                         ttatGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGc                                                                                            >  6:217465/1‑127 (MQ=255)
                                                                         ttatGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGc                                                                                            <  5:217465/127‑1 (MQ=255)
                                                                         ttatGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGc                                                                      <  7:66611/149‑1 (MQ=255)
                                                                           atGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCt                                                                     <  8:307771/148‑1 (MQ=255)
                                                                              gTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCTAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCa                                                                 >  4:53684/1‑149 (MQ=255)
                                                                                   ttAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGa                                                            <  2:299647/149‑1 (MQ=255)
                                                                                        tGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTgg                                                                   <  2:337102/137‑1 (MQ=255)
                                                                                        tGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTgg                                                                   >  1:337102/1‑137 (MQ=255)
                                                                                            gtggtgGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTAt                                                   >  4:256483/1‑149 (MQ=255)
                                                                                             tggtgtTTACCGAAGAGGAAACTCCGCATAACGAATTTACCCATAAAATCTGCACTCTCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGATTGTGCCTGCTGCGTACCg                                                                               <  6:202283/120‑1 (MQ=255)
                                                                                             tggtggTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCACTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCg                                                                               >  5:197267/1‑120 (MQ=255)
                                                                                             tggtggTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCACTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCg                                                                               <  6:197267/120‑1 (MQ=255)
                                                                                             tggtggTTACCGAAGAGGAACCGCCGAATAACGAATTTACACATAAACTATGCACTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCg                                                                               >  5:202283/1‑120 (MQ=255)
                                                                                               gtggtTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCATGCTGCGTCCCGGTGTAGTGCTGGCACAGGAAGGCGGTATtct                                                >  1:28511/1‑149 (MQ=255)
                                                                                                tggtTACCGAAGAGGCACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATtctc                                               >  1:224446/1‑149 (MQ=255)
                                                                                                tggtTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATtctc                                               <  6:310080/149‑1 (MQ=255)
                                                                                                ctatctccGAAGAGGAAGCGCCGCATAACGAATTTACCCATAAACTCTGCTCTCGCTGGGAAGAAATTGCCTTCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATtctc                                               <  7:163622/143‑1 (MQ=255)
                                                                                                  tttACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCGTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCg                                              <  1:121602/147‑1 (MQ=255)
                                                                                                  ggtACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGAGACAAAACGCGAGTGTGCCTGCTGCGTACCCGTGTAGTGCTGGAACCGGATGGCGGTATTCTCgg                                             >  7:74385/3‑149 (MQ=255)
                                                                                                    tACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCAGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGt                                            >  5:212827/1‑148 (MQ=255)
                                                                                                     aCCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTa                                           <  3:162576/148‑1 (MQ=255)
                                                                                                       cGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATtctc                                               >  8:163622/1‑142 (MQ=255)
                                                                                                       cGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTaaaa                                        >  7:259139/1‑149 (MQ=255)
                                                                                                       cGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTaaaa                                        >  3:61704/1‑149 (MQ=255)
                                                                                                            agGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAAtgctg                                   <  1:48454/149‑1 (MQ=255)
                                                                                                              gAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTgcc                                 >  5:333889/1‑149 (MQ=255)
                                                                                                                  cgccgcATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGccgccg                             >  6:148422/1‑149 (MQ=255)
                                                                                                                      gcATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGt                                                     <  8:211886/121‑1 (MQ=255)
                                                                                                                      gcATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGt                                                     >  7:211886/1‑121 (MQ=255)
                                                                                                                         tAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTc                      >  8:174587/1‑149 (MQ=255)
                                                                                                                          aaCGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCt                     >  4:185933/1‑149 (MQ=255)
                                                                                                                            cGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTg                   >  1:28913/1‑149 (MQ=255)
                                                                                                                             gAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCGGCTGCGTACCGGGGTAGTGCTGGCACCGGATGGCGGTCTTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTgg                  >  4:285834/1‑149 (MQ=255)
                                                                                                                              aaTTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGc                 >  1:55914/1‑149 (MQ=255)
                                                                                                                                tttACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGt                            >  1:317715/1‑136 (MQ=255)
                                                                                                                                tttACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGt                            <  2:317715/136‑1 (MQ=255)
                                                                                                                                    cccATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGGGCCTGCTGCGTACCGGGGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTgggc           <  2:141061/149‑1 (MQ=255)
                                                                                                                                     ccATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTgggcg          <  5:148422/149‑1 (MQ=255)
                                                                                                                                        tAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTgggcgggc       <  7:33789/149‑1 (MQ=255)
                                                                                                                                         aaaCTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGGTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTgggcgggc       <  8:287720/148‑1 (MQ=255)
                                                                                                                                               tGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGGGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGAtt  <  4:41056/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
AAAAACTGGTCGATTTGATTAATGCCAGCGACACGCCACCGTCGGTACTCATTTCCGGCTCGGCAACGGGCTATTATGGTGACTTAGGTGAAGTGGTGGTTACCGAAGAGGAACCGCCGCATAACGAATTTACCCATAAACTCTGCGCTCGCTGGGAAGAAATTGCCTGCCGGGCGCAAAGTGACAAAACGCGAGTGTGCCTGCTGCGTACCGGTGTAGTGCTGGCACCGGATGGCGGTATTCTCGGTAAAATGCTGCCGCCGTTTCGTCTTGGCCTGGGCGGGCCGATT  >  NC_000913/2421985‑2422274

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