Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,156,591 T→G 7.6% intergenic (+73/‑32) yqhD → / → dkgA aldehyde reductase, NADPH‑dependent/2,5‑diketo‑D‑gluconate reductase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,156,5910TG 7.6% 149.0 / 2.2 61intergenic (+73/‑32)yqhD/dkgAaldehyde reductase, NADPH‑dependent/2,5‑diketo‑D‑gluconate reductase A
Reads supporting (aligned to +/- strand):  ref base T (29/26);  new base G (2/3);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

CATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAG  >  NC_000913/3156446‑3156737
                                                                                                                                                 |                                                                                                                                                  
cATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCCCCGGTGCCTGcg                                                                                                                                                 >  1:95033/1‑148 (MQ=255)
 aTGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTc                                                                                                                                                <  2:136769/149‑1 (MQ=255)
 aTGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACAACTCACCGGAGCAGGCTc                                                                                                                                                <  2:43725/149‑1 (MQ=255)
     cccAACTGGGAGAAAATCATGACATTACGTTGGATGTCAGCAGCCGTATATACGAAGCCGCCCGCTAAGCTTATTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGt                                                                                                                                            <  1:905745/149‑1 (MQ=255)
        aaCTGGGCGAAAATGATGACATTTCGTTGGATGTCAGCCGCCGTATATTCGAAGCCGCTCGCTATGGTTTTTACGCCTCAATTATTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACATCAGCTCACCGGCGCCTGCTCCGGTGAg                                                                                                                                         <  2:930689/149‑1 (MQ=255)
        aaCTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAACCTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCAAAACACCTCACCGGAGCCTGCTCCGGTGAg                                                                                                                                         >  1:949970/1‑149 (MQ=255)
        aaCTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAg                                                                                                                                         <  1:901137/149‑1 (MQ=255)
         aCTGGGCGAAAATCATGACATTAGGTTGGATGGCAGCCGCCGTATATAGGAAGCCGCCCGCTAAGATTTTTACGCCTCAAACTTTCGTTTTCGGGGATTTCGTCCAGAGTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGt                                                                                                                                        <  1:726220/149‑1 (MQ=255)
            gggCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCa                                                                                                                                     <  1:892578/149‑1 (MQ=255)
              gCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTCCGCCTCAAACTTTCGTTTTCGGGCTTTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCata                                                                                                                                   >  2:63426/1‑149 (MQ=255)
                gAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCAtata                                                                                                                                 >  1:247514/1‑149 (MQ=255)
                 aaaaTCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATaa                                                                                                                                >  1:518427/1‑149 (MQ=255)
                   aaTCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAag                                                                                                                              <  1:301997/149‑1 (MQ=255)
                      cATGACATTACGTTCGATCTCAGCCGCCGTATATACGAAGCCGCCCGCAAAGCTTTTTACGCCTCAAACTTTCGTATTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGCTGAGTTCATATAAaggag                                                                                                                           >  2:866919/1‑149 (MQ=255)
                          aCATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTATGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAAc                                                                                                                       <  1:1007590/149‑1 (MQ=255)
                                cGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATgg                                                                                                                 <  1:925201/149‑1 (MQ=255)
                                    ggATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTaa                                                                                                             <  1:427557/149‑1 (MQ=255)
                                      aTGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCTAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGCGCCTGCTCCGGTGAGTTCATATTAAGGAGGAACGTATGGcgtttc                                                                                                           >  2:765670/1‑144 (MQ=255)
                                       tGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATcc                                                                                                          >  1:861206/1‑149 (MQ=255)
                                        gTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGGTAATccc                                                                                                         >  1:848094/1‑148 (MQ=255)
                                          cAGCCGCCGTATATAAGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAAATTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAAc                                                                                                       <  1:346482/149‑1 (MQ=255)
                                                        aCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTAc                                                                                         >  2:1004239/1‑149 (MQ=255)
                                                        aCGAAGCCGCACGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTAc                                                                                         <  1:299816/149‑1 (MQ=255)
                                                        aCGAAGACGCCCGCTAAGATTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTAc                                                                                         <  1:211200/149‑1 (MQ=255)
                                                          gAAGCCGCCCGCTAAGCTTTTTACGCCTCACACTTTCGTTGTCGGGCATTTCGTCCAGACTTAAGTTCCCAACACCTCACCGGAGCCTGCTCCGGTGCGTTCATCTAAAGGAGGAACGTATGGCTAATCCCACCGTTGATATGCTAccg                                                                                       >  1:147051/1‑147 (MQ=255)
                                                          gAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAAcc                                                                                                      >  1:1000728/1‑134 (MQ=255)
                                                          gAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAAcc                                                                                                      <  2:1000728/134‑1 (MQ=255)
                                                          gAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCTTTATTAAGCTACAg                                                                                       >  2:57453/1‑149 (MQ=255)
                                                               cgcccgcTAAGCTTTTTACGCATCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTTCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATgg                                                                                  <  1:422732/149‑1 (MQ=255)
                                                                 cccgcTAAGCTTTTTACGCCTCACACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCAGCTCCGGTGAGTTCATATAAAGGAGGATCGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCa                                                                                >  2:823817/1‑149 (MQ=255)
                                                                    gcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAAAGTTTGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATg                                                                             >  1:324740/1‑149 (MQ=255)
                                                                         gCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATTAAGGAGGAACGTATGGCTAATCCAACCGTGATTAAGCTACCGGATGGCAATGTCATg                                                                        >  2:174026/1‑149 (MQ=255)
                                                                          cTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGc                                                                       >  1:1009534/1‑149 (MQ=255)
                                                                          cTTTTTACGCCTAAAACTTTCGTTTTCGGGCATTTCTTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGc                                                                       <  2:924465/149‑1 (MQ=255)
                                                                           tttttACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGcc                                                                      >  1:279105/1‑149 (MQ=255)
                                                                                    cTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCGAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTggg                                                             >  2:934955/1‑149 (MQ=255)
                                                                                      cAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACACCACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc                                                           <  2:5908/149‑1 (MQ=255)
                                                                                            ttCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGTGCg                                                     >  2:875357/1‑149 (MQ=255)
                                                                                            ttCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAAAGTCATGCCCCAGCTGGGACTGGGCg                                                     >  1:780532/1‑149 (MQ=255)
                                                                                                ttttCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTg                                                 <  1:407086/149‑1 (MQ=255)
                                                                                                 tttCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATTAACGAGGACCGTATTGCTAATCCATCCGTTATTAAGCTCGAGGATGGCAATGTCATGCACCAGCTGGGAGTGGGCGTATgg                                                >  2:719480/1‑149 (MQ=255)
                                                                                                    cGGGCATTTCGTCCAGACTTAATTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgcaa                                             <  1:322324/149‑1 (MQ=255)
                                                                                                    cGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGGGAGTTCATATAAAGGAGGACAGTATGGCTAATCCAACCGTGATTACGCTACAGGATGGCAATGTCATGCCCCAGCTTGGACTGGGCGTCTGgcaa                                             >  2:956920/1‑149 (MQ=255)
                                                                                                        cATTTCGGCCAGACTTAAGTACACAACACCGCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGCAACGTCTGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgcaagcaa                                         <  2:129576/149‑1 (MQ=255)
                                                                                                         aTTTCGTCCAGACTTAAGTTCACAACACCTCACAGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGcaagcaag                                        <  2:900949/149‑1 (MQ=255)
                                                                                                           ttCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAACGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTa                                      <  1:787421/149‑1 (MQ=255)
                                                                                                                   gACTTAAGTTCACAACACCTCACCGGAGGCGGCTCCGGTGAGTTCTTATAAAGGAGGAACGTATGGAGAATCCAACCGGTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGaa                              <  2:95033/149‑1 (MQ=255)
                                                                                                                   gACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGaa                              <  1:272679/149‑1 (MQ=255)
                                                                                                                      ttAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATTGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTa                           <  1:900280/149‑1 (MQ=255)
                                                                                                                      ttAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGAAATCAAGTAATGAGGAAGTa                           >  2:786609/1‑149 (MQ=255)
                                                                                                                          gTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCa                       >  1:302955/1‑149 (MQ=255)
                                                                                                                           ttCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATATTGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGCCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAAAGAGGAAGTAATCAc                      >  1:13865/1‑149 (MQ=255)
                                                                                                                              acaacaCCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGc                   >  2:139983/1‑149 (MQ=255)
                                                                                                                              acaacaCCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGc                   <  2:815365/149‑1 (MQ=255)
                                                                                                                                aacaCCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTATGCTACAGGATGGCGATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGCGGAAGTAATCACCGCCa                 >  1:343125/1‑149 (MQ=255)
                                                                                                                                  cacCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGCTGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCAtt               <  2:861206/149‑1 (MQ=255)
                                                                                                                                     cTCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATTGCAATGTCATGCCACAGCTGGGACTGGTCGTCTGGAAATCAAGTAATGACGAAGTAATCACCGACATTCaa            >  2:159573/1‑149 (MQ=255)
                                                                                                                                     cTCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCAAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGGACTCTCCGCCAGTCaa            >  2:133135/1‑149 (MQ=255)
                                                                                                                                          cGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATTGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGAAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGc       <  1:875357/149‑1 (MQ=255)
                                                                                                                                          cGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGc       >  1:297507/1‑149 (MQ=255)
                                                                                                                                               ccTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCAATTAAAAAAGCGTTat  >  1:868785/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CATGACCCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAGCGTTAG  >  NC_000913/3156446‑3156737

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: ‑