Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 995,030 G→A A320V (GCC→GTC)  ssuD ← FMNH2‑dependent alkanesulfonate monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913995,0300GA97.4% 110.3 / NA 38A320V (GCC→GTC) ssuDFMNH2‑dependent alkanesulfonate monooxygenase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (25/12);  minor base C (0/1);  total (25/13)
Fisher's exact test for biased strand distribution p-value = 3.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

AAACCAGGGGGCAACGCGCAATAACCACTTCTTCACTGGCGTTGCCATTATTTCTTCCTTAGCTTTGCGCGACTTTACGGGGGATAAAATCATTCGCCACCGCTTCGCCTTGCGGATTCAGCGGCTGCGGCTGGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAACTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGCCAAGCGCGGCATATTCGTTGATTCGCGCAGCGACCGTAGGACCATCGCCCACCAGCGCCGTCCCGGCACCGCCGCGCACTAAGCCAACGCCCGCCCATAAATTGGGGCTGATCTCCAGATTGTCGCGCTTGCCGTTATGTAACGCCGCCATTCGCTGTTGCCCTACGGAATCCGTCCGGGCGAATGCGGCCTGTGCTTTGGCGATAGTTTCATCATCAAGATGCGAGATTAACCGCTCGGCGGCCTGCC  >  NC_000913/994785‑995270
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                
aaaCCAGGGGGCAACGCGCAATAACCACTTCTTCACTGGCGTTGCCATTATTTCTTCCTTAGCTTTGCGCGACTTTACGGGGGATAAAATCATTCGCCACCGCTTCGCCTTGCGGATTCAGCTGCTGCGGCTGGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAACTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGCCAAGCGCGACatat                                                                                                                                                                                                                                             >  1:303052/1‑251 (MQ=255)
              cgcgCAATAACCACTTCTTCACTGGCGTTGCCATTATTTCTTCCTTAGCTTTGCGCGACTTTACGGGGGATAAAATCATTCGCCACCGCTTCGCCTTGCGGATTCAGCGGCTGCGGCTGGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAACTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGCCAAGCGCGACATATTCGTTGATTCgcgc                                                                                                                                                                                                                               >  1:170825/1‑251 (MQ=255)
              cgcgCAATAACCACTTCTTCACTGGCGTTGCCATTATTTCTTCCTTAGCTTTGCGCGACTTTACGGGGGATAAAATCATTCGCCACCGCTTCGCCTTGCGGATTCAGCGGCTGCGGCTGGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAACTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGCCAAGCGCGACATATTCGTTGATTCgcgc                                                                                                                                                                                                                               >  1:49887/1‑251 (MQ=255)
                            ttcttcACTGGCGTTGCCATTATTTCTTCCTTAGCTTTGCGCGACTTTACGGGGGATAAAATCATTCGCCACCGCTTCGCCTTGCGGATTCAGCGGCTGCGGCTGGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAACTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGCCAAGCGCGAAATATTCGTTGATTCGCGCAGCGCCTGTAGGAc                                                                                                                                                                                                                 >  2:95683/1‑251 (MQ=255)
                                 cACTGGCGTTGCCATTATTTCTTCCTTAGCTTTGCGCGACTTTACGGGGGATAAAATCATTCGCCACCGCTTCGCCTTGCGGATTCAGCGGCTGCGGCTGGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAACTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGCCAAGCGCGACATATTCGTTGATTCGCGCAGCGACCGTAGGACCATc                                                                                                                                                                                                             <  1:137291/250‑1 (MQ=255)
                                                           tAGCTTTGCGCGACTTTACGGGGGATAAAATCATTCGCCACCGCTTCGCCTTGCGGATTCAGCGGCTGCGGCTGGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAACTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGCCAAGCGCGACATATTCGTTGATTCGCGCAGCGACCGTAGGACCATCGCCCACCAGCGCCGTCCCGGCAccgcc                                                                                                                                                                                  >  1:23858/1‑251 (MQ=255)
                                                                 tGCGCGACTTTACGGGGGATAAAATCATTCGCCACCGCTTCGCCTTGCGGATTCAGCGGCTGCGGCTGGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAACTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGCCAAGCGCGACATATTCGTTGATTCGCGCAGCGACCGTAGGACCATCGCCCACCAGCACCGTCCCGGCACCGCCGCGCAc                                                                                                                                                                            >  2:58018/1‑251 (MQ=255)
                                                                        cTTTACGGGGGATAAAATCATTCGCCACCGCTTCGCCTTGCGGATTCAGCGGCTGCGGCTGGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAACTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGCCAAGCGCGACATATTCGTTGATTCGCGCAGCGACCGTAGGACCATCGCCCACCAGCGCCGTCCCGGCACCGCCGCGCACTAAGccc                                                                                                                                                                     >  2:158052/1‑250 (MQ=255)
                                                                            aCGGGGGATAAAATCATTCGCCACCGCTTCGCCTTGCGGATTCAGCGGCTGCGGCTGGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAACTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGCCAAGCGCGACAGATTCGTTGATTCGCGCAGCGACCGTAGGACCATCGCCCACCAGCGCCGTCCCGGCACCGCCGCGCACTAAGCCAAcg                                                                                                                                                                  <  2:249828/250‑1 (MQ=255)
                                                                                  gATAAAATCATTCGCCACCGCTTCGCCTTGCGGATTCAGCGGCTGCGGCTGGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAACTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGCCAAGCGCGACATATTCGTTGATTCGCGCAGCGACCGTAGGACCATCGCCCACCAGCGCCGTCCCGGCACCGCCGCGCACTAAGCCAACgcccgccc                                                                                                                                                           >  1:349868/1‑251 (MQ=255)
                                                                                    aaaaaTCATTCGCCACCGCTTCGCCTTGCGGATTCAGCGGCTGCGGCTGGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAACTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGCCAAGCGCGACATATTCGTTGATTCGCGCAGCGACCGTAGGACCATCGCCCACCAGCGCCGTCCCGGCCCCGCCGCGCACTAAGCCAACGCCCGCCCAt                                                                                                                                                         <  2:244592/250‑1 (MQ=255)
                                                                                     aaaaTCATTCGCCACCGCTTCGCCTTGCGGATTCAGCGGCTGCGGCTGGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAACTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGCCAAGCGCGACATATTCGTTGATTCGCGCAGCGACCGTAGGACCATCGCCCACCAGCGCCGTCCCGGCACCGCCGCGCACTAAGCCAACGCCCGCCCATa                                                                                                                                                        <  1:133356/251‑1 (MQ=255)
                                                                                                           acTTGCGGATTCAGCGGCTGCGGCTGGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAACTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGCTGCCAAACGCGACATATTCGTTGATTCGCGCAGAGACCGTAGGACCATAGCCCACCAGCGCCGGCCCGGAACAGCAGCGCACTAAGCCAACGCCCTCCCCTAAATTGGGGCTGATCTCCCGAtt                                                                                                                                  >  2:11715/2‑251 (MQ=255)
                                                                                                                  gATTCAGCGGCTGCGGATGGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAACTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGACCGAAAGCACAAAACTGTCGATGCCAAGCGCTACAAATTCGTTGATTCGCGCAGCGACCGTAGGACCATCGCCCACCAGCGCCGTCCCGGCACCGCCGCGCACTAAGCCAACGCCCGCCCATAAATTGGGGCTGATCTCCCGATTGTCgcgc                                                                                                                           >  2:290306/1‑251 (MQ=255)
                                                                                                                    ttCAGCGGCTGCGGCTGGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAACTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGCCAAGCGCGACATATTCGTTGATTCGCGCAGCGACCGTAGGACCATCGCCCACCAGCGCCGTCCCGGCACCGCCGCGCACTAAGCCAACGCCCGCCCATAAATTGGGGCTGATCTCCAGATTGTCGCGCtt                                                                                                                         >  1:118384/1‑251 (MQ=255)
                                                                                                                             tgcggctgGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAACTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGCCAAGCGCGACATATTCGTTGATTCGCGCAGCGACCGTAGGACCATCGCCCACCAGCGCCGTCCCGGCACCGCCGCGCACTAAGCCAACGCCCGCCCATAAATTGGGGCTGATCTCCAGATTGTCGCGCTTGCCGTTATg                                                                                                                >  1:332412/1‑251 (MQ=255)
                                                                                                                                ggctgGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAACTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGCCAAGCGCGACATATTCGTTGATTCGCGCAGCGACCGTAGGACCATCGCCCACCAGCGCCGTCCCGGCACCGCCGCGCACTAAGCCAACGCCCGCCCATAAATTGGGGCTGATCTCCAGATTGTCGCGCTTGCCGTTATGTa                                                                                                              >  2:211559/1‑250 (MQ=255)
                                                                                                                                 gctgGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAACTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGCCAAGCGCGACATATTCGTTGATTCGCGCAGCGACCGTAGGACCATCGCCCACCAGCGCCGTCCCGGCACCGCCGCGCACTAAGCCAACGCCCGCCCATAAATTGGGGCTGATCTCCAGATTGTCGCGCTTGCCGTTATGTAAc                                                                                                            <  2:202584/251‑1 (MQ=255)
                                                                                                                                           tgccGGGAGGGCGACATCCAGAAGCGGAAACAGCAACTCGCCAACCCGATACGCTTCTCCCAGAGGCGGAGAGCCCGAAAGCACAAAACGGTCGATGCCAAGCGCGCCATATTCGTTGATTCGCGCAGCGACCGTAGGCCCATCCCCCACCAGCGCCGCCCGGGCACCGCCGCGCACTAAGCCAACGCCCTCCCATAAATTGGGGCTGATCTCCAGATTGTCGCGCTTGCCGTTATGTAACTCCGCCAtt                                                                                                   <  2:368928/248‑1 (MQ=255)
                                                                                                                                                  aTGGCGACATCCAGAAGCGGGAACAGCAACTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGCCAAGCGCGACATATTCGTTGATTCGCGCAGCGACCGTAGGACCATCGCCCACCAGCGCCGTCCCGGCACCGCCGCGCACTAAGCCAACGCCCGCCCATAAATTGGGGCTGATCTCCAGATTGTCGCGCTTGCCGTTATGTAACGCCGCCATTCGCTGTTg                                                                                           <  1:371843/251‑1 (MQ=255)
                                                                                                                                                    ggCGACATCCAGAAGCGGGAACAGCAACTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGCCAAGCGCGACATATTCGTTGATTCGCGCAGCGACCGTAGGACCATCGCCCACCAGCGCCGTCCCGGCACCGCCGCGCACTAAGCCAACGCCCGCCCATAAATTGGGGCTGATCTCCAGATTGTCGCGCTTGCCGTTATGTAACGCCGCCATTCGCTGTTGcc                                                                                         <  1:291515/251‑1 (MQ=255)
                                                                                                                                                             cAGAAGCGGGAACAGCAACTCGCCAACCCGCTCCGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGCCAAGCGCGACCTATTCGTTGATTCGCGCAGCGACCGTAGGACCATCGCCCACCAGCGCCGTCCCGGCACCGCCGCGCACTAAGCCAACGCCCGCCCATAAATTGGGGCTGATCTCCAGATTGTCGCGCTTGCCGGTATGTAACGCCGCCATACGCTGTTGCCCtccgcaat                                                                                >  2:229679/1‑244 (MQ=255)
                                                                                                                                                                    gggAACAGCAACTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGCCAAGCGCGACATATTCGTTGATTCGCGCAGCGACCGTAGGACCATCGCCCACCAGCGCCGTCCCGGCACCGCCGCGCACTAAGCCAACGCCCGCCCATAAATTGGGGCTGATCTCCAGATTGTCGCGCTTGCCGTTATGTAACGCCGCCATTCGCTGTTGCCCTACGGAAtccgtccg                                                                         >  1:59857/1‑251 (MQ=255)
                                                                                                                                                                      gAACAGCAACTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGCCAAGCGCGACATATTCGTTGATTCGCGCAGCGACCGTAGGACCATCGCCCACCAGCGCCGTCCCCGCACCGCCGCGCACTAAGCCAACGCCCGCCCATAAATTGGGGCTGATCTCCAGATTGTCGCGCTTGCCGTTATGTAACGCCGCCATTCGCTGTTGCCCTACGGAATCCGTCCggg                                                                       >  2:271510/1‑251 (MQ=255)
                                                                                                                                                                               cTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGCCAAGCGCGACATATTCGTTGATTCGCGCAGCGACCGTAGGACCATCGCCCACCAGCGCCGTCCCGGCACCGCCGCGCACTAAGCCAACGCCCGCCCATAAATTGGGGCTGATCTCCAGATTGTCGCGCTTGCCGTTATGTAACGCCGCCATTCGCTGTTGCCCTACGGAATCCGTCCGGGCGAATGCgg                                                              >  2:227466/1‑251 (MQ=255)
                                                                                                                                                                                 cGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGCCAAGCGCGACATATTCGTTGATTCGCGCAGCGACCGTAGGACCATCGCCCACCAGCGCCGTCCCGGCACCGCCGCGCACTAAGCCAACGCCCGCCCATAAATTGGGGCTGATCTCCAGATTGTCGCGCTTGCCGTTATGTAACGCCGCCATTCGCTGTTGCCCTACGGAATCCGTCCGGGCGAATGCgg                                                              <  1:49897/249‑1 (MQ=255)
                                                                                                                                                                                         cGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGCCAAGCGCGACATATTCGTTGATTCGCGCAGCGACCGTAGGACCATCGCCCACCAGCGCCGTCCCGGCACCGCCGCGCACTAAGCCAACGCCCGCCCATAAATTGGGGCTGATCTCCAGATTGTCGCGCTTGCCGTTATGTAACGCCGCCATTCGCTGTTGCCCTACGGAATCCGTCCGGGCGAATGCGGCCTGTGCttt                                                    >  1:136327/1‑251 (MQ=255)
                                                                                                                                                                                                     tcCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGCCAAGCGCGACATATTCGTTGATTCGCGCAGCGACCGTAGGACCATCGCCCACCAGCGCCGTCCCGGCACCGCCGCGCACTAAGCCAACGCCCGCCCATAAATTGGGGCTGATCTCCAGATTGTCGCGCTTGCCGTTATGTAACGCCGCCATTCGCTGTTGCCCTACGGAATCCGTCCGGGCGAATGCGGCCTGTGCTTTGGCGATAGTTtc                                        >  1:73475/1‑251 (MQ=255)
                                                                                                                                                                                                      ccAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGCCAAGCGCGACATATTCGTTGATTCGCGCAGCGACCGTAGGACCATCGCCCACCAGCGCCGTCCCGGCACCGCCGCGCACTAAGCCAACGCCCGCCCATAAATTGGGGCTGATCTCCAGATTGTCGCGCTTGCCGTTATGTAACGCCGCCATTCGCTGTTGCCCTACGGAATCCGTCCGGGCGAATGCGGCCTGTGCTTTGGCGATAGTTtca                                       <  2:68270/251‑1 (MQ=255)
                                                                                                                                                                                                       cAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGCCAAGCGCGACATATTCGTTGATTCGCGCAGCGACCGTAGGACCATCGCCCACCAGCGCCGTCCCGGCACCGCCGCGCACTAAGCCAACGCCCGCCCATAAATTGGGGCTGATCTCCAGATTGTCGCGCTTGCCGTTATGTAACGCCGCCATTCGCTGTTGCCCTACGGAATCCGTCCGGGCGAATGCGGCCTGTGCTTTGGCGATAGTTtcat                                      <  2:170825/251‑1 (MQ=255)
                                                                                                                                                                                                         gATGCGGATAGCCCGAAAGCACAACCCTGTCGATGCCAAGCGCGACATATTCGTTGATTCGCGCAGCGACCGTAGGACCATCGCCCACCAGCGCCGTCCCGGCCCCGCCGCGCACTAAGCCAACGCCCGCCCATAAATTGGGGCTGATCTCCAGATTGTCGCGCTTGCTGTTATGTAACGCCGCCATTCGCTGTTGCCCTACGGCCTCCGTCCGGGCGAATGCGGCCTGTGCTTTGGCGATAGTTtcatcc                                    >  2:123307/1‑250 (MQ=255)
                                                                                                                                                                                                                 aaGCCCGAAAGACCAAAACTGTCGATGCCAAGCGCGACATATGCGATGATTCGCGCAGCGACCGGAGGACCATCGCCCACCAGCGCCGTCCCGGCACCGCCGCGCACTAAGCCAACGCCCGCCCATAAATTGGGGCTGATCTCCAGATTGTCGCGCTTGCCGTTATGTAACGCCGCCATTCGCTGTTGCCCTACGGAATCCGTCCGGGCGAATGCGGCCTGTGCTTTGGCGATAGTTTCATCATCAAGATg                            <  2:98556/250‑1 (MQ=255)
                                                                                                                                                                                                                      cGAAAGCACAAAACTGTCGATGCCAAGCGCGACATATTCGTTGATTCGCGCAGCGACCGTAGGACCATCGCCCACCAGCGCCGTCCCGGCACCGCCGCGCACTAAGCCAACGCCCGCCCATAAATTGGGGCTGATCTCCAGATTGTCGCGCTTGCCGTTATGTAACGCCGCCATTCGCTGTTGCCCTACGGAATCCGTCCGGGCGAATGCGGCCTGTGCTTTGGCGATAGTTTCATCATCAAGATGCg                          >  2:75223/1‑248 (MQ=255)
                                                                                                                                                                                                                                aaaCTGTCGATGCCAAGCGCGACATATTCGTTGATTCGCGCAGCGACCGTAGGACCATCGCCCACCAGCGCCGTCCCGGCACCGCCGCGCACTAAGCCAACGCCCGCCCATAAATTGGGGCTGATCTCCAGATTGTCGCGCTTGCCGTTATGTAACGCCGCCATTCGCTGTTGCCCTACGGAATCCGTCCGGGCGAATGCGGCCTGTGCTTTGGCGATAGTTTCATCATCAAGATGCGAGATTAACCGCTc             >  2:328375/1‑251 (MQ=255)
                                                                                                                                                                                                                                 aaCTGTCGATGCCAAGCGCGACATATTCGTTGATTCGCGCAGCGACCGTAGGACCATCGCCCACCAGCGCCGTCCCGGCACCGCCGCGCACTAAGCCAACGCCCGCCCATAAATTGGGGCTGATCTCCAGATTGTCGCGCTTGCCGTTATGTAACGCCGCCATTCGCTGTTGCCCTACGGAATCCGTCCGGGCGAATGCGGCCTGTGCTTTGGCGATAGTTTCATCATCAAGATGCGAGATTAACCGCTcg            >  2:115808/1‑251 (MQ=255)
                                                                                                                                                                                                                                       ggATGCCAAGCGCGACATATTCGTTGATTCGCGCAGCGACCGTAGGACCATCGCCCACCAGCGCCGTCCCGGCACCGCCGCGCACTAAGCCAACGCCCGCCCATAAATTGGGGCTGATCTCCAGATTGTCGCGCTTGCCGTTATGTAACGCCGCCATTCGCTGTTGCCCTACGGAATCCGTCCGGGCGAATGCGGCCTGTGCTTTGGCGATAGTTTCATCATCAAGATGCGAGATTAACCGCTCGGCGGcc      <  2:336935/250‑1 (MQ=255)
                                                                                                                                                                                                                                          tGCCAAGCGCGACATATTCGTTGATTCGCGCAGCGACCGTAGGACCATCGCCCACCAGCGCCGTCCCGGCACCGCCGCGCACTAAGCCAACGCCCGCCCATAAATTGGGGCTGATCTCCAGATTGTCGCGCTTGCCGTTATGTAACGCCGCCATTCGCTGTTGCCCTACGGAATCCGTCCGGGCGAATGCGGCCTGTGCTTTGGCGATAGTTTCATCATCAAGATGCGAGATTAACCGCTCGGCGGCCTGc   >  1:286976/1‑251 (MQ=255)
                                                                                                                                                                                                                                           gcaaaGCGCGACATATTCGTTGATTCGCGCAGCGACCGTAGGACCATCGCCCACCAGCGCCGTCCCGGCACCGCCGCGCACTAAGCCAACGCCCGCCCATAAATTGGGGCTGATCTCCAGATTGTCGCGCTTGCCGTTATGTAACGCCGCCATTGGCTGGGGCCCTACGGAATCCGTCCGGGCGACTGCGGCCTGTGCTTTGGCGATAGTTTCATCATCAAGATGCGAGATTAACCGCTCGGCGGCCTGcc  >  2:67886/4‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                
AAACCAGGGGGCAACGCGCAATAACCACTTCTTCACTGGCGTTGCCATTATTTCTTCCTTAGCTTTGCGCGACTTTACGGGGGATAAAATCATTCGCCACCGCTTCGCCTTGCGGATTCAGCGGCTGCGGCTGGGGAATTTCCGGGATGGCGACATCCAGAAGCGGGAACAGCAACTCGCCAACCCGATACGCTTCTTCCAGATGCGGATAGCCCGAAAGCACAAAACTGTCGATGCCAAGCGCGGCATATTCGTTGATTCGCGCAGCGACCGTAGGACCATCGCCCACCAGCGCCGTCCCGGCACCGCCGCGCACTAAGCCAACGCCCGCCCATAAATTGGGGCTGATCTCCAGATTGTCGCGCTTGCCGTTATGTAACGCCGCCATTCGCTGTTGCCCTACGGAATCCGTCCGGGCGAATGCGGCCTGTGCTTTGGCGATAGTTTCATCATCAAGATGCGAGATTAACCGCTCGGCGGCCTGCC  >  NC_000913/994785‑995270

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑