Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 944,703 G→A G224D (GGT→GAT)  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,7030GA98.4% 229.7 / NA 63G224D (GGT→GAT) dmsCdimethyl sulfoxide reductase, anaerobic, subunit C
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (29/33);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGGTGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCTGCGGTTCCGCTACTTTCTGTATCATTCATTCTGCTTCTGGCAGGGGAATTAATCGGTCGCGGCG  >  NC_000913/944558‑944849
                                                                                                                                                 |                                                                                                                                                  
atCGATTGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGTCACTGGTAGTGCGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTTGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATTATgcg                                                                                                                                                >  1:149895/2‑150 (MQ=255)
  cGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCt                                                                                                                                              <  2:585535/150‑1 (MQ=255)
  cGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCt                                                                                                                                              <  1:570610/150‑1 (MQ=255)
  cGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTg                                                                                                                                             >  1:633294/1‑151 (MQ=255)
    aTGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGAt                                                                                                                                           <  2:184414/151‑1 (MQ=255)
          gggCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCt                                                                                                                                      <  1:690415/150‑1 (MQ=255)
             cGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCg                                                                                                                                  >  1:814496/1‑151 (MQ=255)
                   gTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGt                                                                                                                            >  2:53464/1‑151 (MQ=255)
                    tCTGCTTCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGGGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTg                                                                                                                           >  1:59305/1‑151 (MQ=255)
                            cAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGc                                                                                                                   >  1:405925/1‑151 (MQ=255)
                              gCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCg                                                                                                                 >  1:640224/1‑151 (MQ=255)
                              gCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCg                                                                                                                 >  2:5075/1‑151 (MQ=255)
                                aaTTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGtt                                                                                                               <  1:354017/151‑1 (MQ=255)
                                  tttCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTg                                                                                                              >  1:225886/1‑150 (MQ=255)
                                    tCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGccc                                                                                                           >  1:469159/1‑151 (MQ=255)
                                       gTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCtgt                                                                                                        <  2:117917/151‑1 (MQ=255)
                                       gTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTAAGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCt                                                                                                          <  2:468095/149‑1 (MQ=255)
                                           tGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCtt                                                                                                    >  1:508831/1‑151 (MQ=255)
                                                   gttgtgagtgGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCTACTCTTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGGCTATGGTGCGCTGATGGTCTGGCGGATCTTGCTTTTGGCCGTTGCCCTGTGCTTGTGGAtt                                                                                             <  1:602952/147‑1 (MQ=255)
                                                    tAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGc                                                                                           >  2:598026/1‑151 (MQ=255)
                                                        gagtgGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGcacca                                                                                       <  1:79073/151‑1 (MQ=255)
                                                            ggCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAg                                                                                    <  1:18092/150‑1 (MQ=255)
                                                            ggCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGc                                                                                   <  2:586244/151‑1 (MQ=255)
                                                            ggCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGc                                                                                   >  1:767775/1‑151 (MQ=255)
                                                               gtggtgTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTaa                                                                                <  1:125421/151‑1 (MQ=255)
                                                                 ggtgTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTaaa                                                                               >  2:661456/1‑150 (MQ=255)
                                                                  gtgtCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTaaa                                                                               >  2:216054/1‑149 (MQ=255)
                                                                      cAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGtt                                                                          <  2:802022/150‑1 (MQ=255)
                                                                        gTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATc                                                                       <  2:315835/151‑1 (MQ=255)
                                                                            tGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGcc                                                                   >  1:778320/1‑151 (MQ=255)
                                                                             gCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGcc                                                                   <  1:53442/150‑1 (MQ=255)
                                                                                 ggCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCTGCgg                                                              >  1:65888/1‑151 (MQ=255)
                                                                                       gagCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCTGCGGTTCCGc                                                        <  2:158534/151‑1 (MQ=255)
                                                                                         gCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCTGCGGTTCCGCTa                                                      <  1:234938/151‑1 (MQ=255)
                                                                                               gACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCTGCGGTTCCGCTACTTTCt                                                <  2:237247/151‑1 (MQ=255)
                                                                                                  tATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCTGCGGTTCCGCTACTTTCTGt                                              >  1:842329/1‑150 (MQ=255)
                                                                                                         aGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCTGCGGTTCCGCTACTTTCTGTAtcattca                                      <  1:210719/151‑1 (MQ=255)
                                                                                                         aGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCTGCGGTTCCGCTACTTTCTGTAtcattc                                       >  1:759258/1‑150 (MQ=255)
                                                                                                         aGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCTGCGGTTCCGCTACTTTCTGTAtcattc                                       <  2:56002/150‑1 (MQ=255)
                                                                                                          gTTCTTTGCAGCAGGACGCAGAGCTGGTGACGGACTATGATGCGCTGGTGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGTCTGCGGTTCCGCTACTTTCTGTAtcattcat                                     >  1:791950/1‑151 (MQ=255)
                                                                                                          gTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCTGCGGTTCCGCTACTTTCTGTAtcattcat                                     <  2:640266/151‑1 (MQ=255)
                                                                                                               gtgCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCTGCGGTTCCGCTACTTTCTGTATCATTCATTCt                                  <  2:800683/149‑1 (MQ=255)
                                                                                                               gtgCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCTGCGGTTCCGCTACTTTCTGTATCATTCATTCTg                                 <  2:225910/150‑1 (MQ=255)
                                                                                                                tgCAGCAGGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAATGTTATCAGCCTGCGGTTCCGCTACTTTCTGTATCATTCATTCTg                                 <  1:696840/149‑1 (MQ=255)
                                                                                                                  cagcagGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCTGCGGTTCCGCTACTTTCTGTATCATTCATTCTGCTTc                             <  1:858111/151‑1 (MQ=255)
                                                                                                                  cagcagGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCTGCGGTTCCGCTACTTTCTGTATCATTCATTCTGCTTc                             <  2:763450/151‑1 (MQ=255)
                                                                                                                    gcagGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCTGCGGTTCCGCTACTTTCTGTATCATTCATTCTGCTTCTg                           <  2:363896/151‑1 (MQ=255)
                                                                                                                    gcagGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGGTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCTGCGGTTCCGCTACTTTCTGTATCATTCATTCTGCTTCt                            >  1:251288/1‑150 (MQ=255)
                                                                                                                     cagGCCGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCTGCGGTTCCGCTACTTTCTGTATCATTCATTCTGCTTCTgg                          <  2:778584/151‑1 (MQ=255)
                                                                                                                          cGCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCTGCGGTTCCGCTACTTTCTGTATCATTCATTCTGCTTCTGGCAggg                     <  2:359274/151‑1 (MQ=255)
                                                                                                                           gCAGCGCTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCTGCGGTTCCGCTACTTTCTGTATCATTCATTCTGCTTCTGGCAgggg                    >  1:58180/1‑151 (MQ=255)
                                                                                                                                tcTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCTGCGGTTCCGCTACTTTCTGTATCATTCATTCTGCTTCTGGCAGGGGaa                  >  1:709524/2‑148 (MQ=255)
                                                                                                                                gcTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCTGCGGTTCCGCTACTTTCTGTATCATTCATTCTGCTTCTGGCAGGGGAAtt                >  2:58865/1‑150 (MQ=255)
                                                                                                                                gcTGGTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCTGCGGTTCCGCTACTTTCTGTATCATTCATTCTGCTTCTGGCAGGGGAATTa               <  1:30827/151‑1 (MQ=255)
                                                                                                                                   ggTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCTGCGGTTCCGCTACTTTCTGTATCATTCATTCTGCTTCTGGCAGGGGAATTaa              <  2:684895/149‑1 (MQ=255)
                                                                                                                                   ggTGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCTGCGGTTCCGCTACTTTCTGTATCATTCATTCTGCTTCTGGCAGGGGAATTAAt             <  2:41385/150‑1 (MQ=255)
                                                                                                                                     tGCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCTGCGGTTCCGCTACTTTCTGTATCATTCATTCTGCTTCTGGCAGGGGAATTAATCgg          >  2:509018/1‑151 (MQ=255)
                                                                                                                                      gCCGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCTGCGGTTCCGCTACTTTCTGTATCATTCATTCTGCTTCTGGCAGGGGAATTAATCgg          <  1:210439/150‑1 (MQ=255)
                                                                                                                                       ccGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCTGCGGTTCCGCTACTTTCTGTATCATTCATTCTGCTTCTGGCAGGGGAATTAATCGGTc        >  1:395383/1‑151 (MQ=255)
                                                                                                                                       ccGGACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCTGCGGTTCCGCTACTTTCTGTATCATTCATTCTGCTTCTGGCAGGGGAATTAATCGGTc        >  1:338441/1‑151 (MQ=255)
                                                                                                                                          gACTATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCTGCGGTTCCGCTACTTTCTGTATCATTCATTCTGCTTCTGGCAGGGGAATTAATCGGTcgcg     <  2:469376/151‑1 (MQ=255)
                                                                                                                                             tATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCTGCGGTTCCGCTACTTTCTGTATCATTCATTCTGCTTCTGGCAGGGGAATTAATCGGTCgcggcg  >  1:219360/1‑151 (MQ=255)
                                                                                                                                             tATGATGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCTGCGGTTCCGCTACTTTCTGTATCATTCATTCTGCTTCTGGCAGGGGAATTAATCGGTCgcggcg  >  1:712750/1‑151 (MQ=255)
                                                                                                                                              atgatgCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCTGCGGTTCCGCTACTTTCTGTATCATTCATTCTGCTTCTGGCAGGGGAATTAATCGGTCgcggcg  <  2:759315/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GTCGATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCTGGTGCCGGACTATGGTGCGCTGATGTCCTGGCGGATCGTGCTTTTGGCCGTTGCCCTGTGCTTGTGGATTGCACCACAGCTAAAAGGTTATCAGCCTGCGGTTCCGCTACTTTCTGTATCATTCATTCTGCTTCTGGCAGGGGAATTAATCGGTCGCGGCG  >  NC_000913/944558‑944849

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