Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 944,403 T→C 5.1% V124A (GTT→GCT)  dmsC → dimethyl sulfoxide reductase, anaerobic, subunit C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913944,4030TC 5.1% 237.9 / 2.9 79V124A (GTT→GCT) dmsCdimethyl sulfoxide reductase, anaerobic, subunit C
Reads supporting (aligned to +/- strand):  ref base T (33/40);  new base C (2/2);  total (36/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.82e-01

GCTCAACCGGGTAGGGGCTTCAGCACTCAGTAACGAAATCGCCAGCGGTTCGATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAGTCTG  >  NC_000913/944257‑944551
                                                                                                                                                  |                                                                                                                                                    
gCTCAACCGGGTAGGGGCTTCAGCACTCAGTAACGAAATCGCCAGCGGTTCGATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAACGCTGTCACCTGCATTGCGTACGCTTTGGCTGATAGTGACGATGGTTc                                                                                                                                                    <  2:159933/149‑1 (MQ=255)
 cTCAACCGGGTAGGGGCTTCAGCACTCAGTAACGAAATCGCCAGCGGTTCGATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCt                                                                                                                                                   >  4:396898/1‑149 (MQ=255)
  tCAACCGGGTAGGGGCTTCAGCACTCAGTAACGAAATCGCCAGCGGTTCGATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCtt                                                                                                                                                  >  2:409279/1‑149 (MQ=255)
   cAACCGGGTAGGGGCTTCAGCACTCAGTAACGAAATCGCCAGCGGTTCGATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTg                                                                                                                                                 <  3:257444/149‑1 (MQ=255)
    aaCCGGGTAGGGGCTTCAGCACTCAGTAACGAAATCGCCAGCGGTTCGATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTgg                                                                                                                                                >  8:233340/1‑149 (MQ=255)
        gggTAGGGGCTTCAGCACTCAGTAACGAAATCGCCAGCGGTTCGATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTc                                                                                                                                            <  7:151134/149‑1 (MQ=255)
          gTAGGGGCTTCAGCACTCAGTAACGAAATCGCCAGCGGTTCGATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGCTCTTGGCGTCAt                                                                                                                                          >  4:276280/1‑149 (MQ=255)
                 cTTCAGCACTCAGTAACGAAATCGCCAGCGGTTCGATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTc                                                                                                                                   <  3:403978/149‑1 (MQ=255)
                 cTTCAGCACTCAGTAACGAAATCGCCAGCGGTTCGATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGGACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTc                                                                                                                                   <  8:305718/149‑1 (MQ=255)
                    cAGCACTCAGTAACGAAATCGCCAGCGGTTCGATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCt                                                                                                                                  <  2:20724/147‑1 (MQ=255)
                         cTCAGTAACGAAATCGCCAGCGGTTCGATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGgatgat                                                                                                                           >  1:421388/1‑149 (MQ=255)
                            aGTAACGAAATCGCCAGCGGTTCGATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGatgatg                                                                                                                          <  5:150591/147‑1 (MQ=255)
                            aGTAACGAAATCGCCAGCGGTTCGATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGatgatg                                                                                                                          <  1:332268/147‑1 (MQ=255)
                            aGTAACGAAATCGCCAGCGGTTCGATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGCTCTTGGCGTCATCTTTGTCTGGATGATGGt                                                                                                                        >  2:342266/1‑149 (MQ=255)
                              tAACGAAATCGCCAGCGGTTCGATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCGGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTGGTGGGCGTCAGCTTTGTCTGGATGATGGTGc                                                                                                                      >  4:288939/1‑149 (MQ=255)
                                aCGAAATCGCCAGCGGTTCGATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCgt                                                                                                                    <  3:264529/149‑1 (MQ=255)
                                aCGAAATCGCCAGCGGTTCGATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCgt                                                                                                                    <  5:254123/149‑1 (MQ=255)
                                 cGAAATCGCCAGCGGTTCGATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCgtg                                                                                                                   >  4:75765/1‑149 (MQ=255)
                                  gAAATCGCCAGCGGTTCGATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTg                                                                                                                                                 >  4:247105/1‑118 (MQ=255)
                                  gAAATCGCCAGCGGTTCGATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTg                                                                                                                                                 <  3:247105/118‑1 (MQ=255)
                                  gAAATCGCCAGCGGTTCGATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCgtgt                                                                                                                  >  3:395213/1‑149 (MQ=255)
                                  gAAATCGCCAGCGGTTCGATCTTTTTTGCCGTACGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCAATGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTATGTCTGGATGATGGTGCgtgt                                                                                                                  >  2:334211/1‑149 (MQ=255)
                                        gCCAGCGGTTCGATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATaa                                                                                                            <  5:215729/149‑1 (MQ=255)
                                        gCCAGCGGTTCGATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCAACGGCATTGCGTACGCAGTGGCAGATAGTGACGATGGTTCTTGGAGTCATCTATCTCTGGATGATGGAGCGTGTGTATaa                                                                                                            >  3:292570/1‑149 (MQ=255)
                                          cAGCGGTTCGATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCgtgt                                                                                                                  >  4:38665/1‑141 (MQ=255)
                                          cAGCGGTTCGATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCgtgt                                                                                                                  <  3:38665/141‑1 (MQ=255)
                                              ggTTCGATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCAt                                                                                                      <  3:112016/149‑1 (MQ=255)
                                               gTTCGATCTTTTTTTCCGTTGGCGGCATGTGCTGTTTTCTTGCAATGCTGAATAAGCTGTCTCCGGCATTGCGTACGATGTGGATGATATTGACGATGGTTCTTGGGGTCATCTTTGTCTGGATGATGGTGCGTGTGTTTAACAGCAtt                                                                                                     <  4:292570/149‑1 (MQ=255)
                                                ttCGATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTg                                                                                                    <  2:421388/149‑1 (MQ=255)
                                                ttCGATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTg                                                                                                    <  3:75765/149‑1 (MQ=255)
                                                 tggagcgtttttGCCGTAGGGGGATTCGGGGGGTTGCTTGCAATGCTGAAACAGCTGTCACCGGCATTGCGGAGGGTGTGGCTTATCGTGACGATGGTTCTAGGCGTCATCTTAGTCTGGATGATGGGGCGTGTGTATAACAGCATTGa                                                                                                   <  6:136734/142‑1 (MQ=255)
                                                  cGATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGCTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGa                                                                                                   <  1:342266/148‑1 (MQ=255)
                                                   gATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGAt                                                                                                  <  6:146172/148‑1 (MQ=255)
                                                    aTCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACa                                                                                                          >  5:294409/1‑139 (MQ=255)
                                                    aTCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGGCTGGATGATGGTGCGTGTGTATAACa                                                                                                          <  6:294409/139‑1 (MQ=255)
                                                     tCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCAt                                                                                                                                          >  7:180280/1‑106 (MQ=255)
                                                     tCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCAt                                                                                                                                          <  8:180280/106‑1 (MQ=255)
                                                     tCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATAcc                                                                                               >  7:273497/1‑149 (MQ=255)
                                                     tCTTTTTTGACGTAGGCGGCATCGGCTGGTTGCTGGCAATGATGAAAAATCAGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGGTCATGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTTTAACAGCATTGATAcc                                                                                               >  5:136056/1‑149 (MQ=255)
                                                         ttttGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATAc                                                                                                <  6:88121/144‑1 (MQ=255)
                                                         ttttGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGGGCGTGTGTATAACAGCATTGATACCGtt                                                                                            <  5:425052/148‑1 (MQ=255)
                                                                    gcggcATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTAc                                                                                <  8:273497/149‑1 (MQ=255)
                                                                     cggcATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACa                                                                               <  7:233340/149‑1 (MQ=255)
                                                                     cggcATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACa                                                                               <  3:176871/149‑1 (MQ=255)
                                                                      ggcATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAg                                                                              <  3:288939/149‑1 (MQ=255)
                                                                         aTCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCAt                                                                           >  3:97103/1‑149 (MQ=255)
                                                                          tCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATc                                                                          >  1:6195/1‑149 (MQ=255)
                                                                          tCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATc                                                                          >  2:432481/1‑149 (MQ=255)
                                                                            ggCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTTGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTg                                                                        <  8:180956/149‑1 (MQ=255)
                                                                              cTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTg                                                                                                                       <  6:380259/100‑1 (MQ=255)
                                                                              cTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTg                                                                                                                       >  5:380259/1‑100 (MQ=255)
                                                                                 gTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATAGCGTTCCGACCTGGTACAGCATCTGGACGc                                                                   <  1:147774/149‑1 (MQ=255)
                                                                                 gTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGCTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGc                                                                   <  3:276280/149‑1 (MQ=255)
                                                                                        gCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCATCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATggg                                                            >  3:41480/1‑149 (MQ=255)
                                                                                         cAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGc                                                           >  6:416261/1‑149 (MQ=255)
                                                                                            tGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGcttc                                                        >  7:169794/1‑149 (MQ=255)
                                                                                               tGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACAGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGcttctt                                                      >  2:429911/1‑148 (MQ=255)
                                                                                                 aaaaaGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGGGCTTGGCGTCATCGTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCttcttc                                                     <  1:334211/147‑1 (MQ=255)
                                                                                                        tGTCACCGGAATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAAAAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCCATTGGCTTCTTCATGACGATg                                            >  4:250329/1‑149 (MQ=255)
                                                                                                            aCCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGct                                                          <  4:21013/131‑1 (MQ=255)
                                                                                                            aCCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGct                                                          >  3:21013/1‑131 (MQ=255)
                                                                                                            aCCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTa                                        >  3:365001/1‑149 (MQ=255)
                                                                                                               ggCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATgg                                     <  2:6195/149‑1 (MQ=255)
                                                                                                                        tACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGt                            <  5:416261/149‑1 (MQ=255)
                                                                                                                          cGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGGGTGTATAACAGCTTTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGtt                           <  3:396898/148‑1 (MQ=255)
                                                                                                                               tgGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTggg                     <  7:119185/149‑1 (MQ=255)
                                                                                                                                 gCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGt                    >  5:148845/1‑148 (MQ=255)
                                                                                                                                     aTAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCt               >  7:279186/1‑149 (MQ=255)
                                                                                                                                        gTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCtgtt            >  4:159641/1‑149 (MQ=255)
                                                                                                                                         tGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGtt                           >  7:230944/1‑133 (MQ=255)
                                                                                                                                         tGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGtt                           <  8:230944/133‑1 (MQ=255)
                                                                                                                                         tGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCtgttg           <  1:409279/149‑1 (MQ=255)
                                                                                                                                           aCGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCtgttgt          <  4:395213/148‑1 (MQ=255)
                                                                                                                                             gATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGa       <  4:97103/149‑1 (MQ=255)
                                                                                                                                              aTGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATAc                                                                                                <  4:398709/59‑1 (MQ=255)
                                                                                                                                              aTGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATAc                                                                                                >  3:398709/1‑59 (MQ=255)
                                                                                                                                              aTGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAg      >  3:250800/1‑149 (MQ=255)
                                                                                                                                              aTGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGATTCTTCCTGACGATGTTTATGGGCGGCACGTTACTGGGTTACCTGTTGTTGAg      >  5:76459/1‑149 (MQ=255)
                                                                                                                                                ggTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAGTc    >  7:45607/1‑149 (MQ=255)
                                                                                                                                                  ttCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAGTctg  >  4:270595/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCTCAACCGGGTAGGGGCTTCAGCACTCAGTAACGAAATCGCCAGCGGTTCGATCTTTTTTGCCGTAGGCGGCATCGGCTGGTTGCTGGCAATGCTGAAAAAGCTGTCACCGGCATTGCGTACGCTGTGGCTGATAGTGACGATGGTTCTTGGCGTCATCTTTGTCTGGATGATGGTGCGTGTGTATAACAGCATTGATACCGTTCCGACCTGGTACAGCATCTGGACGCCGATGGGCTTCTTCCTGACGATGTTTATGGGCGGCCCGTTACTGGGTTACCTGTTGTTGAGTCTG  >  NC_000913/944257‑944551

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