Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 758,459 G→A 6.3% intergenic (+54/‑247) sdhB → / → sucA succinate dehydrogenase, FeS subunit/2‑oxoglutarate decarboxylase, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913758,4590GA 6.3% 205.8 / 3.5 63intergenic (+54/‑247)sdhB/sucAsuccinate dehydrogenase, FeS subunit/2‑oxoglutarate decarboxylase, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base G (24/35);  new base A (2/2);  total (26/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

CATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGTTTTTAAAGGTTCCTTCGCGAGCCACT  >  NC_000913/758315‑758605
                                                                                                                                                |                                                                                                                                                  
catGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATgcggc                                                                                                                                                <  3:30695/149‑1 (MQ=255)
  tGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGt                                                                                                                                              <  3:212534/149‑1 (MQ=255)
    aaCTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTg                                                                                                                                             <  2:427332/148‑1 (MQ=255)
      cTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGaa                                                                                                                                           <  1:126064/148‑1 (MQ=255)
      cTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAAc                                                                                                                                          <  4:246953/149‑1 (MQ=255)
       tGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAAc                                                                                                                                          <  1:901287/148‑1 (MQ=255)
             aGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTa                                                                                                                                    >  2:84998/1‑148 (MQ=255)
              gtgtATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                                  >  3:782073/1‑149 (MQ=255)
              gtgtATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                                  >  4:394793/1‑149 (MQ=255)
                gtATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                >  4:509374/1‑149 (MQ=255)
                   tGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                             <  4:25234/149‑1 (MQ=255)
                    gTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                             <  1:87204/148‑1 (MQ=255)
                    gTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCt                                                                                                                            >  3:93280/1‑149 (MQ=255)
                     tCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                                                                                           >  4:44254/1‑149 (MQ=255)
                        gAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGCATccggcactggt                                                                                                                                          <  3:441990/131‑7 (MQ=255)
                        gAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGCATccggcactggt                                                                                                                                          >  4:441990/1‑125 (MQ=255)
                          aGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAg                                                                                                                       <  1:614387/148‑1 (MQ=255)
                                tGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTAc                                                                                                                <  3:694012/149‑1 (MQ=255)
                                  aaCCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTAccc                                                                                                              <  4:782073/149‑1 (MQ=255)
                                      cGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTa                                                                                                           >  2:309896/1‑148 (MQ=255)
                                            gcgcCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGa                                                                                                    <  3:820566/149‑1 (MQ=255)
                                                    ggCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAagcgca                                                                                            <  3:509374/149‑1 (MQ=255)
                                                         tatCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATAAGGCGTTCACGCCGCATCCGGcactg                                                                                                                          >  3:819475/1‑110 (MQ=25)
                                                         tatCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATAAGGCGTTCACGCCGCATCCGGcactg                                                                                                                          <  4:819475/114‑5 (MQ=25)
                                                               gTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACTAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAgtc                                                                                 <  4:650199/149‑3 (MQ=255)
                                                               gTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGc                                                                                 <  3:855295/149‑1 (MQ=255)
                                                                tCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGc                                                                                 >  2:362188/1‑148 (MQ=255)
                                                                    tgttgtTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCAtt                                                                                                                  <  2:765914/111‑1 (MQ=255)
                                                                    tgttgtTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCAtt                                                                                                                  >  1:765914/1‑111 (MQ=255)
                                                                     gttgttGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAc                                                                            >  1:408040/1‑148 (MQ=255)
                                                                     gttgttGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAc                                                                            >  2:735664/1‑148 (MQ=255)
                                                                         ttgCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGa                                                                       <  4:93280/149‑1 (MQ=255)
                                                                          tgCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGaa                                                                      <  3:44254/149‑1 (MQ=255)
                                                                               cGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGTGTAACAAAGAAATGc                                                                  <  1:309896/148‑1 (MQ=255)
                                                                                  aaTGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAgg                                                               >  1:768788/1‑148 (MQ=255)
                                                                                      cGTAAACCGTAGGCCTGATATGACGCGCAAGCGTCGCATCAGGCAACCTGTGCCGGATGCGGCGTGATCGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAAt                                                           <  1:702745/148‑1 (MQ=255)
                                                                                       gTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATc                                                          >  1:833560/1‑148 (MQ=255)
                                                                                        tAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCCGGCGTAACAAAGAAATGCAGGAAATCtt                                                        >  4:697959/1‑149 (MQ=255)
                                                                                        tAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCtt                                                        >  3:237995/1‑149 (MQ=255)
                                                                                        tAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCtt                                                        >  4:507124/1‑149 (MQ=255)
                                                                                         aaaCCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCtt                                                        <  1:84998/148‑1 (MQ=255)
                                                                                             cGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTaaa                                                    <  1:446840/148‑1 (MQ=255)
                                                                                               tAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAAc                                                 <  3:804416/149‑1 (MQ=255)
                                                                                                      gATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGcccc                                           <  1:735664/148‑1 (MQ=255)
                                                                                                       aTAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCt                                          <  2:768788/148‑1 (MQ=255)
                                                                                                        tAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTg                                         <  2:408040/148‑1 (MQ=255)
                                                                                                        tAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTg                                         <  1:32240/148‑1 (MQ=255)
                                                                                                        tAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGa                                        >  4:319683/1‑149 (MQ=255)
                                                                                                           gACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGacac                                     <  4:237995/149‑1 (MQ=255)
                                                                                                           gACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGaca                                      <  1:901682/148‑1 (MQ=255)
                                                                                                                    gCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAg                             >  2:67453/1‑148 (MQ=255)
                                                                                                                     cGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGtt                           <  3:195027/149‑1 (MQ=255)
                                                                                                                       tCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACa                                                                           >  4:447274/1‑99 (MQ=255)
                                                                                                                       tCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACa                                                                           <  3:447274/99‑1 (MQ=255)
                                                                                                                            tCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGTTTTTAAAg                    >  4:430733/1‑149 (MQ=255)
                                                                                                                            tCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGTTTTTAAAg                    <  3:394793/149‑1 (MQ=255)
                                                                                                                              aGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGtt                           >  3:71591/1‑140 (MQ=255)
                                                                                                                              aGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGtt                           <  4:71591/140‑1 (MQ=255)
                                                                                                                              aGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGTTTTTAAAGGt                  >  3:337807/1‑149 (MQ=255)
                                                                                                                                gCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGTTTTTAAAGGtt                 >  1:667858/1‑148 (MQ=255)
                                                                                                                                gCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGTTTTTAAAGGtt                 <  1:292739/148‑1 (MQ=255)
                                                                                                                                gCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGTTTTTAAAGGTTc                >  4:703111/1‑149 (MQ=255)
                                                                                                                                    ccAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGTTTTTAAAGGTTCCtt             <  2:81945/148‑1 (MQ=255)
                                                                                                                                          ccGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGa                                                              >  2:149782/1‑93 (MQ=25)
                                                                                                                                          ccGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGa                                                              <  1:149782/93‑1 (MQ=25)
                                                                                                                                               tGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGTTTTTAAAGGTTCCTTCGCGAGCCACt  <  2:833560/148‑1 (MQ=255)
                                                                                                                                                gcggcgTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCt                                          >  4:103113/1‑107 (MQ=255)
                                                                                                                                                gcggcgTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCt                                          <  3:103113/107‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
CATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGTTTTTAAAGGTTCCTTCGCGAGCCACT  >  NC_000913/758315‑758605

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